using hinit function to determine initial integration step
added experimental data of april measurements