The forward Euler scheme for an ordinary differential equation

Question

Brendan Clark 2021-4-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/798637-the-forward-euler-scheme-for-an-ordinary-differential-equation

回答： Alan Stevens 2021-4-11

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)