PROGRAM KAHAN4 IMPLICIT NONE INTEGER, PARAMETER :: QP = SELECTED_REAL_KIND(33,3000) REAL (KIND=QP) :: A, B, C, D, E A = 1.0_QP/3.0_QP B = 4.0_QP*A - 1.0_QP C = 3.0_QP*B - 1.0_QP D = 0.5_QP*C E = ABS(D) WRITE(*,*) ' my = ', E END PROGRAM KAHAN4