Здравствуйте,
Мне нужно решить уравнение теплопроводности  du/dt=-d(du/dx)/dx+exp(u)-sin(x)  с  граничными условиями u(0,t)=u(Pi,t)=1. Начальный профиль беру u(x,0)=1.  Меня интересует слой, когда u перестанет зависеть от времени.
Я написал следующую программу: 
 
Сетка равномерная, шаг времени - t, координаты - h. u0-старое поколение, u1 - новое.
while нужен, чтобы прекратить работу, когда перестанет меняться u.
Проблема в том, что результат работы этой программы расходится.
Я не смог найти ошибку, по-этому прошу вашей помощи.

Это сообщение отредактировал(а) plunger - 30.1.2012, 11:21

Не знаю, есть ли ошыбки в программе -- но когда я увидел du/dt = ... exp(u) -- я сразу подумал, а не расходится ли? А тем более с начальными условиями u=1 (т.е. exp(u)=2.71), и компенсацыей этой экспоненты через sin(x), т.е. значение меньшэ 1. Тут уж будет удивительно, если это уравнение не разойдётся.

Спасибо, проблема оказалась в методе)

Странное это какое-то уравнение теплопроводности - с отрицательным коэффициентом и экспоненциальным источниковым членом... Метод - методом, но Вы точно не перепутали само уравнение?