Are any of these what you're looking for?
a1=rand(5,1)*3;
a2=rand(5,1)*3;
[x1, x2] = meshgrid(sort(a1),sort(a2)); % sort for for surf plot
ff= -sin(x1).*(sin(x1.^2./pi())).^2 -sin(x2).*(sin(2.*x2.^2./pi())).^2;
subplot(2,2,1);
plot(x1, x2, 'k.');
title('plot');
subplot(2,2,2);
plot3(x1, x2, ff, 'k.');
title('plot3');
subplot(2,2,3);
surf(x1,x2,ff);
title('surf');
subplot(2,2,4);
scatter(x1(:), x2(:), [], ff(:), 'filled');
title('scatter');
