real,allocatable:: a(:,:),b(:,:),c(:,:) open(10,file="test_3.in") read (10,*) n allocate (a(n,n),b(n,n),c(n,n)) read (10,*)a ,b do i=1,n do j=1,n sum=0. do k=1,n sum=sum+a(i,k)*b(k,j) enddo c(i,j)=sum enddo enddo do i=1,n print '(100f10.5)',c(i,:) enddo stop end