small changes to case study code
uses only matrix operations for K time step calc, x100 speed improvement
Used matrix methods for Q time stepping too