StaticODESolver-LorenzExample.py

#!/usr/bin/env python3

import sys
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import axes3d, Axes3D
import numpy as np

plt.rcParams['text.usetex'] = True

f = open( sys.argv[1], 'r' )
x_lst = []
y_lst = []
z_lst = []
for line in f:
    line = line.strip()
    a = line.split()
    x_lst.append( float( a[ 0 ] ) )
    y_lst.append( float( a[ 1 ] ) )
    z_lst.append( float( a[ 2 ] ) )

x = np.array(x_lst)
y = np.array(y_lst)
z = np.array(z_lst)

fig = plt.figure()
ax = Axes3D(fig)
theta = np.linspace(-4 * np.pi, 4 * np.pi, 100)
ax.plot(x, y, z, label='Lorenz attractor')
ax.legend()

plt.savefig( sys.argv[2] )
plt.close(fig)