started updating LQR and LQ classes to take S matrix and deltas in q(t), i.e. jump terms in p(t)