PROGRAM create_input 
      IMPLICIT NONE
      INTEGER :: n, i, j
      REAL*8, allocatable :: A(:,:)
      CHARACTER*80 :: argc1
      REAL*8   Rand

      call getarg(1,argc1)
      read(argc1,*) n
      print *, "Array size: ", n
      allocate(A(n,n))

!     Initialize array using random numbers
      call SRand(1001)
      DO j = 1, n
         DO i = 1, n
            A(i,j) = Rand()
         END DO
      END DO

      open(60, file='input')
      write(60,*) A

      END PROGRAM create_input