what is wrong with this euler method code?

1 次查看（过去 30 天）
Louis
Louis 2014-3-18
回答： Sagar Damle 2014-3-18
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?

请先登录，再回答此问题。

回答（1 个）

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

类别

在 Help Center 和 File Exchange 中查找有关 Mathematics 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by