Nana
North Carolina A&T State University
自 2011 起处于活动状态
Followers: 0 Following: 0
I copied this
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),'.')
and tried to run but it gave error message "??? Input argument "y" is undefined.
Error in ==> rigid at 3
dy(1) = y(2) * y(3);
What should be the problem?
Feeds
提问
Example in your ode45 tutorial is not running
I have the same problem trying to solve 3 nonlinear differential equations for my modeling. I tried to run your own tutorial ex...
13 years 前 | 2 个回答 | 0
