pub fn solve_dyn<T>( A: &DMatrix<T>, b: &DVector<T>, c: &DVector<T>, eps: T, theta: T, gamma: T, max_it: usize ) -> (DVector<T>, DVector<T>)where T: Scalar + RealField + Copy,