Advertisement
Guest User

Untitled

a guest
May 5th, 2019
89
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. PROGRAM EX7_50
  2.    implicit none
  3.  
  4.    integer(4), parameter      :: R_ = 4
  5.  
  6.    character(*), parameter    :: input = "../data/input", output = "../data/output"
  7.    integer(R_)                :: in = 1, out = 2, N, i
  8.  
  9.    integer(R_), allocatable   :: A(:,:), B(:,:)
  10.  
  11.    open(newunit = in, file = input)
  12.       read(in,*) N
  13.       allocate(A(N,N))
  14.       read(in,*) (A(i,:), i = 1, N)
  15.    close(in)
  16.    write(*,*) A
  17.  
  18.    allocate(B(N,N-1))
  19.  
  20.    do concurrent (i = 1:N)
  21.       B(i, i:N-1) = A(i, i+1:N)
  22.       B(i, 1:i-1) = A(i, 1:i-1)
  23.    end do
  24.  
  25.    write(*,*) (B(i,:), i = 1, N)
  26.    
  27. end program
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement