3y(y^2+1)=x^3+3x+1