dx1/dt = -b*x3*x1*Ninv - x1*Ninv*u1;
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
