OTVET
Newbie | Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору
Помогите, пожалуйста, решить систему дифференциальных уравнений с использованием функции numol. При помощи Pdesolve не соблюдается условие электронейтральности (алгебраическое уравнение). Cp(x,t),Ck(x,t),Ca(x,t) – искомые функции. В numol алгебраическое уравнение должно восприниматься вместе с системой дифференциальных уравнений. Кака это правильно записать?
zp = zk = -za = 1
Dp = Dk = Da = 10 -5
Алгебраическое уравнение:
zpCp(x,t) + zkCk(x,t) + zaCa(x,t) = 0
Система дифференциальных уравнений:
Cpt(x,t) = DpCpxx(x,t) + zpDpCpx(x,t)
Ckt(x,t) = DkCkxx(x,t) + zkDkCkx(x,t)
Cat(x,t) = DaCaxx(x,t) + zaDaCax(x,t)
Условия:
Cp(0,t) = 1 Ck(0,t) = 0
Cp(x,0) = 0 Ck(x,0) = 1 Ca(x,0) = 1
Cpx(L,t) = 0 Ckx(L,t) = 0 Cax(L,t) = 0
|
, в том же програмировании индексы вообще не используются и в перемнных никто еще не запутался