Solving with Equality

The code from this video is listed below.

import sympy as sp
x = symbols('x')
f = x ** 2 + 4 * x - 16
sp.solve(sp.Eq(f, 5), x)
sp.plot(f - 5)