Обозначим
всю правую часть первого уравнения f1(x1, x2), а
всю правую часть второго уравнения f2(x1, x2). Пишем соответствующие функции.
Задаем константы (значения приведены для примера):
Код:
dt=0.01; // шаг по t
tmax=1.0; // предельное значение t
Заводим переменные с начальными значениями:
Код:
t:=0; x1:=0; x2:=0;
Дальше в цикле делаем следующее:
Код:
while t<=tmax do
begin
вывод t, x1, x2
x1:=x1+f1(x1, x2)*dt;
x2:=x2+f2(x1, x2)*dt;
t:=t+dt;
end
Это называется
метод Эйлера.