2005
Jul
11th
#!/usr/bin/env python
import Gnuplot
gp = ''
def init_plotter():
global gp
gp = Gnuplot.Gnuplot()
gp('set terminal png size 400, 240')
gp('set output "chaos.png"')
gp('set title "logistic map"')
gp('set data style lines')
gp('set xtics 10')
gp('set grid')
gp('set key left bottom')
def make_data(b):
v = [b]
while len(v) < 31:
next = 3.9 * v[-1] * (1 - v[-1])
v.append(next)
d = Gnuplot.Data(v, title = ('%f' % b), with = 'lines')
return d
if __name__ == '__main__':
init_plotter()
d1 = make_data(0.9)
d2 = make_data(0.91)
gp.plot(d1, d2)
Posted by setomits at 10:24 |
Comments: 0
調子に乗ってもうちょい。