求各位高手修改修改此程序，不知道是那里出问题了

熊猫酒仙

clear;clc;
function dy=dfun(t,y)
dy=zeros(4,1);
r=0.1+0.01*cos(2*t);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(9.8*sin(y(1))+r*y(4).^2*sin(y(1))*cos(y(1))+0.04*y(2)*sin(2*t))/r;
dy(4)=(0.04*y(4)*sin(2*t)*sin(y(1))-2*r*y(1)*y(4)*cos(y(1)))/(r*sin(y(1)));
Y0=[4*pi/3,0,0,0.01];
[T,Y]=ode45('dfun',[16 21],Y0);
plot(T,Y)



不知道哪里出了问题，程序运行不出来，求高人相助。

happy

脚本文件中不能定义函数
可以将函数独立写成函数文件，脚本文件调用
或者把上述脚本文件改写成函数文件

liuyong

好像function的变量不行吧