function [rc,cc] = singspotcent(A)
%SINGSPOTCENT returns centroid location for first image in input 3D matrix
MB=double(A(:,:,1));
[rn cn]=size(MB);
M=MB-mean(mean(min(MB)));
r=1;
c=1;
sumIc=0;
sumIr=0;
while c<(cn+1)
    sumpIr=0;
    i=1;
    while i<(rn+1)
    sumpIr=sumpIr+i*M(i,c)^2;
    i=i+1;
    end
    sumIr = sumIr + sumpIr;
    c=c+1;
end
c=1;
while r<(rn+1)
    sumpIc=0;
    i=1;
    while i<(cn+1)
    sumpIc=sumpIc+i*M(r,i)^2;
    i=i+1;
    end
    sumIc = sumIc + sumpIc;
    r=r+1;
end
sumI=sum(sum(M.^2));
rc=sumIr/sumI;
cc=sumIc/sumI;
end