pure function funA_0_2overR2(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = OverR2 * funA(b, 0.0_dp, 2.0_dp) end function funA_0_2overR2