dx1/dt = -b*x3*x1*Ninv - x1*Ninv*u1(t),
dx2/dt = b*x3*x1*Ninv - a*x2,
dx3/dt = -x3*g + s*x4 + a*x2,
dx4/dt = -e*x3*g + x3*g - s*x4,
y1 = x3*Ninv
