MATLAB code for stiff differential equation with explicit Euler method

조회 수: 16 (최근 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
Shumaila Khan
Shumaila Khan 2017년 4월 27일
Actually, I am working on autonomous case
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 CenterFile Exchange에서 Linear Algebra에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by