plot y^2 arctan{x / (y sqrt(x^2+y^2+1))