from sympy import *
from spb import *
x, y, a, b, c, d, e = symbols("x, y, a, b, c, d, e")
expr = Eq(a * x**2 - b * x + c, d * y + y**2)
graphics(
    implicit_2d(expr, (x, -2, 2), prange(y, -e, e),
        params={
            a: (10, -15, 15),
            b: (7, -15, 15),
            c: (3, -15, 15),
            d: (2, -15, 15),
            e: (10, 1, 15),
        }, n=150),
    ylim=(-10, 10))