ClearAll;Clear[f1];Clear[f2];
$DefaultFont = {"Times-Roman",10};
Needs["Statistics`DataManipulation`"];
Needs["Graphics`Graphics`"];
cumulative = {312/545,413/545,487/545,522/545,537/545,1};
datapoints = {{0,0},{1,312/545},{2,413/545},
			  {3,487/545},{4,522/545},{5,537/545},{6,1}};
density = {312/545,101/545,74/545,35/545,15/545,8/545};
F1[A_,c_,x_] = 4 A c x/(1 + A c^2) - A x^2;
F2[A_,c_,x_] = 1 - ((1 - A c^2)^2 (c/x)^(2A c^2/(1 - A c^2)))/(1 + A c^2);

f1[A_,c_,x_] = 2 A(c-x) + 2 A c(1 - A c^2)/(1 + A c^2);
f2[A_,c_,x_] = (2A c(1 - A c^2) (c/x)^((1 + A c^2)/(1 - A c^2)))/(1 + A c^2);

FF[a_,c_,x_] = Which[x<=1.5,F1[a,c,x],x>1.5,F2[a,c,x]];

error[a_,c_] = Evaluate[Apply[Plus,Table[(FF[a,c,n]-cumulative[[n]])^2,{n,1,6}]]];

minimizer[c_]:= Max[a/.FindMinimum[error[a,c],{a,0.13,0.131,0,0.15}][[2]],0.12];
leasterror[c_]:= error[minimizer[c],c];
				
topleft = Show[Plot[error[a,1.7],{a,0.12,0.15},PlotRange->{{0.119,0.151},{0.006,0.035}},
			   AxesOrigin->{0.118,0.005},AxesLabel->{"  A","z"},
			   Ticks->{{0.12,0.13,0.14,0.15},{0.01,0.02,0.03}},
			   PlotLabel->"(a) z  =  S(A, 1.7)"],
			   Graphics[Disk[{minimizer[1.7],leasterror[1.7]},{0.0004,0.0005}]]];
			   
topright = Show[Plot[error[a,1.8],{a,0.12,0.15},PlotRange->{{0.119,0.151},{0.006,0.035}},
			   AxesOrigin->{0.118,0.005},AxesLabel->{"  A","z"},
			   Ticks->{{0.12,0.13,0.14,0.15},{0.01,0.02,0.03}},
			   PlotLabel->"(b) z  =  S(A, 1.8)"],
			   Graphics[Disk[{minimizer[1.8],leasterror[1.8]},{0.0004,0.0005}]]];
			   
bottomleft = Show[Plot[error[a,1.9],{a,0.12,0.15},PlotRange->{{0.119,0.151},{0.006,0.035}},
			   AxesOrigin->{0.118,0.005},AxesLabel->{"  A","z"},
			   Ticks->{{0.12,0.13,0.14,0.15},{0.01,0.02,0.03}},
			   PlotLabel->"(c) z  =  S(A, 1.9)"],
			   Graphics[Disk[{minimizer[1.9],leasterror[1.9]},{0.0004,0.0005}]]];
			   
bottomright = Show[Plot[error[a,2],{a,0.12,0.15},PlotRange->{{0.119,0.151},{0.006,0.035}},
			   AxesOrigin->{0.118,0.005},AxesLabel->{"  A","z"},
			   Ticks->{{0.12,0.13,0.14,0.15},{0.01,0.02,0.03}},
			   PlotLabel->"(d) z  =  S(A, 2)"],
			   Graphics[Disk[{minimizer[2],leasterror[2]},{0.0004,0.0005}]]];
			   
Show[GraphicsArray[{{topleft,topright},{bottomleft,bottomright}}],
				  PlotLabel->"Figure 24.5"]