Open in CoCalc with one click!
In [1]:
%display latex
In [2]:
var('u v')
(u,v)\left(u, v\right)
In [3]:
def W(x): if x < -1: return RDF(0) return RDF(lambert_w(RDF(x)/RDF(e))+1)
In [4]:
plot(W, (x,-2,5), aspect_ratio=1)
In [5]:
def h(u,v): return sinh(tan(u))*sinh(tan(v))-1
In [6]:
a = 0.99*pi/2
In [7]:
implicit_plot(h, (u,-a,a), (v,-a,a), plot_points=400, fill=True, fillcolor='yellow')
In [8]:
def g(u,v): return 4*W(-sinh(tan(u))*sinh(tan(v)))^2 * cos(u)^2 * cos(v)^2
In [9]:
graph = plot3d(g, (u,-a,a), (v,-a,a), plot_points=200)
In [10]:
graph += plot3d(lambda u,v: 0.001, (u,-a,a), (v,-a,a), color='yellow')
In [11]:
show(graph.rotate((0,0,1), pi/2), viewer='tachyon', view_point=(-3,3,2))
In [12]:
graph