* SAS code to fit transformed data from Table 7.2;
 options linesize=70;

 data table72;
 input x y @@;
    v = 5*((y**0.2) - 1);    * transform Y observations;
    N0tr = 10**1.6;          * transform N0 = 10^8;
 cards;
  0.000   20   0.000  31   0.000  27
  0.002  100   0.002  92   0.002  89
  0.007  147   0.007 123   0.007 178
  0.014  216   0.014 170   0.014 181
  0.020  176   0.020 154   0.020 183
  0.030  154   0.030 153   0.030 149

 proc nlin method=dud ;
 parameters   b0 =  0.00000026
              b1 =  0.00034
              b3 =  89.312 ;
 model v = (5*N0tr*exp(-b3*x)*((1-exp(-b0-b1*x))**.2)) - 5;

 run;