MATLAB code for stiff differential equation with explicit Euler method

14 次查看（过去 30 天）
Shumaila Khan
Shumaila Khan 2017-4-27
评论： Torsten 2017-4-27
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.

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2017-4-27
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 个评论
Torsten
Torsten 2017-4-27
... 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.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Linear Algebra 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by