LmH = 290.75 (slider from 50 to 500, step 0.05)
CnF = 250.25 (slider from 50 to 500, step 0.05)
RkOhm = 2.05 (slider from 0 to 10, step 0.05)
Q: sqrt(LmH / CnF) / RkOhm
ω_{0}: 1000 / sqrt((LmH * CnF))
μ: ω_{0} / ((2 * Q))
ω_p: (ω_{0} * sqrt(1 - 1 / ((4 * Q^(2)))))
E = 3.2 (slider from -5 to 5, step 0.05)
r_1: (-(ω_{0} / ((2 * Q)))) + (ω_{0} * sqrt(1 / ((4 * Q^(2))) - 1))
r_2: (-(ω_{0} / ((2 * Q)))) - (ω_{0} * sqrt(1 / ((4 * Q^(2))) - 1))
u(t) = If[t > 0, If[Q > 0.5, (((-E) * (cos((ω_p * t)) + (μ / ω_p * sin((ω_p * t))))) * ℯ^(((-μ) * t))) + E, ((r_2 / (r_1 - r_2) * E) * ℯ^((r_1 * t))) - ((r_1 / (r_1 - r_2) * E) * ℯ^((r_2 * t))) + E], 0]
u' = Derivative(u)
t_0 = 2.18 (slider from 0 to 15, step 0.02)
M: (t_0, u(t_0))
M': (u(t_0), u'(t_0))
A: (t_0, 0)
a = Line(M, A)
TrajPhase = CurveCartesian(u(t), u'(t), t, 0, 50)