! y = a * x + y
PROGRAM exd_2

  IMPLICIT NONE
  INTEGER, PARAMETER :: n = 5, &         ! ベクトルの配列の大きさ
                        incx = 1, &      ! xの増分，1でよい
                        incy = 1         ! yの増分，1でよい
  INTEGER :: i
  REAL :: a = 3E0, &    ! スカラ
          x( 1 : n ), & ! ベクトル
          y( 1 : n )    ! ベクトル

  x(:) = 0E0
  DO i = 1, n
    x(i) = REAL( i )
  END DO
  y(:) = 0E0
  CALL saxpy( n, a, x, incx, y, incy )   ! ライブラリを呼び出す

  PRINT '( 5F7.3 )', x, y

END PROGRAM exd_2