ClearAll;
$DefaultFont = {"Times-Roman",10};
Needs["Graphics`Graphics`"];

xticks = {0, 0.25, 0.5, 0.75, 1};
yticks = {1, 1.2, 1.4, 1.6};
ph[1,x] = 1;
ph[n_,x_]:= Together[1 + x/ph[n-1,x]];
plotsq:= Plot[(1+Sqrt[1+4x])/2,{x,0,1},
			PlotStyle -> Dashing[{0.01,0.01}]];
plot[n_]:= Plot[Evaluate[ph[n,x]],{x,0,1},
			    Ticks -> {xticks, yticks}];
a11:= Show[plot[1],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 1", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
a12:= Show[plot[2],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 2", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
a21:= Show[plot[3],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 3", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
a22:= Show[plot[4],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 4", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
a31:= Show[plot[5],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 5", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
a32:= Show[plot[6],plotsq, AxesOrigin -> {-0.1,0.9}, Axes->True,
		   PlotLabel -> "n = 6", AxesLabel -> {"x","y"},
		   PlotRange -> {{0,1},{0.95,1.65}}];
Show[GraphicsArray[{{a11,a12},{a21,a22},{a31,a32}}],
		PlotLabel->"                                                    Figure 7.1\n\n\n
		     Oscillatory convergence of Fibonacci rational function sequence"]