what is wrong with this euler method code?

5 vues (au cours des 30 derniers jours)
Louis
Louis le 18 Mar 2014
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?

Réponses (1)

Sagar Damle
Sagar Damle le 18 Mar 2014
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));

Catégories

En savoir plus sur Mathematics dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by