f(x,y) = y^2 - x^3+4x^2+1