%Sample script to graphically find periodic points.
%
%Step 1:  Create a domain for the graph.
xinit=linspace(-1,3,500);  %Creates 500 evenly spaced points
                           %  between -1 and 3.
x=xinit;                       
y1=-(3/2)*x.^2 + (5/2)*x + 1;
x=y1;                       
y1=-(3/2)*x.^2 + (5/2)*x + 1;
x=y1;
y1=-(3/2)*x.^2 + (5/2)*x + 1;

plot(xinit,y1,xinit,xinit);
axis([-1,3,-1.5,2.3]);
title('Third iterate for -3/2x^2+5/2x+1')