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

printf("The Quartic Map: ");
printf("f(x) = lambda x (1 - x) (1 - 2x)^2\n");

DynSys := x -> evalf(lambda * x * (1 - x) * (1 - 2*x)^2):

x0 := 0.94:
lambda := 5.8:
xrange := 0..1:
lrange := [0, 16]:
critical := 0.5 - 0.25*evalf(sqrt(2)):

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