Eigen library replaced with Seldon. Program changed to work with momentum (q) instead of velocity (u).