Adding ab09bd as a reduction function

Updating with `ab09bd` which also reduces the model order, like `ab09ad` but tries to keep better stability by using Singular Perturbation Approximation (SPA).