MATLAB code for stiff differential equation with explicit Euler method

14 vues (au cours des 30 derniers jours)
Shumaila Khan
Shumaila Khan le 27 Avr 2017
Commenté : Torsten le 27 Avr 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.

Réponses (1)

Torsten
Torsten le 27 Avr 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 commentaires
Shumaila Khan
Shumaila Khan le 27 Avr 2017
Actually, I am working on autonomous case
Torsten
Torsten le 27 Avr 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.

Connectez-vous pour commenter.

Community Treasure Hunt

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

Start Hunting!

Translated by