Подскажите, в чем ошибка (дифф.ур)

Форум пользователей пакета Mathematica

Модератор: Admin

asisai
Сообщения: 4
Зарегистрирован: Вт апр 08, 2014 6:16 pm

Подскажите, в чем ошибка (дифф.ур)

Сообщение asisai » Вт апр 08, 2014 8:14 pm

Здравствуйте!
Пожалуйста, помогите найти ошибку в решении уравнения с помощью ф-и NDSolve в Mathematica


Sol = NDSolve[{D[f[x]*Cm[x, t], t] == D[Dch4*f[x]^2*D[Cm[x, t], x], x] - D[f[x]*u*Cm[x, t], x], Cm[0, t] == 0, Сm[x, 0] == 0, Cm[x, 600] == 30}, Cm, {t, 0, 2000}, {x, 0, 600}]

Plot[Evaluate[Cm[x, t] /. sol, {x, 0, 600}]]

Markiyan Hirnyk
Сообщения: 1181
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Re: Подскажите, в чем ошибка (дифф.ур)

Сообщение Markiyan Hirnyk » Ср апр 09, 2014 9:18 am

asisai писал(а):Здравствуйте!
Пожалуйста, помогите найти ошибку в решении уравнения с помощью ф-и NDSolve в Mathematica


Sol = NDSolve[{D[f[x]*Cm[x, t], t] == D[Dch4*f[x]^2*D[Cm[x, t], x], x] - D[f[x]*u*Cm[x, t], x], Cm[0, t] == 0, Сm[x, 0] == 0, Cm[x, 600] == 30}, Cm, {t, 0, 2000}, {x, 0, 600}]

Plot[Evaluate[Cm[x, t] /. sol, {x, 0, 600}]]

NDSolve::ndode: Input is not an ordinary differential equation. В частности, не заданы f[x], u, Dch4.

asisai
Сообщения: 4
Зарегистрирован: Вт апр 08, 2014 6:16 pm

Сообщение asisai » Ср апр 09, 2014 10:44 am

Да, все верно, в сообщении не указала параметры и переменные. Но они у меня заданы:
u = 20;
Dch4 = 314;
f=ff+(f0-ff)*exp(-fa*x);
f0 =0.95;
ff =0.7;
fa =0.0013;

Ошибка остается.(

Markiyan Hirnyk
Сообщения: 1181
Зарегистрирован: Вс дек 04, 2011 11:07 pm

Частные производные

Сообщение Markiyan Hirnyk » Чт апр 10, 2014 1:49 pm

Ваше уравнение вроде бы в частных производных, не так ли?

asisai
Сообщения: 4
Зарегистрирован: Вт апр 08, 2014 6:16 pm

Re: Частные производные

Сообщение asisai » Ср май 14, 2014 2:38 pm

Да в частных, метод решения Lines...
Теперь нужно решить систему и вероятно в том виде, в каком она у меня есть, не решить. Вылезает ошибка Maximum number of 10000 steps reached at the point t == 0.0012830624959653718`. при tmax=2000. Что делать, не знаю(((