Sir Model using Forward Euler
18 次查看(过去 30 天)
显示 更早的评论
% Solve the system x'=y, y'=-x using Forward Euler
clear all;
clc;
h=0.001;
T=10.0;
xs=0.2;
x(1)=xs;
ys=0.3;
y(1)=ys;
ts=0.;
t(1)=ts;
index=1;
while ts<=T
ts
%Forward Euler Step
ts=ts+h;
xshold=xs+h*ys;
ys=ys-h*xs;
xs=xshold;
index=index+1;
x(index)=xs;
y(index)=ys;
t(index)=ts;
end
figure(1)
plot(t,x); hold on;
plot(t,y,'g')
figure(2)
plot(x,y)
Here is the code for Forward Euler that my professor has provided us. The photo below are the given equations for the problem. I need to solve the ODE using forward euler. I am struggling to comprehend even where to start for putting the equations in.
1 个评论
回答(1 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!