coding and plotting of gaussian indicator functions

Hello friends,

I am working on palmprint recognition system using spectral minutiae representation.

Using code available here , I have extracted various minutiae i.e ridge termination points and ridge bifurcation points from a segment of palmprint image

http://in.mathworks.com/matlabcentral/fileexchange/16728-fingerprint-application/content/FingerPrint/html/fingerprint.html

the original image is :

whereas image showing minutiae points is:

CentroidTerm and CentroidBif contains the x and y coordinates of minutiae .

now I want to code following gaussian indicator function

where xi and yi are the x and y coordinates of minutiae points. and z is the number of minutiae points i.e I have to plot the summation of gaussian indicator function at all the locations of minutiae.

I have tried the following code(I have tried only CentroidTerm ,not CentroidBif, in following code):

    X= 1:1:1116;
Y= 1:1:404;
[x,y]= meshgrid(X,Y);  % //mesh
 F  = zeros(size(x));
 for q= 1:size(CentroidTerm,1)
    F= F + (exp(-((x-CentroidTerm(q,1)).^2)+((y-CentroidTerm(q,2)).^2)/(2*.001)))/(2*3.14*.001);
end
surf(X,Y,F); 
shading interp

but iam not getting the required plot..

All I am getting is F as 404X1116 double matrix of infinite values and a 3D grid with no gaussian

what values should be taken for x and y..?

The resulting plot would be a grid of gaussians

Waiting for your response.

Any help would be highly appreciated.