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


Print["The Quartic Map:"];
Print["f[x] = lambda x (1 - x) (1 - 2x)^2"];

DynSys[x_] := N[lambda x (1 - x) (1 - 2x)^2];

x0 = 0.94;
lambda = 5.8;
its = 100;
xrange = {0, 1};

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

critical = 0.5 - 0.25 Sqrt[2];
lrange = {0, 16};