These are fortran codes and gnuplot scripts that i have written for my astrophysics researches which included 4th order runge kutta method with adaptive step size.