我有一个微分方程组
dx1/dt=x3;
dx2/dt=x4;
dx3/dt=10*x3*x4*(x1-x2)+10*x4^2*(x1-x2)-20*x3^2*(x1-x2)+100*cos(20*t+30)-50*cos(20*t+45);
dx4/dt=20*x3^2*(x1-x2)-30*x4^2*(x1-x2)-10*x3*x4*(x1-x2)+25*cos(20*t+45)-20*cos(20*t+30);

x1(0)=0;
x2(0)=0;
x3(0)=0;
x4(0)=0;
dx1/dt(t=0)=dx2/dt(t=0)=dx3/dt(t=0)=dx4/dt(t=0)=0;

因为是个新手，折腾了两天了，看了n个帮助文件了，但还是不知道怎么处理，要是谁能给我指点一下，多谢！！

那就光用这个初始条件：
x1(0)=x2(0)=x3(0)=x4(0)=0;

还是不行啊，斑竹能不能帮我写一个程序框架，多谢！！！