function kernel1 = kernel(X,l,n,q);
%  Gauss Kernel Function;
for i = 1:l
    for j = 1:l
        temp = 0;
        for k = 1:n
            temp = temp + (X(i,k)-X(j,k))^2;
        end;
        kernel1(i,j)=exp(-q*temp);
    end;
end;