function [w,t ] = euler(f,a,b,n,alpha ) w=zeros(n+1,1); t=linspace(a,b,n+1); h=(b-a)/n; w(1)=alpha; for i=1:n, w(i+1)=w(i)+h*f(t(i),w(i)); end end