-- 60-


where


          Ce (n+l) = value of Ce after n + 1 iterations

            Ce (n) = value of Ce after n iterations

            C (i) = initial estimate of C e                                       e

                  dGC
            e         , e C = C (n)
                           e  e

G(Ce(n)) is obtained directly from (4.20). G' (Ce(n)) is obtained by differentiation of (4.20).


              13       ~1i-~
G'(Ce) = 1 +  I   (2-i) (l+Zi)C    - (1-i)Z. C  K'
    e'       i=2 L                                 i


                                                   (4.22)


Using equations (4.20) and (4.22) with (4.21), the electron fraction is then given by C (n+l) = Ce (n)


      Ca(n) +   3          2l+i)C (n) - Z   (n) Kii=2
        13
   1 +   n) I 2-i) (l+zi)Ci(n) -      iZ C e-i(n)K K'
       i=                                          i(4.23)

Beginning with an assumed average electron concentration of Ce (1) = 0.5 as the initial estimate, expression (4.23)