Euler's method

Answers (2)

Steve Areola
Steve Areola on 5 Aug 2023

0 votes

clc
a=0;b=2;N=10; alph=0.5;
h=(b-a)/N; t(1)=a; w(1)=alph;
f = @(t,y) y-t^2+1;
for i=2:N+1
t(i)=a+(i+1)*h;
w(i)=w(i-1)+h*f(t(i-1),w(i-1));
end
disp("t w")
x=[t;w];
fprintf("%5.9f %5.9f\n",x)

Categories

Find more on Programming in Help Center and File Exchange

Tags

Asked:

on 13 Jan 2021

Answered:

on 5 Aug 2023

Community Treasure Hunt

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

Start Hunting!