Using Sympy Equations for Plotting

MANA624 picture MANA624 · Feb 14, 2016 · Viewed 18.4k times · Source

What is the best way to create a Sympy equation, do something like take the derivative, and then plot the results of that equation?

I have my symbolic equation, but can't figure out how to make an array of values for plotting. Here's my code:

from sympy import symbols
import matplotlib.pyplot as mpl

t = symbols('t')
x = 0.05*t + 0.2/((t - 5)**2 + 2)

nums = []
for i in range(1000):
    nums.append(t)
    t += 0.02

plotted = [x for t in nums]

mpl.plot(plotted)
mpl.ylabel("Speed")
mpl.show()

In my case I just calculated the derivative of that equation, and now I want to plot the speed x, so this is fairly simplified.

Answer

user picture user · Feb 14, 2016

You can use numpy.linspace() to create the values of the x axis (x_vals in the code below) and lambdify().

from sympy import symbols
from numpy import linspace
from sympy import lambdify
import matplotlib.pyplot as mpl

t = symbols('t')
x = 0.05*t + 0.2/((t - 5)**2 + 2)
lam_x = lambdify(t, x, modules=['numpy'])

x_vals = linspace(0, 10, 100)
y_vals = lam_x(x_vals)

mpl.plot(x_vals, y_vals)
mpl.ylabel("Speed")
mpl.show()

(improvements suggested by asmeurer and MaxNoe)

enter image description here

Alternatively, you can use sympy's plot():

from sympy import symbols
from sympy import plot

t = symbols('t')
x = 0.05*t + 0.2/((t - 5)**2 + 2)

plot(x, (t, 0, 10), ylabel='Speed')