function C=matrixmult(A,B)
[n1,m1]=size(A);
[n2,m2]=size(B);
C=zeros(n1,m2);
if m1~=n2
    error('Matrizen sind nicht kompatibel:')
end
for i=1:n1
    for j=1:m2
        C(i,j)=0;
        for k=1:m1
            C(i,j)=C(i,j)+A(i,k)*B(k,j);
        end
    end
end


end