function M=gaussbf(X,C,v)

%
% Gaussian Basis Function
%

[n,nc]=size(C);
k=size(X,2);

for xi=1:k
    for i=1:nc
        M(i,xi)=exp(-(1/2)*(norm((X(:,xi)-C(:,i))./v)^2));
    end
    % M(:,i)=M(:,i)/sum(M(:,i));
end