Recent content by redgun

I decided to use gsl methods to because rk4 might cause some problems of stability. Now I try to solve Peebles ODE with this code : int peebles_ode(double x, const double y[], double f[], void* params) { Recombinaison* rec = (Recombinaison*)params; double a = exp(x); const double...

Hello everyone, I am trying to code the recombination history of the universe in C. More precisely I have difficulties to solve the peebles ode equation. In fact when I am close to the saha limit, I begin the resolution but the rhs of the ode seems to be too big and it makes my result Xe...