plot.py

import numpy
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt

with open("lorenz.dat") as fd:
    N, = numpy.fromfile(fd, numpy.dtype("i8"), 1)
    dtype = numpy.dtype("{N}f8,({N},3)f8".format(N=N))
    (time,state), = numpy.fromfile(fd, dtype, 1)
    x,y,z = [ state[:,i] for i in range(3) ]
fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')
ax.text2D(0.05, 0.95, "Lorenz attractor", transform=ax.transAxes)

ax.set_xlabel('x')
ax.set_ylabel('y')
ax.set_zlabel('z')

ax.plot(x,y,z)
plt.show()