One approach:
x = linspace(-5, 5, 50);
[X,Y] = meshgrid(x);
fcn = @(x,y,a,b) exp(-((x.^2+a) + (y.^2+b))*0.2); % Parameters ‘a’ & ‘b’ Are Offsets
figure(1)
mesh(X,Y,fcn(X,Y,0,0))
hold on
mesh(X+5,Y+5,fcn(X,Y,5,5)*12)
mesh(X+5,Y+10,fcn(X,Y,5,10)*10)
mesh(X+10,Y+5,fcn(X,Y,10,5)*10)
hold off
grid on
Experiment to get the result you want.
