Minor fixes to have make code work
Numerical model for Siudem's analytical model