Parameters:
a=0.05
ae2=0.005
ae5=0.005
ae=0.0033
ah=0.0036
te=12
th=100
tau=1

Dynamics:
dx_0=(-a-ae)*x_0+a*x_1+ae*te
dx_1=(-4*a-ae2-ah*u_0)*x_1+a*x_0+a*x_6+a*x_8+a*x_2+ae2*te+ah*th*u_0
dx_2=(-2*a-ae)*x_2+a*x_1+a*x_3+ae*te
dx_3=(-2*a-ae)*x_3+a*x_2+a*x_4+ae*te
dx_4=(-4*a-ae5-ah*u_1)*x_4+a*x_3+a*x_7+a*x_5+a*x_9+ae5*te+ah*th*u_1
dx_5=(-a-ae)*x_5+a*x_4+ae*te
dx_6=(-a-ae)*x_6+a*x_1+ae*te
dx_7=(-a-ae)*x_7+a*x_4+ae*te
dx_8=(-a-ae)*x_8+a*x_1+ae*te
dx_9=(-a-ae)*x_9+a*x_4+ae*te