求助：ode求解微分方程组

stayrong

用于求解方程如下:
function dx=ehpcs(t,x)
w=15*pi;
k=300;
J1=0.191;
J2=0.097;
dx=zeros(4,1);
T0=40;
dx(1)=x(2);
dx(2)=T0*sin(4*w*t)/J1-(k/J1)*(x(1)-x(3));
dx(3)=x(4);
dx(4)=-(k/J2)*(x(3)-x(1));

%ode求解程序如下：
tspan=[0:0.001:10];
x0=[0;0;0.8;0];
[t x]=ode45('ehpcs',tspan,x0);
plot(t,x(:,1),'r',t,x(:,3),'b-.')
grid

微分方程就是最简单的线性微分方程组，是实际应用的减振器的简化模型，但是求解出来的结果不是围绕x轴的波动数据，本人刚接触ode求解，自己没有检查出程序错误，希望有会的大神帮我看看是哪里出了错，不甚感激！