請教各位大神
Q: 問題描述：傘兵問題，（1）求第十秒的速度，（2）依time stop=0.5 和 time stop=
1 ,繪出0到16秒到時間vs速度圖
以下是小弟的程式，不知道那裡錯，拜託幫忙指導
function value = veiocity3(dt,ti,tf,tv,vi,m,cd)
to(1)=0;vo(1)=0;
g=9.81; t=ti; v=vi;
n=(tf-ti)/dt;
for i=2:n+1
dvdt=deriv(v,m,cd);
t=t+dt;
to(i,1)=t;
vo(i,1)=v;
if t ==tv, value=v; end
end
vexact=sqrt(g*m/cd)*tanh(sqrt(g*cd/m)*to);
plot(to,vo,’b*-‘,to(:,1),vexact,’r-‘)
title(‘comparision of exact and numerical solution.’)
xlabel(‘values of time’)
ylabel(‘values of velocity’)
legend(strcat(‘dt=‘,num2str(dt)),’exact solution’)
function dv=deriv(v,m,cd)
g=9.81;
dv=g-(cd/m)*v^2;