p^2 = z;
2*p*p' = z';
z'*y - z = y^2;
z' - (1/y)*z = y.