Rather than

  (1-wc)*Ma + wc*Cv

we now have

  Ma + Cv

Also print information for time 0

Also change its default from 1e6 to 1e5 and treat it as size_t

The solvers for u0/a0 and v now each have parameters tolerance,
maximumIterations, verbosity