clear t						%clear old time steps and
clear y						%y values from previous runs
a=0;						%Initial time
b=1;						%Final time
N=10;						%Number of time steps
y0=0; 						%Initial value y(a)
h=(b-a)/N; 					%Time step
t(1)=a;
y(1)=y0; 

for n=1:N					%This loop carries out Euler's
	t(n+1)=t(n)+h;				%method. Note that the function
	y(n+1)=y(n)+h*f(t(n),y(n));		%f(t,y)=dy/dt is called
end

plot(t,y)
title(['Euler Method using N=', num2str(N), ' steps, by MYNAME'])

