I have just done it. Thank you.
rogram inverse
real(8), dimension(:,:), allocatable :: a, ID
real(8), dimension(:), allocatable :: b, x, y
real(8) :: s1, su, sm
integer :: i, j, k, n
print *, 'dimensió'
read *, n
allocate (a(n,n))
allocate (ID(n,n))
allocate (b(n))
allocate (x(n))
allocate...