问题描述:
下面程序size函数的错误(matlab)
function [tvals,yvals] = code2(t0,tf,dt,y0)
tvals=(t0:dt:tf);
n=tf/dt;
[t,y]=size(n);
for i=1:n
y(i+1)=y(i)+2*dt*y(i)/t(i)+2*t(i)^3*dt-3*t(i)^2*dt;
end
hold on
yvals=tvals.^4-3*tvals.^3+tvals.^2*(t0^2-3*t0-y0/(t0^2));
legend('Numerical solution','Analytical solution');
title('IVP of question 2')
xlabel('t')
ylabel('y')
function [tvals,yvals] = code2(t0,tf,dt,y0)
tvals=(t0:dt:tf);
n=tf/dt;
[t,y]=size(n);
for i=1:n
y(i+1)=y(i)+2*dt*y(i)/t(i)+2*t(i)^3*dt-3*t(i)^2*dt;
end
hold on
yvals=tvals.^4-3*tvals.^3+tvals.^2*(t0^2-3*t0-y0/(t0^2));
legend('Numerical solution','Analytical solution');
title('IVP of question 2')
xlabel('t')
ylabel('y')
问题解答:
我来补答展开全文阅读