# -*- coding: utf-8 -*-
"""
logistic
Created on Tue Sep 24 14:48:52 2013

@author: coriander

$ ipython
>>> from pylab import *
>>> from logistic import *
>>> plot (xs(3.5))

"""

from pylab import *

def logistic(a, x):
    value = a*x*(1-x)
    return value

def list_xs(a, x0 = 0.1, n=200):
    ''' returns a list of x's'''
    x_list = [x0]
    for i in xrange(n):
        x_val = logistic(a, x_list[-1])
        #print x_val
        x_list.append(x_val)
    for i in xrange(n-16, len(x_list)):
        print "{}: {:.5f}".format(i,x_list[i])
    return x_list

def flt_step_xrange(start, stop, step):
    # From bottom of https://wiki.python.org/moin/ForLoop
    step = float(step)
    while start <= stop:
        yield start # if I understand this right, this is return without an exit
        start += step

list_xs(3.5) 
print "3.56"
list_xs(3.56) 
print "3.57"
list_xs(3.57) 
print "3.58"
list_xs(3.58)    
print "-------"
         
        
#def main():
#    #the_range = (3.48, 3.59)
#    #num_steps = 100*(the_range[1] - the_range[0])
#    a_range = (3.48, 3.59)
#    t_range = (3.5, 3.7)
#    for i in flt_step_xrange(t_range[0], t_range[1], 0.01):
#        print "{:.5f}".format(list_xs(a=i))
#        
#        
#
#main()