function c=martixmult(a,b)
% multiplication of two symetric matrices

N=size(a,1);

c=zeros(N);

for i=1:N
    for j=1:N
        for k=1:N
            c(i,j)=c(i,j)+a(i,k)*b(k,j);
        end
    end
end

