C
do i = 1, N
if(i /= k) then
w = a(i, k)
do j = 1, M
a(i, j) = a(i, j) - w * a(k, j)
end do
end if
end do
write(*, *)
call matprint(a, n, m)
end do
end
C
subroutine matprint(matrix, n, m)
real matrix(n, m)
integer n, m, i, j
C
do i = 1, n
write(*, 600) (matrix(i, j), j= 1, M)
end do
600 format(1h , 10(G15.7:))
end