(*
	file:   Cubic.m
	author: Hinke Osinga, The Geometry Center
	data:   April 15, 1998
*)


Print["The Cubic Map:"];
Print["f[x] = lambda (x - x^3/3)"];

DynSys[x_] := N[lambda (x - x^3/3)];

x0 = 0.123;
lambda = 2.25;
its = 100;
xrange = {-2, 2};

Print["lambda = ", lambda, ", x0 = ", x0, ", its = ", its];
Print[""];

critical = 1;
lrange = {0, 3Sqrt[3]/2};