function z= euler_putcha(a,b,y0,h,maxits) t(1)=a; w(1)=y0; k=a; while k<=b && k<=maxits w(k)=w(k-1)+h*fun(t,w(k-1)); t(k)=a+k*h; end z = [t;w]