#!/usr/bin/env python
#version from November 11, 2007. Altered to be compatible vith vplotsvg
import sys, string
import vplot
infilename='lorenz.dat'
inf=open(infilename,'r')
n=0
a=vplot.eps_class()
a.scale(xmin=-15,xmax=15,ymin=0,ymax=30)
a.xaxis(dx=5)
a.yaxis(dy=5)
a.linewidth(10L)
print "reading lines from "+infilename+" ...be patient"
while 1:
	aline=inf.readline()
	if not aline: break
	if aline[0]=='#': continue
	t,x,y,z=map(eval,string.split(aline))
	s=.2*y
	red =max(0.,min(1.,.5+s))
	green =0.
	blue=max(0.,min(1.,.5-s))
	if n==0:
#		a.moveto(x,z)
		xl,zl=x,z
		n=n+1
		continue
	a.color(red,green,blue)
#	a.lineto(x,z) #works in postscript only
	a.draw((xl,zl,x,z)) #works for both postscript and svg
	xl,zl=x,z
	n=n+1
a.color(0,0,0)
a.text(.3j,.95j,0.,24,'Lorenz attractor')
a.text(.02j,.57j,0.,14,'z')
a.text(.5j,.02j,0.,14,'x')
a.close()