function [y,d]=sum1(x)
[r,c]=size(x);
if r==1
    a=0;
    for i=1:c
        a=a+x(i)^2;
    end
    y=a;
else
    for j=1:c
        a(j)=0;
        for i=1:r
            a(j)=a(j)+x(i,j)^2;
        end
        y(j)=a(j);
    end
end