na-lab-eulermodified

code

clc;

f=@(t,y)(-y+2*cos(t));

t0=input("enter the initial value");

b=input("enter at which value is to be evaluated");

h=input("input step size");

y0=input("the value of function at t0");

n=(b-t0)/h;

for i=1:n

k1=h*f(t0,y0);

k2=h*f(t0+h,y0+k1);

y1=y0+(k1+k2)/2;

t0=t0+h;

y0=y1;

end

disp(y1);