请教

请教一下，如何实现下面程序
假设某一模型方程为:
  dx1/dt =3x1-2x2+4    (1)
  dx2/dt=-x1+x2+7      (2)
需要控制x1,设x1的设定值为xs,偏差信号为e=xs-x1,偏差信号的变化率为ec,控制信号u=f(e,ec),设仿真时间为t=0:3600,请问ec怎么表示