pure function funA_2_1_overR2(b) result(ret) real(dp), intent(in) :: b real(dp) :: ret ret = OverR2 * funA_2_1(b) end function funA_2_1_overR2