Significant change in the direct scheme; mu and spr now evolve using the rates that are used in the whole function, instead of evolving by a similar proportion every time.