a = 0.0
  b = 1.6
  h = ( b - a ) / n
  quad = 0.0

  do i = 1, n

    x = ( ( n - i ) * a + ( i - 1 ) * b ) / ( n - 1 )  

    sump = (((((    &
            pp(1)   &
      * x + pp(2) ) &
      * x + pp(3) ) &
      * x + pp(4) ) &
      * x + pp(5) ) &
      * x + pp(6) ) &
      * x + pp(7)

    sumq = (((((    &
        x + qq(1) ) &
      * x + qq(2) ) &
      * x + qq(3) ) &
      * x + qq(4) ) &
      * x + qq(5) ) &
      * x + qq(6)
 
    fx = sump / sumq

    quad = quad + h * fx

  end do