Plot3D deboxed

edit

This code works nicely

var("x,y,z")
f(x,y) = (.5*x^0.4+.5*y^0.4)
z1, z2 = 0, 1.3
ni = 5
dz = (z2-z1)/ni
levels = [z1,z1+dz..z2]
curves = contour_plot(f(x,y),(x,.001,2), (y,.001,2), contours=levels, 
                      cmap="gist_earth")
def col2(x,y):
    return float(0.5*dz+floor((f(x,y)-z1)/dz)/ni)
surf = plot3d(f(x,y), (x,.001,2), (y,.001,2), color=(col2,colormaps.gist_earth), plot_points=300)
arrow1= arrow3d((0,0,0), (2,0,0), 1)
arrow2= arrow3d((0,0,0), (0,2,0), 1)
arrow3= arrow3d((0,0,0), (0,0,2), 1)
t1=text3d("x", (2.15,0,0), color=(0.5,0,0))
t2= text3d("y", (0,2.15,0), color=(0.5,0,0))
t3= text3d("EU", (0,0,3.), color=(0.5,0,0))
offset = 0
plane = plot3d(offset, (x,.001,2), (y,.001,2), color=(col2,colormaps.gist_earth), plot_points=300,viewer='threejs')
show(surf+plane+arrow1+arrow2+arrow3+t1+t2+t3)

I know that I could not have presently LaTeX label but I wonder if there is a way to debox it.