help with euler code
6 views (last 30 days)
Show older comments
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 Comments
Answers (1)
Walter Roberson
on 9 Feb 2014
Edited: Walter Roberson
on 9 Feb 2014
After the plot() call, add
hold on
Change your sprintf() line to
fprintf('%.14f %.14f\n', t, y);
2 Comments
Walter Roberson
on 9 Feb 2014
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.
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!