#pragma once
void FindNeutralsEffect(double f, double theta, int sigma, double f_p, double f_B, double T0, double n_e, double n_H, double n_He,
                        double *j, double *k);