In [1]:
import sympy as s
In [ ]:
In [15]:
s.init_printing()
In [16]:
x = s.Symbol('x')
y = s.Symbol('y')
In [17]:
2 * x + y + y
Out[17]:
$$2 x + 2 y$$
In [18]:
t = (x + y)**17
In [19]:
t
Out[19]:
$$\left(x + y\right)^{17}$$
In [20]:
t.expand()
Out[20]:
$$x^{17} + 17 x^{16} y + 136 x^{15} y^{2} + 680 x^{14} y^{3} + 2380 x^{13} y^{4} + 6188 x^{12} y^{5} + 12376 x^{11} y^{6} + 19448 x^{10} y^{7} + 24310 x^{9} y^{8} + 24310 x^{8} y^{9} + 19448 x^{7} y^{10} + 12376 x^{6} y^{11} + 6188 x^{5} y^{12} + 2380 x^{4} y^{13} + 680 x^{3} y^{14} + 136 x^{2} y^{15} + 17 x y^{16} + y^{17}$$
In [23]:
s.printing.ccode(t.expand())
Out[23]:
'pow(x, 17) + 17*pow(x, 16)*y + 136*pow(x, 15)*pow(y, 2) + 680*pow(x, 14)*pow(y, 3) + 2380*pow(x, 13)*pow(y, 4) + 6188*pow(x, 12)*pow(y, 5) + 12376*pow(x, 11)*pow(y, 6) + 19448*pow(x, 10)*pow(y, 7) + 24310*pow(x, 9)*pow(y, 8) + 24310*pow(x, 8)*pow(y, 9) + 19448*pow(x, 7)*pow(y, 10) + 12376*pow(x, 6)*pow(y, 11) + 6188*pow(x, 5)*pow(y, 12) + 2380*pow(x, 4)*pow(y, 13) + 680*pow(x, 3)*pow(y, 14) + 136*pow(x, 2)*pow(y, 15) + 17*x*pow(y, 16) + pow(y, 17)'
In [24]:
t
Out[24]:
$$\left(x + y\right)^{17}$$
In [32]:
c = t.subs(x, 3).subs(y, -1/2)
In [33]:
c
Out[33]:
$$5820766.09134674$$
In [34]:
type(c)
Out[34]:
sympy.core.numbers.Float
In [31]:
type(int(c))
Out[31]:
int
In [35]:
float(c)
Out[35]:
$$5820766.091346741$$
In [37]:
c.evalf(50)
Out[37]:
$$5820766.09134674072265625$$
In [38]:
p = s.pi
In [39]:
p
Out[39]:
$$\pi$$
In [40]:
float(p)
Out[40]:
$$3.141592653589793$$
In [42]:
p.evalf(500)
Out[42]:
$$3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337867831652712019091456485669234603486104543266482133936072602491412737245870066063155881748815209209628292540917153643678925903600113305305488204665213841469519415116094330572703657595919530921861173819326117931051185480744623799627495673518857527248912279381830119491$$
In [43]:
from sympy import limit, sin, oo
In [44]:
oo
Out[44]:
$$\infty$$
In [47]:
limit(1/x, x, 0)
Out[47]:
$$\infty$$
In [49]:
limit(sin(x)/x**2, x, 0)
Out[49]:
$$\infty$$
In [50]:
from sympy import series
In [52]:
series(sin(x)/x, x, 0)
Out[52]:
$$1 - \frac{x^{2}}{6} + \frac{x^{4}}{120} + \mathcal{O}\left(x^{6}\right)$$
In [53]:
from sympy import integrate
In [54]:
a, b = s.symbols('a, b')
In [59]:
integrate(2*x, (x, 0.1, 2))
Out[59]:
$$3.99$$
In [64]:
r = (x + 2)*(x - a)
In [65]:
r.expand()
Out[65]:
$$- a x - 2 a + x^{2} + 2 x$$
In [66]:
s.solve(r, x)
Out[66]:
$$\left [ -2, \quad a\right ]$$
In [ ]: