f=21.5+x1.*sin(4*pi*x1)+x2.*sin(20*pi*x2);
mesh(x1,x2,f),grid on;xlabel('x1'),ylabel('x2'),zlabel('f(x1,x2)');

数组的乘、除、乘方要用点运算,以此与矩阵运算区别。
plot是二维数据绘图,mesh和surf是三维数据绘图,区别主要是网格线颜色和其之间是否填色的不同。
grid on 是生成坐标网格,易于你读取坐标数据。
你还是参考help文档吧,这些实在太基础没啥好说的,你通过简单的试验一看就懂!