Sage Interactions - Dynamical Systems

goto interact main page

Cobweb diagrams on [0,1]

by Marshall Hampton.

cobweb.png

Cythonized Logistic Orbit Map

By Marshall Hampton

%cython
cpdef double logorb(double k,long N,double x0):
    cdef double x = x0
    cdef long i 
    for i from 1 <= i <= N:
        x = k*x*(1-x)
    return x

cpdef logtraj(double k,long N, double x0):
    cdef double x = x0
    xvals = []
    cdef long i 
    for i from 1 <= i <= N:
        x = k*x*(1-x)
        xvals.append(x)
    return xvals

logisticorbits.png