大家好
小弟我寫了一些程式，但是現在有個問題
我大致上利用ode45給的範例
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);
我想問說
假設我在function那邊增加一個變數a
a在做微分方程的時候可以求得
但我想要在執行完畢的時候抓出a這個變數的值
有人知道怎麼用嗎
目前我是想說將第一行改成
function [dy,a] = rigid(t,y)
但是ode45那邊我試著改成[T,Y,a] = ode45(省略)
好像會有問題
請各位大大幫個忙，謝謝(有人懂我的問題嗎~)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.39.44