n=500;
X = rand(n,1)*20 - 10;
Y= rand(n,1)*20 - 10;
R = sqrt(X.^2 + Y.^2) + eps;
Z = sin(R)./R;
Vx = linspace(min(X),max(X),40);
Vy = linspace(min(Y),max(Y),40);
[X2,Y2] = meshgrid(Vx,Vy);
Z2 = griddata(X,Y,Z,X2,Y2,'cubic');
figure(1)
mesh(X2,Y2,Z2) 
axis tight; hold on
plot3(X,Y,Z,'*') 
hold off








