请分析下doc ode45 中示例程序产生错误的问题
function haclear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')
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);
我发现好怪了,一下想不通。这是ode45帮助文档里的例子程序,我把它分开存放运行可以做出图形,而同时存放 在一个m文件里面却出错这是为什么?我觉得它应该可以的啊 原帖由 chinamiracle 于 2008-5-13 03:01 发表 http://www.chinavib.com/forum/images/common/back.gif
function ha
clear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')
function dy = ri ...
在我电脑上运行没有报错,能说明一下报的什么错吗? 原帖由 chinamiracle 于 2008-5-13 03:01 发表 http://www.chinavib.com/forum/images/common/back.gif
function ha
clear all,clc
options = odeset('RelTol',1e-4,'AbsTol',);
= ode45(@rigid,,,options);
plot(t,y(:,1),'-',t,y(:,2),'-.',t,y(:,3),'.')
function dy = ri ... 发帖前请多看看置顶帖,清晰地描述你遇到的问题 哦呵呵,是我的问题,昨天晚上傻了哈一下没注意 ,函数名没一致
页:
[1]