Mercurial > repos > rliterman > csp2

annotate CSP2/CSP2_env/env-d9b9114564458d9d-741b3de822f2aaca6c6caa4325c4afce/lib/python3.8/turtledemo/chaos.py @ 69:33d812a61356

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 2e9511a184a1ca667c7be0c6321a36dc4e3d116d
author jpayne
date Tue, 18 Mar 2025 17:55:14 -0400
parents
children
rev   line source
jpayne@69 1 # File: tdemo_chaos.py
jpayne@69 2 # Author: Gregor Lingl
jpayne@69 3 # Date: 2009-06-24
jpayne@69 4
jpayne@69 5 # A demonstration of chaos
jpayne@69 6
jpayne@69 7 from turtle import *
jpayne@69 8
jpayne@69 9 N = 80
jpayne@69 10
jpayne@69 11 def f(x):
jpayne@69 12 return 3.9*x*(1-x)
jpayne@69 13
jpayne@69 14 def g(x):
jpayne@69 15 return 3.9*(x-x**2)
jpayne@69 16
jpayne@69 17 def h(x):
jpayne@69 18 return 3.9*x-3.9*x*x
jpayne@69 19
jpayne@69 20 def jumpto(x, y):
jpayne@69 21 penup(); goto(x,y)
jpayne@69 22
jpayne@69 23 def line(x1, y1, x2, y2):
jpayne@69 24 jumpto(x1, y1)
jpayne@69 25 pendown()
jpayne@69 26 goto(x2, y2)
jpayne@69 27
jpayne@69 28 def coosys():
jpayne@69 29 line(-1, 0, N+1, 0)
jpayne@69 30 line(0, -0.1, 0, 1.1)
jpayne@69 31
jpayne@69 32 def plot(fun, start, color):
jpayne@69 33 pencolor(color)
jpayne@69 34 x = start
jpayne@69 35 jumpto(0, x)
jpayne@69 36 pendown()
jpayne@69 37 dot(5)
jpayne@69 38 for i in range(N):
jpayne@69 39 x=fun(x)
jpayne@69 40 goto(i+1,x)
jpayne@69 41 dot(5)
jpayne@69 42
jpayne@69 43 def main():
jpayne@69 44 reset()
jpayne@69 45 setworldcoordinates(-1.0,-0.1, N+1, 1.1)
jpayne@69 46 speed(0)
jpayne@69 47 hideturtle()
jpayne@69 48 coosys()
jpayne@69 49 plot(f, 0.35, "blue")
jpayne@69 50 plot(g, 0.35, "green")
jpayne@69 51 plot(h, 0.35, "red")
jpayne@69 52 # Now zoom in:
jpayne@69 53 for s in range(100):
jpayne@69 54 setworldcoordinates(0.5*s,-0.1, N+1, 1.1)
jpayne@69 55 return "Done!"
jpayne@69 56
jpayne@69 57 if __name__ == "__main__":
jpayne@69 58 main()
jpayne@69 59 mainloop()