After solving the O.D.E., a finite number of unknown integration
constants remain. In practical applications, these integration
constants are typically found from initial conditions
at a starting point such as t=0 or from boundary conditions
at the end points a and b of an interval