a, z, y1, y2, 
y1-2/3*a*(-z^4+z)
y2*z^2-1/2*a*(z^4-1)
