solving ivp method