help with euler code
8 次查看(过去 30 天)
显示 更早的评论
i wrote the following matlab( for a euler) code and i can't get it to plot, or provide a figure w. the results. can anyone help?
function [dy] = y(n,to,t,y0)
h=[0.1]
t0=0
y0=1
n=50
h=0.1
t=t0
y=y0
for i=1:n
k=2-exp(-4*t)-2*y
y=y+h*k
t=t+h
plot (t,y)
sprintf('t',' y')
end
0 个评论
回答(1 个)
Walter Roberson
2014-2-9
编辑:Walter Roberson
2014-2-9
After the plot() call, add
hold on
Change your sprintf() line to
fprintf('%.14f %.14f\n', t, y);
2 个评论
Walter Roberson
2014-2-9
If you add the "hold on" after the plot() call then the plot is not blank, but the points are not very visible. Try using
plot(t, y, 'r*')
to make the points more visible.
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!