x= -1:0.05:1;
y=0:0.05:1;
[X,Y]=meshgrid(x,y);
z=4.*sin(2.*pi.*X).*cos(1.5.*pi.*Y).*(1-X.^2).*Y(1-Y);
>> figure
>> subplot(2,2,1),plot3(X,Y,z),title('Plot3')
>> subplot(2,2,2),mesh(X,Y,z),title('mesh')
>> subplot(2,2,3),surf(X,Y,z),title('surf')