x, y , z = var('x,y,z')

plot3d(787 -x^2 - y^2, (x,-500, 500), (y,-500, 500))