In [1]:

%matplotlib inline

import numpy as np
import matplotlib.pyplot as plt
import sympy as sym

In [3]:

from sympy.abc import x

In [4]:

sym.init_printing()

In [5]:

f = 1 + 3*x + x**2

In [6]:

g = 1 + 4*x + 3*x**2

In [7]:

f*g

Out[7]:

$$\left(x^{2} + 3 x + 1\right) \left(3 x^{2} + 4 x + 1\right)$$

In [9]:

sym.expand(f*g)

Out[9]:

$$3 x^{4} + 13 x^{3} + 16 x^{2} + 7 x + 1$$

In [10]:

h = f*g

In [13]:

type(h)

Out[13]:

sympy.core.mul.Mul

In [16]:

hc = sym.collect(h, x)

In [19]:

type(hc)

Out[19]:

sympy.core.mul.Mul

In [25]:

h.subs(x,3)

Out[25]:

$$760$$

In [26]:

sym.factorial(5)

Out[26]:

$$120$$

In [ ]: