#	file:   Cubic.mpl
#	author: Hinke Osinga, University of Bristol
#	date:   January 2001
#
# 

printf("The Cubic Map: ");
printf("f(x) = lambda (x - x^3/3)\n");

DynSys := x -> evalf(lambda * (x - x^3 / 3)):

x0 := 0.94:
lambda := 2.25:
xrange := -2..2:
lrange := [0, 3*evalf(sqrt(3))/2]:
critical := 1:

printf("lambda = %f, x0 = %f\n\n", lambda, x0);