%Sample run of y=double1(x) to see if its correct
echo on  %This command echos each line to the screen

x=linspace(0,1,200);
y=double1(x);
plot(x,y,'.',x,x)
title('The doubling function');

pause   %This command pauses the script until a key
        % is pressed
%What does the second iterate look like?
y1=double1(y);
plot(x,y1,'.',x,x);
title('The second iterate of D(x)');
pause

y2=double1(y1);
plot(x,y2,'.',x,x);
title('The third iterate of D(x)');

echo off