function C=mult(A,B)

%This program computes the product of two
%square matrices

[m n]=size(A);

for j=1:n
    for k=1:n
        C(j,k)=0;
        for m=1:n
            C(j,k)=C(j,k)+A(j,m)*B(m,k);
        end
    end
end