MATLAB code for stiff differential equation with explicit Euler method

14 Ansichten (letzte 30 Tage)
Shumaila Khan
Shumaila Khan am 27 Apr. 2017
Kommentiert: Torsten am 27 Apr. 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.

Antworten (1)

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

Melden Sie sich an, um zu kommentieren.

Community Treasure Hunt

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

Start Hunting!

Translated by