dx1/dt = p2*u1(t) - p1*x1,
dx2/dt = -p3*x2 + p4*u1(t), 
dx3/dt = p2*x2*u1(t) - p3*x3 + p4*u1(t)*x1 - x3*p1, 
y1 = x3
