%Script to look at fourth iterate

x=linspace(-2,2,100);
c=linspace(1/4,-1.7,80);
for j=1:length(c)
    fprintf('Value of c is: %f\n',c(j));
    y1=x.^2+c(j);
    y2=y1.^2+c(j);
    y3=y2.^2+c(j);
    y4=y3.^2+c(j);
    figure(1)
    plot(x,x,x,y1,'y-',x,y2,'g-');
    axis([-2 2 -2 2]);
    title('Plot of Q_c and Q_c^2')
    figure(2)
    plot(x,x,x,y4,'k-',x,y2,'g-',x,y1,'y-');
    axis([-2,2,-2,2]);
    pause
end