原方程组如下：dx(t)/dt =rx(t）-ax(t)y(t)
              dy(t)/dt=-sx(t)+bx(t)y(t)   
其中：
     r=0.001,  a=2x10^(-6),  x(0)=12000
     s=0.01,   b=1X10^(-6),  y(0)=600
编写仿真程序用于仿真。
请高手多多帮助.