u0, u1, 
u0^2+2*u1^2-u0
u0+2*u1-1
