// Solve linear system
//////////
-ex lsolve(const ex &eqns, const ex &symbols)
+ex lsolve(const ex &eqns, const ex &symbols, unsigned options)
{
// solve a system of linear equations
if (eqns.info(info_flags::relation_equal)) {
matrix solution;
try {
- solution = sys.solve(vars,rhs);
+ solution = sys.solve(vars,rhs,options);
} catch (const std::runtime_error & e) {
// Probably singular matrix or otherwise overdetermined system:
// It is consistent to return an empty list
/** Order term function (for truncated power series). */
DECLARE_FUNCTION_1P(Order)
-ex lsolve(const ex &eqns, const ex &symbols);
+ex lsolve(const ex &eqns, const ex &symbols, unsigned options = determinant_algo::automatic);
/** Check whether a function is the Order (O(n)) function. */
inline bool is_order_function(const ex & e)