[交流] MATLAB程序总是出错? 已有1人参与

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%function dx=myfun(t,x)
syms a e f g a1 a2 e1 e2 f1 f2 g1 g2
dx = zeros(5,1);

dx(1)=x(3)-x(2)+x(4)*x(2)-2*x(5)*x(2)-x(2)*x(2)-x(5)-e-f-g;
dx(2)=x(1)-x(3)+x(4)*x(2)-x(2)*x(2)-x(5)*x(5)-x(5)*x(2)-e*f-g;
dx(3)=x(2)-x(1)+x(5)*x(1)-x(4)*x(3)+x(3)+x(4)+e*g+f;
dx(4)=x(5)*a+x(3)*x(1)+x(3)+e*f*g*x(6);
dx(5)=x(4)*a+x(2)*x(1)-e*f*x(1)-x(2)*x(2)*x(5)-x(2)*x(5)*x(5)-e*f-f*g;

a=a1+a2;
e=e1+e2;
f=f1+f2;
g=g1+g2;

a1=1;
a2=1;
e1=9;
e2=7;
f1=56;
f2=98;
g1=76;
g2=665;

dx=dx(:);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
clc

x0=[1,1,1,1,1,1];
t0=0:0.1:2;
[t,x]=ode45('myfun',[0,10],x0); %?%ode45会自动调整步长
plot(t,x)                         % ?
legend('x(1)','x(2)','x(3)','x(4)','x(5)','x(6)')

a1=1;
a2=1;
e1=9;
e2=7;
f1=56;
f2=98;
g1=76;
g2=665;
a=a1+a2;
e=e1+e2;
f=f1+f2;
g=g1+g2;
% syms a e f g a1 a2 e1 e2 f1 f2 g1 g2