from sympy import *
from spb import *
x = symbols("x")
c = S(2) / 10
p = graphics(
    line(cos(x) * exp(-c * x), (x, 0, 10), label="oscillator"),
    line(exp(-c * x), (x, 0, 10), label="upper limit",
        rendering_kw={"linestyle": ":"}),
    line(-exp(-c * x), (x, 0, 10), label="lower limit",
        rendering_kw={"linestyle": ":"}),
    grid=False
)