function [sqeuclid] = sqeuclidean(data)

[n,p] = size(data);

sqeuclid = zeros(n,n);

for i = 1:n
    for j = 1:n
        
        for k = 1:p
            
            sqeuclid(i,j) = sqeuclid(i,j) + (data(i,k) - data(j,k))^2;
            
        end
    end
end