[fgabertb]

作者: 闫世军 neu

函数fun1.m文件为
function dy=fun1(t,y)
u=[t/1000+0.1 t/1000+0.1];
dy=zeros(6,1);
m1=0.2;m2=0.000034;m3=0.0337;m4=0.0386;
w1=0.3;w2=0;w3=0.1;
p=0;
k1...

dy(2) 算出來應該是一個值

但是因為 u(1:366) 這個地方會使得 dy(2) 算出來是一個向量 所以不對了

這邊修一下吧:smile: