cvode-rust-wrap/examples/plot_oscillator.py
2021-06-10 19:37:09 +02:00

8 lines
299 B
Python

import pandas as pd
import sys
import matplotlib.pyplot as plt
df = pd.read_csv(sys.stdin,names=['t','x',r'\dot{x}',r"dx_dx0", r"d\dot{x}_dx0", r"dx_d\dot{x}0", r"d\dot{x}_d\dot{x}0", r"dx_dk", r"d\dot{x}_dk"],index_col='t')
ax = df.plot(subplots=True)
plt.suptitle(r"\dotdot{x} = -k*x")
plt.show()