1 | initial version |
I do not think there is a specific command. Here is a potential work around:
var('x y z')
f = z==3*(1-x)^2*exp(-(x^2) - (y+1)^2) - 10*(x/5 - x^3 - y^5)*exp(-x^2-y^2)- 1/3*exp(-(x+1)^2 - y^2)
levels=srange(-4,4,0.5)
p=Graphics()
for zlevel in levels:
p+=implicit_plot3d(f.subs(z==zlevel), (x,-4,4),(y,-4,4),(z,zlevel-0.05,zlevel+0.05),thickness=3)
show(p)