L'implementazione (i,k,j)

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