Ucostante(x,y):=10;

U(x,y):=2*(x+1)^0.3*(y+1)^0.5;

plot3d([[Ucostante, [x,0,20], [y,0,20]],
U, [x,0,20], [y,0,20]], 
[plot_format,xmaxima])$