dx1/dt = p2*u1 - p1*x1;
dx2/dt = -p3*x2 + p4*u1; 
dx3/dt = p2*x2*u1 - p3*x3 + p4*u1*x1 - x3*p1; 
y1 = x3
