SF.net SVN: matplotlib:[7978] trunk/py4science/examples/logistic/exercise02 .py

[<fer...@us...>]

Revision: 7978
          http://matplotlib.svn.sourceforge.net/matplotlib/?rev=7978&view=rev
Author:   fer_perez
Date:     2009-11-21 00:30:37 +0000 (Sat, 21 Nov 2009)

Log Message:
-----------
Make invariant density histogram clearer.

Modified Paths:
--------------
    trunk/py4science/examples/logistic/exercise02.py

Modified: trunk/py4science/examples/logistic/exercise02.py
===================================================================
--- trunk/py4science/examples/logistic/exercise02.py	2009-11-21 00:26:41 UTC (rev 7977)
+++ trunk/py4science/examples/logistic/exercise02.py	2009-11-21 00:30:37 UTC (rev 7978)
@@ -66,9 +66,8 @@
     return f
 
 
-def invariant_density(mu, x0,cycles=1000000,ret_all=False):
+def invariant_density(mu, x0,cycles=1000000,ret_all=False, bins=500):
     transients = 1000
-    bins = 500
     f = plt.figure()
     ax = f.add_subplot(111)
     logmap = Logistic(mu)
@@ -94,7 +93,7 @@
         return 1./(np.pi * np.sqrt( x*(1.-x)))
 
     # Don't start from 0.5, which is a fixed point!
-    f = invariant_density(1.0,0.567)
+    f = invariant_density(1.0, 0.567, bins=100)
     ax = f.gca()
     # avoid the edges: rho(x) is singular at 0 and 1!
     x0 = np.linspace(0.001, 0.999, 1000)


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.