from sympy import *
x, y, z = symbols("x y z")
eq= Eq(x**2+y**2+z**2-16,0)
s=solve(eq,z)
plotting.plot3d(*s)
The above code produces the following output.
How can I remove all the zeros plotted where z coordinate value is not defined?
In my actual program, s is a list of solutions obtained from sympy.solve so manually specifying parametric surface is not an option.