Photon (eq, radiation)
n = 3.03963550927e+19, rho = 4.11233516712e+26
--------------------------------------------------------------------------------
Electron (eq, intermediate)
n = 4.5581332389e+19, rho = 7.19109974007e+26
--------------------------------------------------------------------------------
Electron neutrino (non-eq, non-equilibrium)
n = 2.28363298548e+19, rho = 3.59827971959e+26
--------------------------------------------------------------------------------
Muon neutrino (non-eq, non-equilibrium)
n = 2.28363298548e+19, rho = 3.59827971959e+26
--------------------------------------------------------------------------------
Tau neutrino (non-eq, non-equilibrium)
n = 2.28363298548e+19, rho = 3.59827971959e+26
--------------------------------------------------------------------------------
Neutrino species scattering
	ν_e + ν_e ⟶  ν_e + ν_e	(K1=2.0 (0, 1, 2, 3), K1=4.0 (0, 3, 1, 2))
Neutrino species scattering
	ν_e + ν_μ' ⟶  ν_e + ν_μ'	(K1=1.0 (0, 3, 1, 2), K1=1.0 (0, 1, 2, 3))
	ν_μ + ν_e' ⟶  ν_e' + ν_μ	(K1=1.0 (0, 2, 1, 3), K1=1.0 (0, 1, 2, 3))
	ν_μ + ν_μ' ⟶  ν_e + ν_e'	(K1=1.0 (0, 3, 1, 2))
	ν_e + ν_e' ⟶  ν_μ + ν_μ'	(K1=1.0 (0, 3, 1, 2))
Neutrino species scattering
	ν_τ + ν_e ⟶  ν_e + ν_τ	(K1=1.0 (0, 1, 2, 3), K1=1.0 (0, 2, 1, 3))
	ν_e + ν_e' ⟶  ν_τ + ν_τ'	(K1=1.0 (0, 3, 1, 2))
	ν_e + ν_τ' ⟶  ν_e + ν_τ'	(K1=1.0 (0, 3, 1, 2), K1=1.0 (0, 1, 2, 3))
	ν_τ + ν_τ' ⟶  ν_e + ν_e'	(K1=1.0 (0, 3, 1, 2))
Neutrino species scattering
	ν_μ + ν_μ' ⟶  ν_μ + ν_μ'	(K1=4.0 (0, 3, 1, 2), K1=2.0 (0, 1, 2, 3))
Neutrino species scattering
	ν_τ + ν_μ ⟶  ν_μ + ν_τ	(K1=1.0 (0, 1, 2, 3), K1=1.0 (0, 2, 1, 3))
	ν_τ + ν_τ' ⟶  ν_μ + ν_μ'	(K1=1.0 (0, 3, 1, 2))
	ν_μ + ν_μ' ⟶  ν_τ + ν_τ'	(K1=1.0 (0, 3, 1, 2))
	ν_μ + ν_τ' ⟶  ν_μ + ν_τ'	(K1=1.0 (0, 3, 1, 2), K1=1.0 (0, 1, 2, 3))
Neutrino species scattering
	ν_τ + ν_τ' ⟶  ν_τ + ν_τ'	(K1=4.0 (0, 3, 1, 2), K1=2.0 (0, 1, 2, 3))
Neutrino pair annihilation into lepton pair
	ν_e + e' ⟶  e' + ν_e	(K1=2.35242752 (0, 2, 1, 3), K1=2.35242752 (0, 1, 2, 3), K2=-1.35242752 (1, 2, 0, 3))
	e + e' ⟶  ν_e + ν_e'	(K1=2.13861376 (0, 3, 1, 2), K1=0.21381376 (1, 3, 0, 2), K2=0.67621376 (0, 1, 2, 3))
	e + ν_e ⟶  e + ν_e	(K1=2.35242752 (0, 1, 2, 3), K1=2.35242752 (1, 2, 0, 3), K2=-1.35242752 (0, 2, 1, 3))
	ν_e + ν_e' ⟶  e + e'	(K1=2.13861376 (0, 3, 1, 2), K1=0.21381376 (0, 2, 1, 3), K2=0.67621376 (2, 3, 0, 1))
Neutrino pair annihilation into lepton pair
	ν_μ + e' ⟶  e' + ν_μ	(K1=0.50282752 (0, 2, 1, 3), K1=0.50282752 (0, 1, 2, 3), K2=0.49717248 (1, 2, 0, 3))
	e + ν_μ ⟶  e + ν_μ	(K1=0.50282752 (0, 1, 2, 3), K1=0.50282752 (1, 2, 0, 3), K2=0.49717248 (0, 2, 1, 3))
	e + e' ⟶  ν_μ + ν_μ'	(K1=0.28901376 (0, 3, 1, 2), K1=0.21381376 (1, 3, 0, 2), K2=-0.24858624 (0, 1, 2, 3))
	ν_μ + ν_μ' ⟶  e + e'	(K1=0.28901376 (0, 3, 1, 2), K1=0.21381376 (0, 2, 1, 3), K2=-0.24858624 (2, 3, 0, 1))
Neutrino pair annihilation into lepton pair
	ν_τ + e' ⟶  e' + ν_τ	(K1=0.50282752 (0, 2, 1, 3), K1=0.50282752 (0, 1, 2, 3), K2=0.49717248 (1, 2, 0, 3))
	ν_τ + ν_τ' ⟶  e + e'	(K1=0.28901376 (0, 3, 1, 2), K1=0.21381376 (0, 2, 1, 3), K2=-0.24858624 (2, 3, 0, 1))
	e + e' ⟶  ν_τ + ν_τ'	(K1=0.28901376 (0, 3, 1, 2), K1=0.21381376 (1, 3, 0, 2), K2=-0.24858624 (0, 1, 2, 3))
	e + ν_τ ⟶  e + ν_τ	(K1=0.50282752 (0, 1, 2, 3), K1=0.50282752 (1, 2, 0, 3), K2=0.49717248 (0, 2, 1, 3))
dx = 0.0102542192752 MeV
#step	Time, s	aT, MeV	T, MeV	scale factor	dx, MeV
#0 	t = 0.0 	aT = 1.0 	T = 5.0 	a = 0.2 	dx = 0.0102542192752