program matrix
implicit none

! program oblicza elementy macierzy 3x3 zdefiniowanej jako A = 1 + 2B
! gdzie 1 - macierz jednostkowa
!       B - dowolna macierz 3x3

!deklaracja zmiennych

real, dimension(3,3) :: a, b, identity
integer :: i, j

!wczytywanie elementow macierzy B

do i = 1, 3
   do j = 1, 3

      write(*,*) 'Wpisz', i, j, 'element macierzy B'
      read(*,*) b(i, j)
      
   end do
end do

!tworzenie macierzy jednostkowej

do i = 1, 3
   do j = 1,3
      if (i == j) then

         identity(i, j) = 1.0

      else
         
         identity(i, j) = 0.0

      end if
   end do
end do


! obliczanie macierzy A

do i = 1, 3
   do j = 1, 3

      a (i, j) = identity(i, j) + (2 * b(i, j))

   end do
end do

! wypisanie macierzy A

do i = 1, 3

   write(*,*) (a(i, j), j = 1,3)

end do

end