y' = x(ρ−z)−y
z' = xy−βz
where σ=10, ρ=28 and β=8/3 using Euler's method.
The output of the program was graphed with gnuplot using the commands
$ gnuplot
G N U P L O T
Version 4.6 patchlevel 0 last modified 2012-03-04
Build System: Linux x86_64
Copyright (C) 1986-1993, 1998, 2004, 2007-2012
Thomas Williams, Colin Kelley and many others
gnuplot home: http://www.gnuplot.info
faq, bugs, etc: type "help FAQ"
immediate help: type "help" (plot window: hit 'h')
Terminal type set to 'x11'
gnuplot> set style data lines
gnuplot> plot "output.dat" using 2:3:4
gnuplot> quit
A screen shot of the resulting 3D graph is
Publication quality graphs can be made with the gnuplot script plot.plt using the commands
$ gnuplot plot.plt $ epstopdf --autorotate=All plot.eps