program matrix
      
C************************************************************
C  Create 2 5x3 matrices, Find sum, use series/equation to file
C************************************************************

C-----SPECIFICATIONS
      real matrixA(5,3), matrixB(5,3), matrixSum(5,3)
      
      
C-----FILL ARRAYS
      call fillArray(matrixA,matrixB)
      
      do 100 nrow = 1, 5
         print *, (matrixA(nrow,i),i=1,3),
     +            '   ', (matrixB(nrow,i),i=1,3)
 100  end do
 
C-----SUM UP MATRICES
      do 200 nrow = 1, 5
         do 210 ncol = 1, 3
	     matrixSum(nrow,ncol) = 
     +          matrixA(nrow,ncol) + matrixB(nrow,ncol)
 210     end do
 200  end do
 
      print *, ' '
      do 300 nrow = 1, 5
           print *, (matrixSum(nrow,i),i=1,3)
 300  end do
      
      stop
      end
C---------------------------------------------------------------
      subroutine fillArray(matrixA, matrixB)
      
C-----A series:  1, 2, 3, 5, 8, 11, ...

      real matrixA(5,3), matrixB(5,3), numA, i, j, numB
      
      numA = 1.0
      numB = 1.0
      do 100 nrow = 1, 5
         do 110 ncol = 1, 3
	 
	    matrixA(nrow,ncol) = numA
	    
	    if(nrow.eq.1.and.ncol.eq.1) then
	       numA = numA + numA
	    else
	       j = nrow
	       i = ncol
	       i = i - 1
	       if (i.eq.0) then
	          i = 3
		  j = j - 1
	       end if
	       numA = matrixA(nrow,ncol) + matrixA(j,i)
	    end if
	    
     
            matrixB(nrow,ncol) = numB
	    numB = numB + 0.3 * (numB * 2)
	 
	    
 110     end do
 100  end do
      
      return
      end