[問題] 用ODE45解ODE加入FOR迴圈問題 已回收
解一ODE方常式如下:dy/dt=(1-y)/h+(1+y^2) 其中h是常數給一個定值隨便
code如下:
function dy=func_1(t,y);
dy=(1-y)/h+(1+y^2);
回command window:
[t,y]=ode45(@func_1,[0,10],[0])
plot(t,y)
以上可以跑出y對時間t的變化作圖
但若今天我想將求出的y再依序帶回function，方程式變的如下:
dy(i)/dt=(y(i-1)-y(i))/h+(1+y(i)^2); i=0:100;
我code該如何寫呢???題目也就是利用y(1)去求y(2),y(2)去求y(3)....以此類推
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.125.43.59
推
11/23 11:52, , 1F
11/23 11:52, 1F
→
11/23 16:24, , 2F
11/23 16:24, 2F