! Initializes arrays a and b.  
        ! Adds the arrays using 1. DO loop and 2. direct addition (a + b)
        PROGRAM addarrays
	IMPLICIT NONE
        INTEGER :: i
        REAL, DIMENSION(4) :: a = (/ 1., 2., 3., 4./)
        REAL, DIMENSION(4) :: b = (/ 1., 2., 3., 4./)
        REAL, DIMENSION(4) :: c,d
	DO i = 1, 4
	   c(i) = a(i) + b(i)
	END DO

	d = a + b
        WRITE(*,100) 'c', c
	WRITE(*,100) 'd', d
100     FORMAT (' ', A, ' = ', 5(F6.1,1X))
        END PROGRAM