仅供参考，没有看到ode45函数
syms y1(t) y2(t) y3(t) y4(t) y5(t) y6(t)
kappa=20;
eqns=[diff(y1,t)==-kappa*(3*y1-y2-y6-1/6*(4200*sin(t+pi)+4200)),...
diff(y2,t)==-kappa*(2*y2-y1-y3), diff(y3,t)==-kappa*(2*y3-y2-y4),...
diff(y4,t)==-kappa*(2*y4-y3-y5), diff(y5,t)==-kappa*(2*y5-y4-y6),...
diff(y6,t)==-kappa*(2*y6-y5-y1)];
cond=[y1(0)==0 y2(0)==0 y3(0)==0 y4(0)==0 y5(0)==0 y6(0)==0];
[y1sol(t),y2sol(t),y3sol(t),y4sol(t),y5sol(t),y6sol(t)]=dsolve(eqns,cond);
t=0:12;
plot(t,y1sol(t))
