Solving for conformal time. Need to work out how the RK4-solver behaves. Current task: correct subroutine input arguments