from sympy import *
init_printing(use_unicode=True, wrap_line=False, no_global=True)
import matplotlib.pyplot as plt
%matplotlib notebook
from sympy.abc import x,y
denominator = x**4 + 2*x**2 + 1
factor(denominator)
plot_implicit(Eq(denominator, y), (x, -0.6, 0.6), (y, 0, 2))
<sympy.plotting.plot.Plot at 0xc5c1908>
plot_implicit(Eq(denominator, y), (x, -6, 6), (y, 0, 1400))
<sympy.plotting.plot.Plot at 0xc668ef0>
denominator = 1 / (x**4 + 2*x**2 + 1)
factor(denominator)
plot_implicit(Eq(denominator, y), (x, -6, 6), (y, 0, 1))
<sympy.plotting.plot.Plot at 0xc605358>
numerator = x**3 * sin(x)**2
numerator
plot_implicit(Eq(numerator, y), (x, -6, 6), (y, -120, 120))
<sympy.plotting.plot.Plot at 0xcb09eb8>
wifieqn = numerator*denominator
factor(wifieqn)
plot_implicit(Eq(wifieqn, y), (x, -6, 6), (y, -0.4, 0.4))
<sympy.plotting.plot.Plot at 0xc64b898>
integ = integrate(wifieqn, (x, -5, 5))
integ
integ.doit()