test_analytical_formula Subroutine

  SUBROUTINE test_analytical_formula()
    IMPLICIT NONE
    REAL(rp) :: z
    INTEGER :: ii

    do ii=1_idef,100
       z = 2.5_rp + REAL((ii-1_idef),rp)*0.5_rp
       write(6,*) IntK(5.0_rp/3.0_rp,z)
    end do
  END SUBROUTINE test_analytical_formula