what is wrong with this euler method code?

x=0;

y=0;

a=0;

b=5;

h=.01;

n=(b-a)/h;

for i=1:n

x(i)=a+(i-1)*h;

x(i+1)=a+i*h;

f = y(i)-5-exp(x(i)/2)sin*(5*x(i))+5*exp(x(i)/2)cos*(5*x(i));

y(i+1)=y(i)+h*f;

end

plot(x,y)

hold on

BTW how do I find exact solution and percent error?

Answer 1

[Sagar Damle]

I think there is an error in line

f = y(i)-5-.........

Try this line instead-

f = y(i)-5*exp(x(i)/2)*sin(5*x(i))+5*exp(x(i)/2)*cos(5*x(i));