      function txfit       (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff(  4)
      data ncoeff/  3/
      data avdat/ -0.2957127E-02/
      data xmin/
     1 -0.57695E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.50054E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.33542139E-02, 0.89052327E-01,-0.11083363E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
c          set up monomials   functions
      x11 = x1
      x12 = x11*x1
c
c                  function
c
      txfit       =avdat
     1  +coeff(  1)    
     2  +coeff(  2)*x11
     3  +coeff(  3)*x12
c
      return
      end
      function delta       (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 51)
      data ncoeff/ 50/
      data avdat/ -0.1282164E-02/
      data xmin/
     1 -0.69571E+00,-0.28888E-01,-0.42594E-01,-0.42817E-01,-0.99954E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.54342E+00, 0.24888E-01, 0.42594E-01, 0.42817E-01, 0.99924E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.49380390E-02, 0.50776087E-01, 0.41714744E-02,-0.16485738E-02,
     +  0.40860069E-02,-0.13720992E-02, 0.76295994E-03, 0.19616097E-03,
     + -0.72629337E-03,-0.11135841E-02, 0.98215812E-03, 0.10337732E-03,
     +  0.19684956E-03, 0.38250367E-03, 0.28657675E-03,-0.57894050E-03,
     +  0.66203420E-03, 0.63566654E-03,-0.92776940E-03, 0.68306108E-03,
     + -0.39594181E-03,-0.32945023E-04,-0.10452494E-03, 0.46513884E-04,
     + -0.17427166E-03,-0.72052964E-03,-0.52282802E-03,-0.74086856E-03,
     +  0.89874033E-04, 0.45724670E-03,-0.46269043E-03, 0.14136000E-03,
     + -0.27232111E-03,-0.34566221E-03, 0.25535608E-03,-0.68936625E-03,
     + -0.64657157E-03,-0.64245665E-04, 0.23122663E-03,-0.22399647E-03,
     +  0.11333229E-03, 0.13813306E-04, 0.13819366E-03,-0.45124142E-03,
     +  0.20185210E-03, 0.32523990E-03,-0.42648180E-03,-0.43995929E-03,
     + -0.34055795E-03,-0.76545676E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x12 = x11*x1
      x13 = x12*x1
      x14 = x13*x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
c
c                  function
c
      delta       =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)*x11                
     3  +coeff(  3)*x12                
     4  +coeff(  4)    *x21            
     5  +coeff(  5)*x11*x21            
     6  +coeff(  6)    *x22            
     7  +coeff(  7)        *x32        
     8  +coeff(  8)                *x51
      delta       =delta       
     9  +coeff(  9)    *x21*x31*x41    
     1  +coeff( 10)    *x21    *x42    
     2  +coeff( 11)        *x31*x41    
     3  +coeff( 12)            *x42    
     4  +coeff( 13)*x13                
     5  +coeff( 14)*x12*x21            
     6  +coeff( 15)*x11*x22            
     7  +coeff( 16)    *x21*x32        
     8  +coeff( 17)*x11    *x31*x41    
      delta       =delta       
     9  +coeff( 18)    *x24            
     1  +coeff( 19)    *x22    *x42    
     2  +coeff( 20)*x11        *x42    
     3  +coeff( 21)*x11*x21*x31*x41    
     4  +coeff( 22)    *x23            
     5  +coeff( 23)*x11    *x32        
     6  +coeff( 24)    *x21        *x51
     7  +coeff( 25)*x11*x23            
     8  +coeff( 26)    *x22*x32        
      delta       =delta       
     9  +coeff( 27)*x11*x22*x31*x41    
     1  +coeff( 28)*x11*x22    *x42    
     2  +coeff( 29)*x14*x22            
     3  +coeff( 30)    *x24    *x42    
     4  +coeff( 31)    *x23*x34        
     5  +coeff( 32)*x12*x22            
     6  +coeff( 33)    *x22*x31*x41    
     7  +coeff( 34)*x12*x23            
     8  +coeff( 35)*x11*x24            
      delta       =delta       
     9  +coeff( 36)    *x21*x33*x41    
     1  +coeff( 37)*x11*x21*x31*x43    
     2  +coeff( 38)*x13*x21            
     3  +coeff( 39)*x11*x21*x32        
     4  +coeff( 40)*x12    *x31*x41    
     5  +coeff( 41)*x12        *x42    
     6  +coeff( 42)    *x22        *x51
     7  +coeff( 43)*x13*x22            
     8  +coeff( 44)    *x23*x32        
      delta       =delta       
     9  +coeff( 45)*x12*x21*x31*x41    
     1  +coeff( 46)    *x23    *x42    
     2  +coeff( 47)    *x21*x32*x42    
     3  +coeff( 48)*x11*x23*x31*x41    
     4  +coeff( 49)*x11*x23    *x42    
     5  +coeff( 50)*x11*x21*x32*x42    
c
      return
      end
      function theta       (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 61)
      data ncoeff/ 60/
      data avdat/  0.1547449E-02/
      data xmin/
     1 -0.69571E+00,-0.28888E-01,-0.42594E-01,-0.42817E-01,-0.99954E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.54342E+00, 0.24888E-01, 0.42594E-01, 0.42817E-01, 0.99924E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.32031753E-02,-0.52152958E-03,-0.63572347E-01, 0.86818938E-04,
     +  0.84566902E-02,-0.49368520E-02, 0.62443910E-03, 0.19738826E-02,
     + -0.11215348E-02, 0.17805009E-03, 0.72298275E-03,-0.90115861E-03,
     + -0.41810339E-03,-0.21269202E-03,-0.17416036E-03, 0.87057572E-03,
     + -0.18210530E-02,-0.10068862E-02, 0.26722773E-03,-0.41330786E-03,
     + -0.90603440E-04,-0.36462740E-03,-0.44847617E-03, 0.14445485E-02,
     + -0.92949427E-03,-0.10356495E-03,-0.12155218E-02,-0.97944438E-04,
     +  0.18603702E-03, 0.53253165E-03, 0.25406176E-04,-0.20909202E-02,
     + -0.48517045E-02,-0.49993483E-03, 0.78767986E-03,-0.68630646E-02,
     + -0.12530059E-02, 0.33226880E-03, 0.34881558E-03, 0.12953716E-03,
     +  0.15649994E-03, 0.80768816E-03,-0.40171237E-03, 0.67872490E-04,
     +  0.46987530E-04,-0.49193110E-04,-0.11205965E-02,-0.26063268E-04,
     +  0.28118910E-04,-0.92174613E-03,-0.55402652E-04,-0.91498229E-03,
     +  0.15436711E-02,-0.13693257E-02,-0.41131405E-03, 0.17153226E-03,
     +  0.11307178E-03, 0.28696402E-02, 0.32500261E-02, 0.21205033E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x12 = x11*x1
      x13 = x12*x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x51 = x5
c
c                  function
c
      theta       =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)*x11                
     3  +coeff(  3)    *x21            
     4  +coeff(  4)*x12                
     5  +coeff(  5)*x11*x21            
     6  +coeff(  6)    *x22            
     7  +coeff(  7)                *x51
     8  +coeff(  8)*x11*x22            
      theta       =theta       
     9  +coeff(  9)    *x21*x32        
     1  +coeff( 10)        *x31*x41    
     2  +coeff( 11)            *x42    
     3  +coeff( 12)    *x23            
     4  +coeff( 13)    *x21    *x42    
     5  +coeff( 14)*x13                
     6  +coeff( 15)*x11        *x42    
     7  +coeff( 16)*x11*x23            
     8  +coeff( 17)*x11*x21*x31*x41    
      theta       =theta       
     9  +coeff( 18)    *x23*x31*x41    
     1  +coeff( 19)*x12*x21    *x42    
     2  +coeff( 20)*x12*x21            
     3  +coeff( 21)*x11    *x32        
     4  +coeff( 22)    *x21*x31*x41    
     5  +coeff( 23)*x12*x22            
     6  +coeff( 24)*x11*x21*x32        
     7  +coeff( 25)    *x22*x32        
     8  +coeff( 26)*x11            *x51
      theta       =theta       
     9  +coeff( 27)*x11*x23    *x42    
     1  +coeff( 28)*x12        *x42    
     2  +coeff( 29)*x11*x21    *x42    
     3  +coeff( 30)    *x22    *x42    
     4  +coeff( 31)    *x21        *x51
     5  +coeff( 32)    *x23*x32        
     6  +coeff( 33)    *x21*x33*x41    
     7  +coeff( 34)*x11*x22    *x42    
     8  +coeff( 35)*x11*x21*x33*x41    
      theta       =theta       
     9  +coeff( 36)*x11*x21*x32*x42    
     1  +coeff( 37)*x11*x22*x31*x43    
     2  +coeff( 38)        *x32        
     3  +coeff( 39)*x11    *x31*x41    
     4  +coeff( 40)*x12    *x31*x41    
     5  +coeff( 41)    *x22*x31*x41    
     6  +coeff( 42)*x12*x21*x31*x41    
     7  +coeff( 43)*x13        *x42    
     8  +coeff( 44)    *x22        *x51
      theta       =theta       
     9  +coeff( 45)        *x32    *x51
     1  +coeff( 46)*x11*x23        *x51
     2  +coeff( 47)*x11*x22*x31*x41    
     3  +coeff( 48)    *x21*x32*x42    
     4  +coeff( 49)        *x31*x41*x51
     5  +coeff( 50)    *x22*x33*x41    
     6  +coeff( 51)*x11*x22        *x51
     7  +coeff( 52)*x12*x23*x32        
     8  +coeff( 53)*x11*x22*x33*x41    
      theta       =theta       
     9  +coeff( 54)*x12*x21*x31*x43    
     1  +coeff( 55)*x13*x23*x31*x41    
     2  +coeff( 56)    *x23*x32    *x51
     3  +coeff( 57)*x11*x22    *x42*x51
     4  +coeff( 58)*x12*x23*x33*x41    
     5  +coeff( 59)*x12*x23*x32*x42    
     6  +coeff( 60)        *x33*x41    
c
      return
      end
      function phi         (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 61)
      data ncoeff/ 60/
      data avdat/  0.0000000E+00/
      data xmin/
     1 -0.69571E+00,-0.28888E-01,-0.42594E-01,-0.42817E-01,-0.99954E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.54342E+00, 0.24888E-01, 0.42594E-01, 0.42817E-01, 0.99924E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.13402536E-01,-0.27650431E-01, 0.82578193E-02, 0.45992858E-02,
     + -0.40542646E-02, 0.71242359E-02, 0.43149018E-02,-0.32907617E-02,
     + -0.34422223E-02,-0.12285090E-02,-0.96058490E-03, 0.17237580E-03,
     + -0.63819590E-03,-0.12440885E-02,-0.53748151E-03, 0.54805703E-03,
     + -0.81267784E-03, 0.13568046E-02,-0.10177636E-02,-0.21962095E-02,
     +  0.20590706E-02, 0.24526690E-02, 0.23862300E-02,-0.17261941E-02,
     +  0.68260604E-04, 0.91567769E-03,-0.59411541E-03, 0.47777031E-03,
     +  0.45710840E-03, 0.94540053E-04, 0.61446073E-04,-0.17115258E-02,
     + -0.36682002E-03, 0.41913049E-03,-0.29439939E-03, 0.12338144E-03,
     +  0.50812092E-03,-0.41950928E-03,-0.54420455E-03, 0.27981692E-03,
     + -0.25018239E-02,-0.14759736E-02, 0.25802916E-04,-0.54292358E-03,
     +  0.78980118E-03, 0.38464030E-03, 0.58107089E-04, 0.13236063E-05,
     + -0.55479775E-02, 0.16117278E-02,-0.15015166E-02,-0.84824208E-03,
     +  0.66479953E-03,-0.61519574E-02,-0.38071294E-03,-0.40244058E-03,
     + -0.50501147E-03, 0.60866791E-03,-0.50770090E-03, 0.49054885E-03,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x12 = x11*x1
      x13 = x12*x1
      x14 = x13*x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x45 = x44*x4
      x51 = x5
c
c                  function
c
      phi         =avdat
     1  +coeff(  1)        *x31        
     2  +coeff(  2)            *x41    
     3  +coeff(  3)*x11    *x31        
     4  +coeff(  4)    *x21*x31        
     5  +coeff(  5)*x11        *x41    
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)*x12        *x41    
     8  +coeff(  8)    *x22*x31        
      phi         =phi         
     9  +coeff(  9)*x11*x22    *x41    
     1  +coeff( 10)        *x32*x41    
     2  +coeff( 11)*x12    *x31        
     3  +coeff( 12)*x11*x21    *x41    
     4  +coeff( 13)*x11    *x32*x41    
     5  +coeff( 14)*x11    *x31*x42    
     6  +coeff( 15)        *x33        
     7  +coeff( 16)*x11*x21*x31        
     8  +coeff( 17)*x12*x21*x31        
      phi         =phi         
     9  +coeff( 18)*x11*x22*x31        
     1  +coeff( 19)    *x23*x31        
     2  +coeff( 20)*x13        *x41    
     3  +coeff( 21)*x12*x21    *x41    
     4  +coeff( 22)*x14        *x41    
     5  +coeff( 23)*x12*x22    *x41    
     6  +coeff( 24)*x13*x21    *x43    
     7  +coeff( 25)*x13*x23    *x41    
     8  +coeff( 26)*x13*x22    *x43    
      phi         =phi         
     9  +coeff( 27)        *x31*x42    
     1  +coeff( 28)*x13    *x31        
     2  +coeff( 29)    *x23    *x41    
     3  +coeff( 30)            *x45    
     4  +coeff( 31)    *x22*x33        
     5  +coeff( 32)*x11*x22*x31*x42    
     6  +coeff( 33)    *x22    *x41    
     7  +coeff( 34)*x11    *x33        
     8  +coeff( 35)    *x21*x33        
      phi         =phi         
     9  +coeff( 36)    *x21*x31*x42    
     1  +coeff( 37)*x12    *x32*x41    
     2  +coeff( 38)*x11*x21*x32*x41    
     3  +coeff( 39)    *x22*x32*x41    
     4  +coeff( 40)    *x22*x31*x42    
     5  +coeff( 41)*x13*x21    *x41    
     6  +coeff( 42)*x11*x23    *x41    
     7  +coeff( 43)        *x31    *x51
     8  +coeff( 44)*x11*x22*x32*x41    
      phi         =phi         
     9  +coeff( 45)*x13        *x43    
     1  +coeff( 46)*x14*x21*x31        
     2  +coeff( 47)    *x21*x31    *x51
     3  +coeff( 48)*x12    *x34*x41    
     4  +coeff( 49)    *x22*x33*x42    
     5  +coeff( 50)    *x22*x32*x43    
     6  +coeff( 51)    *x24*x33        
     7  +coeff( 52)*x11*x23    *x43    
     8  +coeff( 53)    *x24    *x43    
      phi         =phi         
     9  +coeff( 54)*x11*x24*x32*x43    
     1  +coeff( 55)*x11        *x43    
     2  +coeff( 56)*x11*x21*x31*x42    
     3  +coeff( 57)*x14    *x31        
     4  +coeff( 58)*x13*x21*x31        
     5  +coeff( 59)*x12*x22*x31        
     6  +coeff( 60)*x11*x23*x31        
c
      return
      end
      function y00         (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 81)
      data ncoeff/ 80/
      data avdat/  0.0000000E+00/
      data xmin/
     1 -0.69571E+00,-0.28888E-01,-0.42594E-01,-0.42817E-01,-0.99954E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.54342E+00, 0.24888E-01, 0.42594E-01, 0.42817E-01, 0.99924E-03,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.48885833E-01, 0.35042211E-01,-0.18846972E-01,-0.26392015E-01,
     + -0.32463990E-01,-0.15710082E-01, 0.15118452E-02, 0.13637152E-02,
     +  0.90423375E-02,-0.91108903E-02,-0.15958060E-01, 0.19199982E-01,
     +  0.14078059E-02, 0.28932428E-02, 0.83648609E-02, 0.18956250E-02,
     + -0.62049564E-03, 0.32438783E-02, 0.20524459E-02, 0.10909408E-02,
     +  0.13056819E-02,-0.17379763E-02,-0.20552431E-02,-0.37553260E-03,
     +  0.92766166E-03, 0.13767019E-02, 0.57283719E-02, 0.22233014E-02,
     + -0.37635001E-02, 0.28313589E-02,-0.86744496E-03, 0.51456457E-03,
     + -0.13745476E-02,-0.10265069E-03, 0.68503488E-02, 0.37133126E-02,
     +  0.30031172E-03, 0.98112819E-03,-0.11317794E-02, 0.29400224E-02,
     +  0.18111353E-02, 0.65214554E-03,-0.17529501E-02,-0.14002695E-02,
     + -0.11329388E-02, 0.49193790E-02,-0.18903877E-02, 0.44999062E-03,
     + -0.22409034E-02, 0.42555346E-02, 0.33099856E-02, 0.16840221E-02,
     + -0.10843820E-03, 0.34134302E-02, 0.21102938E-02,-0.66344778E-03,
     +  0.14085282E-02, 0.75809518E-03,-0.29004863E-02,-0.47186748E-03,
     +  0.32748117E-02, 0.79825812E-03,-0.39457594E-03, 0.74100861E-03,
     + -0.26531448E-02,-0.84257108E-03, 0.66547317E-03, 0.12885126E-02,
     + -0.10343362E-02, 0.88993516E-02,-0.10669783E-02, 0.23255276E-02,
     + -0.11499356E-02,-0.31322809E-02,-0.13955836E-02, 0.18164153E-02,
     + -0.62069083E-02, 0.77746748E-02,-0.87935288E-04, 0.46355920E-02,
     +      0.      /
      data ientry/0/
c
      if (ientry.ne.0) go to 10
      ientry=1
      do 5 i=1,m
      if (xmin(i).eq.xmax(i)) go to 5
      scale(i)=2./(xmax(i)-xmin(i))
   5  continue
c
  10  continue
c      normalize variables between -1 and +1
      x1 =1.+(x(  1)-xmax(  1))*scale(  1)
      x2 =1.+(x(  2)-xmax(  2))*scale(  2)
      x3 =1.+(x(  3)-xmax(  3))*scale(  3)
      x4 =1.+(x(  4)-xmax(  4))*scale(  4)
      x5 =1.+(x(  5)-xmax(  5))*scale(  5)
c          set up monomials   functions
      x11 = x1
      x12 = x11*x1
      x13 = x12*x1
      x14 = x13*x1
      x15 = x14*x1
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*x2
      x25 = x24*x2
      x31 = x3
      x32 = x31*x3
      x33 = x32*x3
      x34 = x33*x3
      x35 = x34*x3
      x41 = x4
      x42 = x41*x4
      x43 = x42*x4
      x44 = x43*x4
      x51 = x5
c
c                  function
c
      y00         =avdat
     1  +coeff(  1)        *x31        
     2  +coeff(  2)            *x41    
     3  +coeff(  3)*x11    *x31        
     4  +coeff(  4)    *x21*x31        
     5  +coeff(  5)*x11        *x41    
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31*x42    
     8  +coeff(  8)*x11*x21*x31        
      y00         =y00         
     9  +coeff(  9)    *x22*x31        
     1  +coeff( 10)*x12        *x41    
     2  +coeff( 11)*x11*x21    *x41    
     3  +coeff( 12)    *x22    *x41    
     4  +coeff( 13)*x11        *x43    
     5  +coeff( 14)    *x23*x31        
     6  +coeff( 15)    *x23    *x41    
     7  +coeff( 16)        *x33        
     8  +coeff( 17)*x13        *x41    
      y00         =y00         
     9  +coeff( 18)*x15        *x41    
     1  +coeff( 19)*x12    *x31        
     2  +coeff( 20)*x11    *x32*x41    
     3  +coeff( 21)*x11    *x31*x42    
     4  +coeff( 22)*x11*x22*x33        
     5  +coeff( 23)*x15*x21    *x41    
     6  +coeff( 24)    *x21    *x43    
     7  +coeff( 25)*x13    *x31        
     8  +coeff( 26)*x11*x22    *x41    
      y00         =y00         
     9  +coeff( 27)    *x22*x32*x41    
     1  +coeff( 28)*x14        *x41    
     2  +coeff( 29)*x12*x22    *x41    
     3  +coeff( 30)*x11*x23    *x41    
     4  +coeff( 31)    *x21*x33*x42    
     5  +coeff( 32)*x13*x22*x31        
     6  +coeff( 33)*x11*x24*x31        
     7  +coeff( 34)            *x41*x51
     8  +coeff( 35)*x12*x22*x31*x42    
      y00         =y00         
     9  +coeff( 36)*x11*x24*x31*x42    
     1  +coeff( 37)*x14*x23*x31        
     2  +coeff( 38)            *x43    
     3  +coeff( 39)        *x32*x43    
     4  +coeff( 40)*x12        *x43    
     5  +coeff( 41)    *x22    *x43    
     6  +coeff( 42)*x12*x22*x31        
     7  +coeff( 43)*x11*x23*x31        
     8  +coeff( 44)    *x24    *x41    
      y00         =y00         
     9  +coeff( 45)    *x21*x35        
     1  +coeff( 46)*x11    *x32*x43    
     2  +coeff( 47)*x15    *x31        
     3  +coeff( 48)    *x25*x31        
     4  +coeff( 49)    *x25    *x41    
     5  +coeff( 50)    *x22*x35        
     6  +coeff( 51)    *x24*x33        
     7  +coeff( 52)*x12*x24    *x41    
     8  +coeff( 53)    *x21*x31    *x51
      y00         =y00         
     9  +coeff( 54)*x15        *x43    
     1  +coeff( 55)        *x32*x41    
     2  +coeff( 56)*x11    *x33        
     3  +coeff( 57)    *x21*x33        
     4  +coeff( 58)*x12*x21*x31        
     5  +coeff( 59)*x11*x22*x31        
     6  +coeff( 60)        *x34*x41    
     7  +coeff( 61)        *x33*x42    
     8  +coeff( 62)        *x31*x44    
      y00         =y00         
     9  +coeff( 63)*x11*x21*x33        
     1  +coeff( 64)*x11*x21*x32*x41    
     2  +coeff( 65)*x12    *x31*x42    
     3  +coeff( 66)*x14    *x31        
     4  +coeff( 67)*x13*x21*x31        
     5  +coeff( 68)    *x24*x31        
     6  +coeff( 69)*x12*x21*x32*x41    
     7  +coeff( 70)*x11*x22*x32*x41    
     8  +coeff( 71)    *x23*x32*x41    
      y00         =y00         
     9  +coeff( 72)*x12*x21*x31*x42    
     1  +coeff( 73)    *x23*x31*x42    
     2  +coeff( 74)*x13        *x43    
     3  +coeff( 75)*x13*x22    *x41    
     4  +coeff( 76)*x11*x24    *x41    
     5  +coeff( 77)*x12*x22*x32*x41    
     6  +coeff( 78)*x11*x23*x32*x41    
     7  +coeff( 79)    *x21    *x41*x51
     8  +coeff( 80)    *x23*x35        
c
      return
      end
