MATLAB code for stiff differential equation with explicit Euler method

15 visualizaciones (últimos 30 días)
Shumaila Khan
Shumaila Khan el 27 de Abr. de 2017
Comentada: Torsten el 27 de Abr. de 2017
I wrote a following program for of 2 differential equations for stiff differential equation with explicit Euler method. But getting the following error. Any idea? because I don't have, I am new at MATLAB.

Respuestas (1)

Torsten
Torsten el 27 de Abr. de 2017
y(n+1) = y(n)+h*f([x(n) y(n)]);
x(n+1) = x(n)+h*g([x(n) y(n)]);
if your ODE system reads
y' = 98x+198y
x' = -99x-199y
Best wishes
Torsten.
  2 comentarios
Shumaila Khan
Shumaila Khan el 27 de Abr. de 2017
Actually, I am working on autonomous case
Torsten
Torsten el 27 de Abr. de 2017
... and give x and y an initial value x(1) and y(1) before you enter the for-loop:
f=@(y)...;
g=@(x)...;
h=10;
x(1)=2;
y(1)=3;
for n=1:10
v=[x(n) y(n)];
y(n+1)=y(n)+h*f(v);
x(n+1)=x(n)+h*g(v);
end
Best wishes
Torsten.

Iniciar sesión para comentar.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by