# qromb_p (ul_physlib)

value.r4.v = qromb_p(func.f.v, int_lower_bnd.r4.v,
int_upper_bnd.r4.v)
This routine performs a Romberg integration of a function
specified by passed func. It returns the value of the
integral integrated from integral_lower_bnd to integral_upper_bnd.
It is identical to the QROMB.C routine from Numerical Recipes in C
except that the requirements on the fractional accuracy desired
has been increased.
func function containing the integrand
(passed to trapzd function)
int_lower_bnd lower limit of integration
int_upper_bnd upper limit of integration
This function returns values as follows:
0.0 error occurred
otherwise integral value
This function requires the following include files:
ul_physlib_h:physlib_h
Related functions:
calc_sigma, emitt_dp, calc_lum