SymPy recently got a new Linear system solver: linsolve in sympy.solvers.solveset, you can use that as follows: In [38]: from sympy import * In [39]: from sympy.solvers.solveset import linsolve In [40]: x, y, z = symbols(‘x, y, z’) List of Equations Form: In [41]: linsolve([x + y + z – 1, x + y + … Read more
(Full disclosure: I am the lead developer of SymPy) The first thing you should understand is that SymPy and Sage are not quite the same thing. SymPy is a pure Python library, that does computer algebra. Sage is a collection of open source mathematical software. Sage tries to gather together all the major open source … Read more
You can use: f.free_symbols which will return a set of all free symbols. Example: >>> import sympy >>> x, y, z = sympy.symbols(‘x:z’) >>> f = sympy.exp(x + y) – sympy.sqrt(z) >>> f.free_symbols set([x, z, y])
all_coeffs() can be sometime better than using coeffs() for a Poly. The difference lies in output of these both. coeffs() returns a list containing all coefficients which has a value and ignores those whose coefficient is 0 whereas all_coeffs() returns all coefficients including those whose coefficient is zero. >>> a = Poly(x**3 + a*x**2 – … Read more
you need to use display: from IPython.display import display display(yourobject) It will choose the appropriate representation (text/LaTex/png…), in recent enough version of IPython (6.0+) display is imported by default, still we recommend to explicitly import it.
for numerical solution, you can use fsolve: http://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.fsolve.html#scipy.optimize.fsolve from scipy.optimize import fsolve import math def equations(p): x, y = p return (x+y**2-4, math.exp(x) + x*y – 3) x, y = fsolve(equations, (1, 1)) print equations((x, y))