Solving DE's with ode45 and ode23