The forward Euler scheme for an ordinary differential equation
1 次查看(过去 30 天)
显示 更早的评论
I'm having a lot of trouble with a homework problem that requires me to plot a forward euler scheme. I feel like I'm either close to getting somewhere, or, I've gone about solving this problem entirely wrong.
The problem is:
The code I have thus far is:
clear variables
clc
n=100;
h=0.1;
k=0;
x(k+1)=0;
y(k+1)=0;
for k=1
x(k+1)=x*k+h;
y(k+1)=y*k+h*(9.8-0.196*y1);
end
for k=2
x(k+1)=x*k+h;
y(k+1)=y*k+h*(9.8-0.196*y2);
end
for k=3
x(k+1)=x*k+h;
y(k+1)=y*k+h*(9.8-0.196*y3);
end
for k=4
x(k+1)=x*k+h;
y(k+1)=y*k+h*(9.8-0.196*y4);
end
plot(x,y)
I appreciate any pointers that anyone may have to help get me moving in the right direction.
0 个评论
采纳的回答
Alan Stevens
2021-4-11
More like this
n=100;
h=0.1;
k=0;
x(1)=0;
y(1)=48;
for k=1:n
x(k+1)=x(k)+h;
y(k+1)=y(k)+h*(9.8-0.196*y(k));
end
plot(x,y)
0 个评论
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!