about to rejig RK4IP method to match Balac and Mahe 2012 paper