Master

Implemented rk4 solver, to be able to use higher than O(2) solver for simulations with temperature. 

It compiles and I tested it in simulations with fixed timestep, adaptive timestep and with temperature and it gives the same results as the other solvers.