作者samm3320 (sam)
看板MATLAB
標題[討論] 解微分方程
時間Fri Dec 16 14:59:45 2011
各位好最近需要用Matlab解微分方程
我依照
http://www.mathworks.com/help/techdoc/ref/ode23.html
網站上的範例
寫入以下code
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
但是matlab卻告訴我
??? Input argument "y" is undefined.
Error in ==> ode1 at 3
dy(1) = y(2) * y(3);
這到底是怎麼回事呢@@
範例都沒辦法解出來
一直卡在這裡沒辦法解我要的方程式
感謝大家的幫忙
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.114.14.86
1F：→ ejialan:你必須把前五行存成rigid.m 再另外執行後三行 12/16 15:58
2F：→ samm3320:成功了!感謝~ 12/16 18:18