大家好 小弟我写了一些程式，但是现在有个问题 我大致上利用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