|
From: <fer...@us...> - 2009-11-21 00:30:44
|
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.
|