var('u v')
b=0.4
r=1-b^2
w=sqrt(r)
def dihalica(u,v):
xk=-u+(2*r*cosh(b*u)*sinh(b*u))/(b*((w*cosh(b*u))^2 + (b*sin(w*v))^2))
yk=(2*w*cosh(b*u)*(-(w*cos(v)*cos(w*v))-sin(v)*sin(w*v)))/(b*((w*cosh(b*u))^2 + (b*sin(w*v))^2))
zk=(2*w*cosh(b*u)*(-(w*sin(v)*cos(w*v))+cos(v)*sin(w*v)))/(b*((w*cosh(b*u))^2 + (b*sin(w*v))^2))
return yk,zk,xk
slika1=parametric_plot3d(dihalica(u,v),(u,-10,10), (v,-37.4,37.4),plot_points=140, color='lemonchiffon')
slika1.show(viewer='threejs',online=True,frame=False)