[x,y]=meshgrid([-5:0.5:5],[0:0.5:10]);
z=x.^3*y+sqrt(y);
surf(x,y,z)
xlim([-5,5]);
ylim([0,10]);