Calculates thermonuclear reaction rates $\langle\sigma v\rangle$ for D-T, D-D, D-${}^3$He, and p-p fusion reactions. Implements three methods: - Numerical integration of the Gamow peak integrand (most general) - Analytic saddle-point approximation (good for constant S-factor) - Bosch-Hale parameteri