Page 412
% make an m-file, retinox0.m, with
%   function x0=retinox0(p)
%   for i=1:22
%   x0(i)=(1-p)^(2^i-1);
%   end
% make an m-file, retinox1.m, with
%   function x1=retinox1(p,x0)
%   x1(1)=p;
%   for i=2:22
%      x1(i)=x0(i-1)*(2*(1-p)*x1(i-1)+p);
%   end
% make an m-file, retinox2.m, with
%   function x2=retinox2(p,x0,x1)
%   x2(1)=0;
%   for i=2:22
%      x2(i)=(1-p)*(2*x0(i-1)*x2(i-1)+...
%         x1(i-1)^2)+p*x1(i-1);
%   end
% make an m-file, retinox3.m, with
%   function x3=retinox3(p,x0,x1,x2)
%   x3(1)=0;
%   for i=2:22
%      x3(i)=2*(1-p)*(x1(i-1)*x2(i-1)+...
%         x0(i-1)*x3(i-1))+p*x2(i-1);
%   end
   p=.0000007;
   x0=retinox0(p)
   x1=retinox1(p,x0)
   x2=retinox2(p,x0,x1)
   x3=retinox3(p,x0,x1,x2)