Several ways: mesh(), surf(), contoruf(), contour(), pcolor(), imagesc()
[X, Y] = meshgrid(linspace(-1, 1, 50));
Z = exp(X.^2+Y.^2-4*X);
f = figure();
subplot(2,2,1);
surf(X, Y, Z);
shading interp
subplot(2,2,2);
mesh(X, Y, Z);
subplot(2,2,3);
contourf(X, Y, Z, 10);
subplot(2,2,4);
pcolor(X, Y, Z);
shading interp
Or fsurf(), fmesh(), fcontourf()
fun = @(X, Y) exp(X.^2+Y.^2-4*X);
f = figure();
subplot(2,2,1);
fsurf(fun, [-1 1 -1 1]);
shading interp
subplot(2,2,2);
fmesh(fun, [-1 1 -1 1]);
subplot(2,2,3);
fcontour(fun, [-1 1 -1 1], 'Fill', 'on');
subplot(2,2,4);
fcontour(fun, [-1 1 -1 1]);