Support derivatives with respect to `u0` in the forward modes

This isn't too hard: set the initial conditions in the forward propagation or use a larger set of dual numbers in the ForwardDiff form. It just needs to get setup.