请分析下doc ode45 中示例程序产生错误的问题

chinamiracle

function ha
clear all,clc
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),'.')

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文件里面却出错这是为什么？我觉得它应该可以的啊

ch_j1985

原帖由 chinamiracle 于 2008-5-13 03:01 发表

                               
登录/注册后可看大图

function ha
clear all,clc
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),'.')

function dy = ri ...

在我电脑上运行没有报错，能说明一下报的什么错吗？

eight

原帖由 chinamiracle 于 2008-5-13 03:01 发表

                               
登录/注册后可看大图

function ha
clear all,clc
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),'.')

function dy = ri ...

发帖前请多看看置顶帖，清晰地描述你遇到的问题

chinamiracle

哦呵呵，是我的问题，昨天晚上傻了哈一下没注意 ，函数名没一致