func = lambdify(t, equation.rhs, 'numpy')
xvals = np.arange(0, 10, 0.1)
yvals = func(xvals)
fig, ax = plt.subplots(1, 1, subplot_kw=dict(aspect='equal'))
ax.plot(xvals, yvals)
ax.set_xlabel('t')
ax.set_ylabel('v(t)')
plt.show()
Plotting Sympy Result to Particular Solution of Differential Equation