      function x_x2_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.1921339E-01/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.24732618E-01, 0.62007349E-01, 0.16160300E-04, 0.29503512E+00,
     +  0.33592474E-01,-0.18481193E-01, 0.12189335E-01,-0.34367957E-02,
     + -0.90094758E-02,-0.57234867E-02,-0.10617997E-01,-0.44877939E-02,
     + -0.26451452E-02,-0.19768158E-02, 0.81764854E-03,-0.36377122E-02,
     +  0.34601958E-02,-0.35077537E-03, 0.19758109E-02,-0.46431171E-02,
     +  0.60599968E-02, 0.15054771E-02, 0.13421571E-02,-0.20431547E-03,
     +  0.35441015E-03, 0.27333674E-03,-0.85416716E-02,-0.17370562E-03,
     + -0.28273514E-02, 0.53639365E-02,-0.78633591E-03, 0.44080545E-03,
     +  0.36241053E-03,-0.17832190E-02,-0.29944513E-04,-0.18235988E-02,
     + -0.23649146E-02,-0.11383346E-02,-0.16763988E-02,-0.37087894E-02,
     + -0.21462091E-02,-0.14673595E-02,-0.70775009E-03,-0.17202790E-03,
     + -0.26015916E-02,-0.33393959E-03,-0.16173352E-01,-0.10760457E-01,
     + -0.33995609E-02,-0.54657645E-04,-0.10681018E-02, 0.76369097E-03,
     +  0.70984843E-05,-0.68028772E-03,-0.13032008E-02, 0.38799614E-03,
     +  0.47565914E-04,-0.11800259E-02, 0.70468005E-03, 0.27697442E-04,
     +  0.48718706E-04, 0.53181931E-04, 0.19102135E-03,-0.47888450E-03,
     + -0.44033103E-03,-0.22552905E-02,-0.74973553E-02,-0.28336968E-03,
     + -0.91241859E-03, 0.10738690E-03,-0.53190119E-02,-0.84588828E-03,
     + -0.45424258E-04,-0.12064840E-02, 0.38534472E-04,-0.43675984E-04,
     +  0.16844288E-02, 0.38781972E-03, 0.43490398E-03, 0.10380374E-03,
     + -0.11209485E-03,-0.47099101E-03,-0.74967000E-04, 0.24779397E-03,
     + -0.26724183E-02,-0.15160646E-02, 0.12097359E-02, 0.24051622E-02,
     +  0.10826539E-02, 0.59024792E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_x2_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)                *x51
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)                *x52
     7  +coeff(  7)    *x22            
     8  +coeff(  8)*x11                
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff(  9)        *x31*x41    
     1  +coeff( 10)            *x42    
     2  +coeff( 11)    *x21    *x42    
     3  +coeff( 12)    *x21*x31*x41    
     4  +coeff( 13)    *x23*x32        
     5  +coeff( 14)        *x32        
     6  +coeff( 15)    *x23            
     7  +coeff( 16)    *x21*x32        
     8  +coeff( 17)    *x22        *x51
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 18)        *x32    *x51
     1  +coeff( 19)                *x53
     2  +coeff( 20)    *x24            
     3  +coeff( 21)    *x22    *x42    
     4  +coeff( 22)    *x21*x32    *x51
     5  +coeff( 23)        *x31*x41*x51
     6  +coeff( 24)            *x42*x51
     7  +coeff( 25)    *x22*x32        
     8  +coeff( 26)        *x32    *x52
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 27)    *x23*x31*x41    
     1  +coeff( 28)*x11            *x51
     2  +coeff( 29)    *x24        *x51
     3  +coeff( 30)    *x22    *x42*x51
     4  +coeff( 31)    *x21    *x44*x51
     5  +coeff( 32)    *x21*x31*x41*x53
     6  +coeff( 33)    *x21        *x52
     7  +coeff( 34)    *x22*x31*x41    
     8  +coeff( 35)        *x33*x41    
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 36)    *x23        *x51
     1  +coeff( 37)    *x21*x31*x41*x51
     2  +coeff( 38)    *x21    *x42*x51
     3  +coeff( 39)    *x22        *x52
     4  +coeff( 40)    *x23    *x42    
     5  +coeff( 41)    *x23        *x52
     6  +coeff( 42)    *x21*x32    *x52
     7  +coeff( 43)        *x32    *x53
     8  +coeff( 44)*x11*x22            
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 45)    *x24    *x42    
     1  +coeff( 46)                *x54
     2  +coeff( 47)    *x21*x32*x42    
     3  +coeff( 48)    *x21*x31*x43    
     4  +coeff( 49)    *x21    *x44    
     5  +coeff( 50)        *x33*x41*x51
     6  +coeff( 51)        *x32*x42*x51
     7  +coeff( 52)    *x21    *x42*x52
     8  +coeff( 53)            *x43*x52
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 54)    *x22        *x53
     1  +coeff( 55)    *x22*x32*x42    
     2  +coeff( 56)        *x34    *x52
     3  +coeff( 57)    *x23        *x53
     4  +coeff( 58)    *x23*x33*x41    
     5  +coeff( 59)    *x21*x34    *x52
     6  +coeff( 60)            *x41    
     7  +coeff( 61)    *x21*x31        
     8  +coeff( 62)    *x21    *x41    
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 63)        *x31*x43    
     1  +coeff( 64)            *x42*x52
     2  +coeff( 65)    *x21        *x53
     3  +coeff( 66)    *x21*x34        
     4  +coeff( 67)    *x21*x33*x41    
     5  +coeff( 68)        *x34    *x51
     6  +coeff( 69)    *x22*x34        
     7  +coeff( 70)*x11        *x42    
     8  +coeff( 71)    *x22*x31*x43    
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 72)    *x22    *x44    
     1  +coeff( 73)*x11*x21        *x51
     2  +coeff( 74)    *x23    *x42*x51
     3  +coeff( 75)*x11            *x52
     4  +coeff( 76)    *x22*x32    *x52
     5  +coeff( 77)    *x22    *x42*x52
     6  +coeff( 78)    *x21    *x42*x53
     7  +coeff( 79)        *x32    *x54
     8  +coeff( 80)*x11*x23            
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 81)*x11*x21*x32        
     1  +coeff( 82)*x11*x21    *x42    
     2  +coeff( 83)*x11    *x32    *x51
     3  +coeff( 84)    *x24*x32    *x51
     4  +coeff( 85)    *x24    *x42*x51
     5  +coeff( 86)    *x22*x31*x43*x51
     6  +coeff( 87)    *x24*x34        
     7  +coeff( 88)    *x24*x31*x43    
     8  +coeff( 89)    *x24*x31*x41*x52
      x_x2_q3ex   =x_x2_q3ex   
     9  +coeff( 90)*x11*x23    *x42    
c
      return
      end
      function t_x2_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.9425533E-02/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.89198938E-02,-0.36614787E-01, 0.47125241E-05, 0.86585529E-01,
     +  0.33059751E-02,-0.77490099E-02,-0.28406349E-02,-0.28955701E-02,
     + -0.13797076E-02,-0.10528243E-02, 0.79211086E-03,-0.39112664E-03,
     + -0.11167463E-03, 0.14866409E-04,-0.36593695E-03, 0.11894820E-03,
     +  0.19340510E-03,-0.57157152E-03, 0.16038105E-03,-0.71394294E-04,
     +  0.76565752E-03, 0.21310686E-03,-0.13740825E-03,-0.32987536E-03,
     +  0.21934668E-02, 0.38176763E-03,-0.65129623E-03, 0.39923037E-03,
     + -0.23180174E-03, 0.11482339E-03,-0.85406768E-03,-0.11779782E-02,
     + -0.48335473E-03, 0.10738663E-02, 0.13220827E-02,-0.39751557E-03,
     +  0.68234815E-03,-0.20719282E-03, 0.45197882E-03,-0.12491618E-02,
     + -0.96785149E-03,-0.46693732E-03,-0.27621412E-03,-0.31844756E-03,
     +  0.68402092E-05, 0.58933092E-04, 0.25034546E-04, 0.97028336E-04,
     + -0.67247497E-03, 0.29072422E-03,-0.35744699E-03, 0.54175843E-03,
     +  0.42022937E-04, 0.59888983E-03,-0.19073808E-03,-0.20695782E-03,
     + -0.32789318E-03,-0.24868510E-03, 0.12622852E-02, 0.13564499E-02,
     +  0.28831765E-03,-0.35582419E-03,-0.15974186E-03, 0.37316247E-03,
     + -0.16662423E-02, 0.11329796E-02, 0.23574226E-04, 0.13804331E-03,
     +  0.24641748E-04,-0.98293160E-04,-0.87867578E-04, 0.26162959E-04,
     +  0.64681465E-03, 0.58408221E-03,-0.52377389E-04, 0.70808339E-03,
     +  0.28140304E-03, 0.30104447E-04, 0.25368550E-04, 0.28953256E-03,
     + -0.15738982E-03,-0.10823649E-03, 0.74753084E-03, 0.70885353E-03,
     + -0.16330356E-03, 0.36962272E-03, 0.19514210E-03, 0.93129012E-04,
     + -0.17268896E-03, 0.74034539E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      t_x2_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)                *x51
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)                *x52
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)    *x21    *x42    
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff(  9)    *x24            
     1  +coeff( 10)            *x42    
     2  +coeff( 11)                *x53
     3  +coeff( 12)        *x32    *x51
     4  +coeff( 13)    *x21*x34        
     5  +coeff( 14)*x11                
     6  +coeff( 15)        *x32        
     7  +coeff( 16)*x11*x21            
     8  +coeff( 17)    *x23            
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 18)    *x21*x32        
     1  +coeff( 19)    *x21*x31*x41    
     2  +coeff( 20)*x11            *x51
     3  +coeff( 21)        *x31*x41*x51
     4  +coeff( 22)            *x42*x51
     5  +coeff( 23)    *x21        *x52
     6  +coeff( 24)    *x22*x31*x41    
     7  +coeff( 25)    *x22    *x42    
     8  +coeff( 26)    *x21*x32    *x51
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 27)    *x22        *x52
     1  +coeff( 28)        *x32    *x52
     2  +coeff( 29)    *x24        *x52
     3  +coeff( 30)    *x22*x32        
     4  +coeff( 31)    *x23        *x51
     5  +coeff( 32)    *x21*x31*x41*x51
     6  +coeff( 33)    *x21    *x42*x51
     7  +coeff( 34)    *x23    *x42    
     8  +coeff( 35)    *x22    *x42*x51
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 36)    *x22            
     1  +coeff( 37)    *x22        *x51
     2  +coeff( 38)            *x42*x52
     3  +coeff( 39)    *x23*x32        
     4  +coeff( 40)    *x24        *x51
     5  +coeff( 41)    *x23        *x52
     6  +coeff( 42)    *x21*x32    *x52
     7  +coeff( 43)        *x32    *x53
     8  +coeff( 44)    *x24    *x42    
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 45)            *x41    
     1  +coeff( 46)*x11        *x42    
     2  +coeff( 47)        *x32*x42    
     3  +coeff( 48)*x11*x23            
     4  +coeff( 49)        *x32*x42*x51
     5  +coeff( 50)    *x21    *x42*x52
     6  +coeff( 51)    *x22        *x53
     7  +coeff( 52)    *x24*x32        
     8  +coeff( 53)    *x21*x32*x42*x51
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 54)    *x22    *x42*x52
     1  +coeff( 55)    *x23        *x53
     2  +coeff( 56)    *x23*x33*x41    
     3  +coeff( 57)    *x22*x34    *x51
     4  +coeff( 58)    *x24*x31*x41*x51
     5  +coeff( 59)    *x21*x33*x41*x52
     6  +coeff( 60)    *x21*x31*x43*x52
     7  +coeff( 61)    *x21    *x44*x52
     8  +coeff( 62)        *x34    *x53
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 63)    *x23    *x44*x51
     1  +coeff( 64)        *x34    *x54
     2  +coeff( 65)    *x22*x32*x42*x53
     3  +coeff( 66)    *x24*x33*x41*x52
     4  +coeff( 67)*x11    *x32        
     5  +coeff( 68)        *x33*x41    
     6  +coeff( 69)*x11            *x52
     7  +coeff( 70)    *x21        *x53
     8  +coeff( 71)                *x54
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 72)*x11*x21*x31*x41    
     1  +coeff( 73)    *x23*x31*x41    
     2  +coeff( 74)    *x21*x33*x41    
     3  +coeff( 75)*x11*x21    *x42    
     4  +coeff( 76)    *x21*x31*x43    
     5  +coeff( 77)    *x21    *x44    
     6  +coeff( 78)*x11*x22        *x51
     7  +coeff( 79)*x11*x21        *x52
     8  +coeff( 80)        *x31*x41*x53
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 81)            *x42*x53
     1  +coeff( 82)    *x21        *x54
     2  +coeff( 83)    *x24*x31*x41    
     3  +coeff( 84)    *x22*x32*x42    
     4  +coeff( 85)    *x22    *x44    
     5  +coeff( 86)    *x23*x31*x41*x51
     6  +coeff( 87)    *x22*x31*x41*x52
     7  +coeff( 88)        *x31*x41*x54
     8  +coeff( 89)            *x42*x54
      t_x2_q3ex   =t_x2_q3ex   
     9  +coeff( 90)    *x21*x34*x42    
c
      return
      end
      function y_x2_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 92)
      data ncoeff/ 91/
      data avdat/ -0.6881399E-03/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.69755345E-03,-0.12262969E+00, 0.10005265E+00, 0.15231663E-03,
     + -0.38510054E-01,-0.77267632E-01, 0.37237044E-01, 0.69166698E-01,
     + -0.15352550E-01,-0.18436385E-01,-0.42877123E-01,-0.13441253E-01,
     + -0.12031127E-01,-0.76802764E-02,-0.41291099E-02,-0.52013995E-04,
     +  0.13094076E-02,-0.29289417E-02, 0.74968194E-02,-0.48696348E-02,
     + -0.60510258E-02,-0.17007880E-04,-0.15608560E-03, 0.12954762E-02,
     + -0.27817395E-02, 0.74279500E-03, 0.24013345E-02, 0.33280306E-03,
     +  0.17644183E-02,-0.47628203E-03, 0.10131260E-03,-0.49245131E-03,
     +  0.35674339E-02, 0.60201553E-02, 0.20091219E-02,-0.81390369E-03,
     + -0.70623955E-03,-0.11306390E-01,-0.15062844E-01,-0.66433735E-02,
     + -0.15229724E-02,-0.25383339E-02, 0.34136360E-02, 0.20337850E-02,
     + -0.18379791E-03, 0.15806088E-02,-0.20218511E-03, 0.37033830E-03,
     +  0.30465017E-03, 0.63443667E-03, 0.77171106E-04, 0.14313055E-02,
     +  0.89898903E-03,-0.27827015E-02, 0.12944260E-02,-0.31590488E-03,
     + -0.89730660E-03, 0.42624064E-03, 0.59887342E-03, 0.11657162E-03,
     +  0.50776358E-04, 0.39849212E-04, 0.62371517E-03,-0.64672444E-04,
     + -0.44131710E-03, 0.65838516E-03, 0.85716048E-04, 0.36212736E-04,
     + -0.47540246E-03, 0.17016068E-03, 0.36306770E-02, 0.24148193E-02,
     +  0.15275946E-02,-0.25532756E-02, 0.27686261E-03, 0.69563661E-03,
     + -0.61431463E-04,-0.11130027E-03, 0.10253820E-02, 0.37599390E-03,
     + -0.25944770E-03, 0.37296311E-03,-0.16622989E-02,-0.18027540E-02,
     + -0.25449743E-03,-0.11713900E-01,-0.82305716E-02,-0.23869576E-02,
     +  0.86616878E-04,-0.44417506E-03, 0.48629041E-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
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*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
      x45 = x44*x4
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      y_x2_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)    *x21            
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff(  9)        *x31*x42    
     1  +coeff( 10)    *x22*x31        
     2  +coeff( 11)    *x22    *x41    
     3  +coeff( 12)    *x21    *x41*x51
     4  +coeff( 13)        *x32*x41    
     5  +coeff( 14)            *x43    
     6  +coeff( 15)    *x21*x31    *x51
     7  +coeff( 16)        *x35        
     8  +coeff( 17)*x11        *x41    
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 18)            *x41*x52
     1  +coeff( 19)    *x23    *x41    
     2  +coeff( 20)    *x22*x31    *x51
     3  +coeff( 21)    *x22    *x41*x51
     4  +coeff( 22)*x11    *x33        
     5  +coeff( 23)                *x51
     6  +coeff( 24)*x11*x21    *x41    
     7  +coeff( 25)        *x33        
     8  +coeff( 26)*x11    *x31        
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 27)    *x21    *x43    
     1  +coeff( 28)*x11*x21*x31        
     2  +coeff( 29)    *x23*x31        
     3  +coeff( 30)    *x21    *x43*x52
     4  +coeff( 31)    *x22            
     5  +coeff( 32)        *x31    *x52
     6  +coeff( 33)    *x21*x32*x41    
     7  +coeff( 34)    *x21*x31*x42    
     8  +coeff( 35)        *x31*x42*x51
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 36)        *x31    *x53
     1  +coeff( 37)            *x41*x53
     2  +coeff( 38)    *x22*x32*x41    
     3  +coeff( 39)    *x22*x31*x42    
     4  +coeff( 40)    *x22    *x43    
     5  +coeff( 41)    *x21*x32*x41*x51
     6  +coeff( 42)    *x24*x31        
     7  +coeff( 43)    *x23    *x41*x51
     8  +coeff( 44)    *x22    *x41*x52
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 45)    *x21*x35        
     1  +coeff( 46)    *x23*x33        
     2  +coeff( 47)    *x21*x33    *x52
     3  +coeff( 48)        *x32*x41*x53
     4  +coeff( 49)    *x23    *x41*x52
     5  +coeff( 50)    *x22*x35        
     6  +coeff( 51)            *x42    
     7  +coeff( 52)    *x21*x33        
     8  +coeff( 53)    *x21    *x41*x52
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 54)    *x22*x33        
     1  +coeff( 55)    *x21*x31*x42*x51
     2  +coeff( 56)*x11*x22    *x41    
     3  +coeff( 57)    *x24    *x41    
     4  +coeff( 58)        *x31    *x54
     5  +coeff( 59)    *x23    *x45    
     6  +coeff( 60)        *x31*x41    
     7  +coeff( 61)    *x22        *x51
     8  +coeff( 62)        *x34        
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 63)        *x32*x41*x51
     1  +coeff( 64)*x11    *x31    *x51
     2  +coeff( 65)    *x21*x31    *x52
     3  +coeff( 66)        *x31*x42*x52
     4  +coeff( 67)*x11*x21*x31    *x51
     5  +coeff( 68)*x12*x21            
     6  +coeff( 69)    *x21*x34*x41    
     7  +coeff( 70)        *x35    *x51
     8  +coeff( 71)    *x23*x32*x41    
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 72)    *x23*x31*x42    
     1  +coeff( 73)    *x23    *x43    
     2  +coeff( 74)    *x22*x31*x42*x51
     3  +coeff( 75)    *x22    *x43*x51
     4  +coeff( 76)    *x21*x31*x42*x52
     5  +coeff( 77)*x11*x22    *x42    
     6  +coeff( 78)*x11*x23    *x41    
     7  +coeff( 79)    *x22    *x41*x53
     8  +coeff( 80)    *x21*x33    *x53
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 81)            *x43*x54
     1  +coeff( 82)    *x24*x31    *x52
     2  +coeff( 83)    *x24    *x43*x51
     3  +coeff( 84)    *x22*x32*x41*x53
     4  +coeff( 85)*x12    *x31*x41*x52
     5  +coeff( 86)    *x24*x32*x43    
     6  +coeff( 87)    *x24*x31*x44    
     7  +coeff( 88)    *x24    *x45    
     8  +coeff( 89)*x12*x23        *x51
      y_x2_q3ex   =y_x2_q3ex   
     9  +coeff( 90)*x12    *x31*x42*x52
     1  +coeff( 91)    *x23*x34*x43    
c
      return
      end
      function p_x2_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/  0.2195690E-03/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.19273526E-03,-0.31512227E-04, 0.16610067E-01,-0.33561151E-01,
     +  0.59880586E-02, 0.16069550E-01,-0.11128638E-01,-0.10222793E-01,
     +  0.61070374E-02, 0.17313984E-02, 0.82812150E-03, 0.40299143E-02,
     +  0.19031615E-02, 0.34315011E-02, 0.13708723E-02, 0.27628909E-03,
     + -0.35174282E-04, 0.65306731E-03,-0.28817100E-03, 0.22936820E-02,
     + -0.77017868E-03, 0.31028838E-04,-0.10122289E-03, 0.45201700E-03,
     + -0.87302370E-03,-0.63032529E-03, 0.20532850E-02, 0.10918867E-02,
     +  0.20945759E-03,-0.54895191E-03,-0.10976530E-03,-0.15308170E-03,
     + -0.12677560E-02,-0.43000988E-03,-0.25786058E-03, 0.39225123E-02,
     +  0.49953605E-03,-0.57523005E-03, 0.36066329E-04, 0.40381434E-03,
     + -0.10142152E-03,-0.95010514E-03, 0.15235379E-04,-0.35432901E-03,
     + -0.72125852E-03,-0.65311708E-03,-0.60012279E-03,-0.14489674E-03,
     +  0.92831225E-03, 0.69045881E-03, 0.24080505E-02, 0.35941196E-03,
     + -0.49427594E-03, 0.33439617E-03, 0.44165028E-03, 0.54954190E-03,
     +  0.84951846E-03,-0.64648280E-03,-0.44909288E-03, 0.50371062E-04,
     + -0.31917842E-03,-0.15855834E-04,-0.33826631E-03, 0.48323677E-03,
     +  0.88756278E-04, 0.21131093E-05,-0.46962989E-04,-0.29617202E-03,
     +  0.17416292E-03, 0.33825322E-03, 0.43846449E-03, 0.55187318E-03,
     + -0.55457222E-04,-0.83503601E-05,-0.45743600E-05,-0.55270171E-04,
     +  0.41653329E-04,-0.42964937E-04,-0.11933495E-02,-0.58944442E-03,
     + -0.19076827E-03, 0.44397431E-04, 0.80645412E-04,-0.64792641E-03,
     + -0.35455750E-03, 0.58866353E-04, 0.19301870E-03, 0.50943607E-04,
     + -0.15860773E-03, 0.22227698E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      p_x2_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff(  9)    *x22    *x41    
     1  +coeff( 10)        *x31    *x52
     2  +coeff( 11)    *x24*x31        
     3  +coeff( 12)    *x22*x31        
     4  +coeff( 13)            *x41*x52
     5  +coeff( 14)        *x31*x42    
     6  +coeff( 15)            *x43    
     7  +coeff( 16)    *x21*x31    *x51
     8  +coeff( 17)        *x34*x41    
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 18)        *x33        
     1  +coeff( 19)*x11        *x41    
     2  +coeff( 20)        *x32*x41    
     3  +coeff( 21)        *x31    *x53
     4  +coeff( 22)                *x51
     5  +coeff( 23)*x11    *x31        
     6  +coeff( 24)    *x21    *x41*x51
     7  +coeff( 25)    *x23    *x41    
     8  +coeff( 26)            *x41*x53
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 27)    *x22    *x43    
     1  +coeff( 28)        *x31    *x54
     2  +coeff( 29)    *x22    *x43*x52
     3  +coeff( 30)    *x23*x31        
     4  +coeff( 31)    *x21*x33        
     5  +coeff( 32)*x11*x21    *x41    
     6  +coeff( 33)    *x21*x31*x42    
     7  +coeff( 34)    *x21    *x43    
     8  +coeff( 35)        *x33    *x51
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 36)    *x22*x31*x42    
     1  +coeff( 37)    *x21*x31    *x53
     2  +coeff( 38)    *x23*x32*x41    
     3  +coeff( 39)    *x21*x34*x41    
     4  +coeff( 40)    *x24*x32*x41    
     5  +coeff( 41)*x11*x21*x31        
     6  +coeff( 42)    *x21*x32*x41    
     7  +coeff( 43)*x11        *x41*x51
     8  +coeff( 44)    *x22    *x41*x51
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 45)        *x32*x41*x51
     1  +coeff( 46)        *x31*x42*x51
     2  +coeff( 47)            *x43*x51
     3  +coeff( 48)    *x21*x31    *x52
     4  +coeff( 49)    *x22*x33        
     5  +coeff( 50)    *x24    *x41    
     6  +coeff( 51)    *x22*x32*x41    
     7  +coeff( 52)    *x23*x31    *x51
     8  +coeff( 53)    *x21*x32*x41*x51
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 54)    *x21    *x43*x51
     1  +coeff( 55)    *x22*x31    *x52
     2  +coeff( 56)        *x33    *x52
     3  +coeff( 57)        *x31*x42*x52
     4  +coeff( 58)    *x21    *x41*x53
     5  +coeff( 59)            *x41*x54
     6  +coeff( 60)    *x24*x31    *x51
     7  +coeff( 61)        *x33    *x53
     8  +coeff( 62)    *x22            
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 63)    *x22*x31    *x51
     1  +coeff( 64)    *x21*x31*x42*x51
     2  +coeff( 65)    *x22    *x41*x52
     3  +coeff( 66)    *x23*x33        
     4  +coeff( 67)*x11*x22*x31    *x51
     5  +coeff( 68)    *x23    *x41*x52
     6  +coeff( 69)    *x21*x32*x41*x52
     7  +coeff( 70)    *x22*x31    *x53
     8  +coeff( 71)    *x21*x31    *x54
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 72)    *x22*x31*x44    
     1  +coeff( 73)    *x22*x34*x41*x51
     2  +coeff( 74)            *x42    
     3  +coeff( 75)                *x52
     4  +coeff( 76)    *x22*x31*x41    
     5  +coeff( 77)    *x21*x33    *x51
     6  +coeff( 78)    *x24    *x42    
     7  +coeff( 79)    *x23*x31*x42    
     8  +coeff( 80)    *x23    *x43    
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 81)    *x21*x31*x44    
     1  +coeff( 82)    *x24    *x41*x51
     2  +coeff( 83)*x11    *x32*x41*x51
     3  +coeff( 84)    *x22*x32*x41*x51
     4  +coeff( 85)    *x22*x31*x42*x51
     5  +coeff( 86)    *x22    *x43*x51
     6  +coeff( 87)    *x21*x33    *x52
     7  +coeff( 88)*x11        *x41*x53
     8  +coeff( 89)    *x22    *x41*x53
      p_x2_q3ex   =p_x2_q3ex   
     9  +coeff( 90)        *x31*x42*x53
c
      return
      end
      function l_x2_q3ex   (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 51)
      data ncoeff/ 50/
      data avdat/ -0.2328331E-01/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.12743056E-01,-0.32734451E+00,-0.32093614E-01, 0.61825612E-02,
     + -0.37352514E-01,-0.13765091E-01,-0.52409181E-02,-0.27748151E-03,
     + -0.69273114E-02,-0.67291399E-02, 0.12216757E-02, 0.10456796E-01,
     +  0.75787031E-02,-0.28027298E-02, 0.64455564E-02,-0.48650340E-02,
     + -0.12646365E-02, 0.44567636E-02,-0.59576560E-03, 0.16518802E-02,
     +  0.28077343E-02, 0.92854310E-03, 0.55605004E-03, 0.36703228E-03,
     +  0.36586256E-03, 0.40806727E-02, 0.27629323E-02,-0.10134117E-02,
     + -0.88810449E-03,-0.16961226E-02, 0.11722791E-01, 0.98427152E-02,
     +  0.76053751E-03, 0.10035564E-02, 0.18054426E-01, 0.17094130E-01,
     +  0.48521906E-02, 0.76333311E-03,-0.81579667E-03, 0.17810795E-01,
     +  0.18173390E-02, 0.10363160E-01, 0.37488213E-02,-0.24426891E-02,
     +  0.15311138E-01, 0.37680067E-04, 0.33392458E-04, 0.74923685E-03,
     + -0.90082103E-04,-0.44015489E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      l_x2_q3ex   =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)*x11                
     5  +coeff(  5)    *x22            
     6  +coeff(  6)    *x21        *x51
     7  +coeff(  7)                *x52
     8  +coeff(  8)        *x34        
      l_x2_q3ex   =l_x2_q3ex   
     9  +coeff(  9)        *x31*x41    
     1  +coeff( 10)            *x42    
     2  +coeff( 11)*x11*x21            
     3  +coeff( 12)    *x21*x31*x41    
     4  +coeff( 13)    *x21    *x42    
     5  +coeff( 14)    *x22        *x51
     6  +coeff( 15)    *x23*x32        
     7  +coeff( 16)        *x32        
     8  +coeff( 17)    *x23            
      l_x2_q3ex   =l_x2_q3ex   
     9  +coeff( 18)    *x21*x32        
     1  +coeff( 19)*x11            *x51
     2  +coeff( 20)        *x32    *x51
     3  +coeff( 21)        *x31*x41*x51
     4  +coeff( 22)                *x53
     5  +coeff( 23)            *x42*x51
     6  +coeff( 24)    *x21        *x52
     7  +coeff( 25)*x11*x22            
     8  +coeff( 26)    *x22*x31*x41    
      l_x2_q3ex   =l_x2_q3ex   
     9  +coeff( 27)    *x22    *x42    
     1  +coeff( 28)    *x23        *x51
     2  +coeff( 29)    *x21*x32    *x51
     3  +coeff( 30)    *x21*x31*x41*x51
     4  +coeff( 31)    *x23*x31*x41    
     5  +coeff( 32)    *x23    *x42    
     6  +coeff( 33)    *x24*x32        
     7  +coeff( 34)    *x24            
     8  +coeff( 35)    *x21*x32*x42    
      l_x2_q3ex   =l_x2_q3ex   
     9  +coeff( 36)    *x21*x31*x43    
     1  +coeff( 37)    *x21    *x44    
     2  +coeff( 38)    *x22*x31*x43    
     3  +coeff( 39)        *x31*x41*x54
     4  +coeff( 40)    *x23*x33*x41    
     5  +coeff( 41)    *x21*x34    *x52
     6  +coeff( 42)    *x21*x33*x41*x52
     7  +coeff( 43)    *x21*x32*x42*x52
     8  +coeff( 44)    *x21*x31*x41*x54
      l_x2_q3ex   =l_x2_q3ex   
     9  +coeff( 45)    *x23*x34*x42    
     1  +coeff( 46)        *x31        
     2  +coeff( 47)            *x41    
     3  +coeff( 48)    *x22*x32        
     4  +coeff( 49)*x11*x21        *x51
     5  +coeff( 50)            *x42*x52
c
      return
      end
      function x_x2_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.4634684E-01/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.50390247E-01,-0.43253645E-01, 0.30287520E-04, 0.54434544E+00,
     +  0.43222692E-01,-0.40699638E-01, 0.11065912E-01,-0.16875524E-01,
     + -0.84506152E-02,-0.19969482E-01,-0.33707195E-02,-0.63309129E-02,
     + -0.51632733E-02, 0.51939590E-02,-0.12209265E-02, 0.42039468E-02,
     + -0.85830726E-02,-0.70822731E-04, 0.11517089E-01, 0.11484432E-02,
     +  0.37151466E-02, 0.37791880E-04, 0.25112766E-02,-0.24025473E-02,
     + -0.29012773E-02,-0.10821061E-03,-0.50338828E-02,-0.27310706E-02,
     +  0.79615560E-03,-0.88823779E-03, 0.20918139E-03,-0.60451701E-02,
     + -0.42009878E-03,-0.60467473E-02, 0.98623596E-02,-0.20375010E-02,
     + -0.32222306E-05,-0.78066376E-04,-0.55246940E-02,-0.45002280E-02,
     + -0.34954303E-02,-0.47329585E-02,-0.36799912E-02, 0.44355918E-04,
     + -0.14622423E-02,-0.22239746E-02, 0.21107064E-02,-0.17111229E-02,
     + -0.17814325E-03, 0.21939816E-03,-0.16717099E-02, 0.11735458E-03,
     +  0.34171401E-02, 0.60343364E-03, 0.12481645E-02,-0.11256186E-02,
     + -0.40048868E-02,-0.37086741E-05, 0.85230498E-02,-0.10031913E-03,
     + -0.94633870E-03,-0.38791072E-03,-0.71899488E-03,-0.14166668E-02,
     + -0.43275352E-02,-0.12330211E-01,-0.56000114E-02,-0.13857855E-02,
     + -0.87493157E-03, 0.13029572E-02,-0.40894040E-03, 0.24453297E-02,
     +  0.43919832E-02, 0.34968758E-02, 0.60233427E-02, 0.22228502E-02,
     + -0.66558551E-03, 0.90180431E-03, 0.17110087E-02, 0.14437335E-02,
     + -0.12313259E-02,-0.23139888E-03, 0.15806497E-03,-0.37063113E-02,
     +  0.12836579E-01,-0.10568101E-02, 0.65387132E-04, 0.72307077E-04,
     +  0.26469481E-04,-0.37140836E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      x_x2_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)                *x51
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)                *x52
     7  +coeff(  7)    *x22            
     8  +coeff(  8)        *x31*x41    
      x_x2_fp     =x_x2_fp     
     9  +coeff(  9)            *x42    
     1  +coeff( 10)    *x21    *x42    
     2  +coeff( 11)*x11                
     3  +coeff( 12)    *x21*x32        
     4  +coeff( 13)    *x21*x31*x41    
     5  +coeff( 14)    *x22        *x51
     6  +coeff( 15)        *x32    *x51
     7  +coeff( 16)                *x53
     8  +coeff( 17)    *x24            
      x_x2_fp     =x_x2_fp     
     9  +coeff( 18)        *x34        
     1  +coeff( 19)    *x22    *x42    
     2  +coeff( 20)    *x23            
     3  +coeff( 21)        *x31*x41*x51
     4  +coeff( 22)            *x42*x51
     5  +coeff( 23)    *x21*x32    *x51
     6  +coeff( 24)    *x21*x32    *x52
     7  +coeff( 25)        *x32        
     8  +coeff( 26)    *x22*x32        
      x_x2_fp     =x_x2_fp     
     9  +coeff( 27)    *x21*x31*x41*x51
     1  +coeff( 28)    *x21    *x42*x51
     2  +coeff( 29)        *x32    *x52
     3  +coeff( 30)    *x21        *x53
     4  +coeff( 31)*x11*x21            
     5  +coeff( 32)    *x23*x31*x41    
     6  +coeff( 33)*x11            *x51
     7  +coeff( 34)    *x24        *x51
     8  +coeff( 35)    *x22    *x42*x51
      x_x2_fp     =x_x2_fp     
     9  +coeff( 36)        *x32    *x53
     1  +coeff( 37)    *x23*x32    *x51
     2  +coeff( 38)    *x21        *x52
     3  +coeff( 39)    *x22*x31*x41    
     4  +coeff( 40)    *x23        *x51
     5  +coeff( 41)    *x22        *x52
     6  +coeff( 42)    *x23        *x52
     7  +coeff( 43)    *x24    *x42    
     8  +coeff( 44)            *x41    
      x_x2_fp     =x_x2_fp     
     9  +coeff( 45)            *x42*x52
     1  +coeff( 46)        *x32*x42*x51
     2  +coeff( 47)    *x21    *x42*x52
     3  +coeff( 48)    *x22        *x53
     4  +coeff( 49)*x11*x22            
     5  +coeff( 50)*x11        *x42    
     6  +coeff( 51)    *x23    *x42*x51
     7  +coeff( 52)*x11            *x52
     8  +coeff( 53)    *x22    *x42*x52
      x_x2_fp     =x_x2_fp     
     9  +coeff( 54)*x11*x23            
     1  +coeff( 55)    *x23*x32*x42    
     2  +coeff( 56)    *x21*x34*x42    
     3  +coeff( 57)    *x21*x33*x43    
     4  +coeff( 58)    *x21*x32*x42*x52
     5  +coeff( 59)        *x33*x43*x53
     6  +coeff( 60)        *x33*x41    
     7  +coeff( 61)        *x32*x42    
     8  +coeff( 62)            *x44    
      x_x2_fp     =x_x2_fp     
     9  +coeff( 63)                *x54
     1  +coeff( 64)    *x21*x34        
     2  +coeff( 65)    *x21*x33*x41    
     3  +coeff( 66)    *x21*x32*x42    
     4  +coeff( 67)    *x21*x31*x43    
     5  +coeff( 68)    *x21    *x44    
     6  +coeff( 69)        *x34    *x51
     7  +coeff( 70)    *x21*x31*x41*x52
     8  +coeff( 71)    *x21        *x54
      x_x2_fp     =x_x2_fp     
     9  +coeff( 72)    *x24*x32        
     1  +coeff( 73)    *x24*x31*x41    
     2  +coeff( 74)    *x22*x33*x41    
     3  +coeff( 75)    *x22*x32*x42    
     4  +coeff( 76)    *x21*x32*x42*x51
     5  +coeff( 77)    *x24        *x52
     6  +coeff( 78)        *x34    *x52
     7  +coeff( 79)    *x22*x31*x41*x52
     8  +coeff( 80)        *x32    *x54
      x_x2_fp     =x_x2_fp     
     9  +coeff( 81)    *x23*x34        
     1  +coeff( 82)*x11*x21    *x42    
     2  +coeff( 83)*x11*x22        *x51
     3  +coeff( 84)    *x24    *x42*x51
     4  +coeff( 85)    *x21*x33*x43*x52
     5  +coeff( 86)    *x22*x34    *x53
     6  +coeff( 87)    *x21*x31        
     7  +coeff( 88)    *x21    *x41    
     8  +coeff( 89)        *x31    *x51
      x_x2_fp     =x_x2_fp     
     9  +coeff( 90)        *x31*x43    
c
      return
      end
      function t_x2_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/ -0.9425379E-02/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.89152008E-02,-0.36627874E-01, 0.41775425E-05, 0.86573504E-01,
     +  0.32525421E-02,-0.77154851E-02,-0.27859453E-02,-0.27777059E-02,
     + -0.14210464E-02,-0.10372612E-02, 0.81692124E-03,-0.21152597E-03,
     + -0.93520735E-04, 0.16309677E-04,-0.35302146E-03, 0.12565397E-03,
     +  0.22632370E-03,-0.53907163E-03, 0.21175301E-03,-0.74789736E-04,
     +  0.81565313E-03, 0.21335397E-03,-0.10171642E-03,-0.74878940E-03,
     +  0.21563980E-02, 0.50602620E-03,-0.91616384E-03,-0.74366812E-03,
     +  0.19912602E-03,-0.13239810E-03,-0.72862825E-03,-0.36878220E-03,
     + -0.78921654E-03, 0.15512814E-02, 0.71356457E-03, 0.26031828E-03,
     +  0.29154524E-03,-0.35034717E-03,-0.12267587E-02,-0.10317194E-02,
     + -0.48522788E-03,-0.43632471E-03,-0.54159237E-03, 0.58638689E-05,
     +  0.67029160E-03,-0.13630714E-03, 0.97622629E-04, 0.31584711E-03,
     + -0.62708117E-04, 0.66629995E-03, 0.11545674E-03,-0.26343064E-03,
     + -0.85022685E-03, 0.30934007E-03, 0.11348944E-02, 0.66262578E-05,
     + -0.84789775E-04,-0.50859677E-03,-0.46379626E-03, 0.96255944E-05,
     +  0.71742793E-03, 0.20951761E-02,-0.22709179E-03, 0.47226631E-03,
     + -0.88956462E-04, 0.53808722E-03, 0.23533798E-04, 0.49284736E-04,
     +  0.23674407E-04,-0.12484554E-03, 0.84311113E-03, 0.52077149E-03,
     +  0.37255097E-04,-0.12664151E-03,-0.22629254E-03, 0.42398460E-03,
     +  0.26559207E-03,-0.18635472E-03,-0.19575357E-03, 0.75662805E-03,
     + -0.45160751E-03, 0.13659429E-02,-0.82002283E-04, 0.21566627E-03,
     +  0.26348332E-03, 0.23273472E-03, 0.31035958E-03, 0.14476529E-02,
     +  0.68099314E-03,-0.48359216E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      t_x2_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)                *x51
     5  +coeff(  5)    *x21        *x51
     6  +coeff(  6)                *x52
     7  +coeff(  7)        *x31*x41    
     8  +coeff(  8)    *x21    *x42    
      t_x2_fp     =t_x2_fp     
     9  +coeff(  9)    *x24            
     1  +coeff( 10)            *x42    
     2  +coeff( 11)                *x53
     3  +coeff( 12)        *x32    *x51
     4  +coeff( 13)    *x21*x34        
     5  +coeff( 14)*x11                
     6  +coeff( 15)        *x32        
     7  +coeff( 16)*x11*x21            
     8  +coeff( 17)    *x23            
      t_x2_fp     =t_x2_fp     
     9  +coeff( 18)    *x21*x32        
     1  +coeff( 19)    *x21*x31*x41    
     2  +coeff( 20)*x11            *x51
     3  +coeff( 21)        *x31*x41*x51
     4  +coeff( 22)            *x42*x51
     5  +coeff( 23)    *x21        *x52
     6  +coeff( 24)    *x22*x31*x41    
     7  +coeff( 25)    *x22    *x42    
     8  +coeff( 26)    *x21*x32    *x51
      t_x2_fp     =t_x2_fp     
     9  +coeff( 27)    *x21*x31*x41*x51
     1  +coeff( 28)    *x22        *x52
     2  +coeff( 29)        *x32    *x52
     3  +coeff( 30)    *x23*x32    *x51
     4  +coeff( 31)    *x23    *x42*x51
     5  +coeff( 32)    *x22            
     6  +coeff( 33)    *x23        *x51
     7  +coeff( 34)    *x22    *x42*x51
     8  +coeff( 35)    *x24*x32        
      t_x2_fp     =t_x2_fp     
     9  +coeff( 36)    *x24        *x53
     1  +coeff( 37)    *x23        *x54
     2  +coeff( 38)            *x42*x52
     3  +coeff( 39)    *x24        *x51
     4  +coeff( 40)    *x23        *x52
     5  +coeff( 41)    *x21*x32    *x52
     6  +coeff( 42)    *x22        *x53
     7  +coeff( 43)        *x32    *x53
     8  +coeff( 44)            *x41    
      t_x2_fp     =t_x2_fp     
     9  +coeff( 45)    *x22        *x51
     1  +coeff( 46)    *x21        *x53
     2  +coeff( 47)*x11*x23            
     3  +coeff( 48)    *x23*x32        
     4  +coeff( 49)*x11*x21    *x42    
     5  +coeff( 50)    *x23    *x42    
     6  +coeff( 51)    *x21    *x44    
     7  +coeff( 52)    *x22*x31*x41*x51
     8  +coeff( 53)        *x32*x42*x51
      t_x2_fp     =t_x2_fp     
     9  +coeff( 54)    *x21    *x42*x52
     1  +coeff( 55)    *x24*x31*x41    
     2  +coeff( 56)*x11    *x32*x42    
     3  +coeff( 57)    *x22    *x44    
     4  +coeff( 58)    *x21*x31*x43*x51
     5  +coeff( 59)    *x21    *x44*x51
     6  +coeff( 60)*x11        *x42*x52
     7  +coeff( 61)    *x22    *x42*x52
     8  +coeff( 62)    *x23*x31*x43    
      t_x2_fp     =t_x2_fp     
     9  +coeff( 63)    *x22*x34    *x51
     1  +coeff( 64)    *x23*x31*x41*x52
     2  +coeff( 65)    *x22*x33*x43    
     3  +coeff( 66)    *x22*x33*x41*x52
     4  +coeff( 67)*x11    *x32        
     5  +coeff( 68)*x11        *x42    
     6  +coeff( 69)*x11            *x52
     7  +coeff( 70)                *x54
     8  +coeff( 71)    *x21*x33*x41    
      t_x2_fp     =t_x2_fp     
     9  +coeff( 72)    *x21*x31*x43    
     1  +coeff( 73)*x11*x22        *x51
     2  +coeff( 74)    *x22*x32    *x51
     3  +coeff( 75)        *x34    *x51
     4  +coeff( 76)    *x21*x31*x41*x52
     5  +coeff( 77)        *x31*x41*x53
     6  +coeff( 78)            *x42*x53
     7  +coeff( 79)    *x21        *x54
     8  +coeff( 80)    *x22*x33*x41    
      t_x2_fp     =t_x2_fp     
     9  +coeff( 81)    *x24    *x42    
     1  +coeff( 82)    *x22*x32*x42    
     2  +coeff( 83)    *x21*x34    *x51
     3  +coeff( 84)        *x34    *x52
     4  +coeff( 85)    *x22*x31*x41*x52
     5  +coeff( 86)        *x31*x43*x52
     6  +coeff( 87)        *x32    *x54
     7  +coeff( 88)    *x23*x32*x42    
     8  +coeff( 89)    *x23    *x44    
      t_x2_fp     =t_x2_fp     
     9  +coeff( 90)    *x24    *x42*x51
c
      return
      end
      function y_x2_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 96)
      data ncoeff/ 95/
      data avdat/ -0.5585809E-04/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.14894304E-03,-0.74744269E-01, 0.33224355E-02, 0.70478265E-04,
     + -0.21307042E-01,-0.30910566E-01, 0.49445168E-02, 0.39757922E-01,
     + -0.53786873E-02,-0.24904139E-01,-0.12211530E-01, 0.43348470E-02,
     + -0.70369090E-02, 0.71503856E-03, 0.23947221E-04,-0.51600658E-02,
     + -0.35125085E-02,-0.71544489E-02, 0.31132845E-02, 0.48342613E-02,
     + -0.56922217E-02, 0.44504611E-03, 0.48295563E-03,-0.34815939E-02,
     +  0.77497843E-03,-0.27466116E-02,-0.23954946E-02,-0.11233514E-02,
     +  0.11679013E-02, 0.49218669E-03,-0.33028279E-02, 0.34967943E-02,
     +  0.83461654E-03, 0.37253043E-02,-0.27716026E-03, 0.24070439E-04,
     + -0.66377615E-04, 0.11754562E-02,-0.18246602E-02, 0.22603462E-03,
     +  0.58901252E-03,-0.45076674E-02,-0.39076605E-02,-0.18457853E-02,
     +  0.16390544E-02, 0.29590505E-02, 0.98695117E-03, 0.18533802E-02,
     + -0.17927826E-02,-0.16925473E-02, 0.10364651E-02, 0.99757139E-03,
     + -0.15856174E-02,-0.46387044E-03,-0.91327951E-04, 0.64659590E-03,
     + -0.20121832E-02,-0.12847641E-02, 0.25359555E-02, 0.76801877E-03,
     + -0.34216102E-03, 0.91426540E-03, 0.16381727E-02, 0.46059542E-03,
     +  0.28788613E-02, 0.27774673E-03, 0.10234929E-02, 0.40244372E-02,
     +  0.46381526E-03, 0.46390731E-04, 0.26937269E-04, 0.56123241E-04,
     +  0.15467550E-02, 0.58822196E-04, 0.64269821E-04,-0.78013261E-04,
     +  0.13261668E-03,-0.35599314E-03,-0.74179831E-03, 0.21816617E-03,
     + -0.19984231E-03, 0.14166208E-02,-0.29940403E-02,-0.28686266E-03,
     +  0.49890205E-04,-0.27700773E-03, 0.92070410E-03, 0.57268696E-03,
     +  0.10781684E-02,-0.99665974E-03,-0.28625017E-03, 0.13134154E-02,
     +  0.74160902E-03, 0.13944325E-02,-0.24010667E-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
      x21 = x2
      x22 = x21*x2
      x23 = x22*x2
      x24 = x23*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
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      y_x2_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)        *x31        
     3  +coeff(  3)            *x41    
     4  +coeff(  4)    *x21            
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      y_x2_fp     =y_x2_fp     
     9  +coeff(  9)        *x31*x42    
     1  +coeff( 10)    *x22    *x41    
     2  +coeff( 11)    *x21    *x41*x51
     3  +coeff( 12)        *x31    *x52
     4  +coeff( 13)    *x22    *x41*x51
     5  +coeff( 14)    *x22*x33        
     6  +coeff( 15)    *x24*x31        
     7  +coeff( 16)        *x32*x41    
     8  +coeff( 17)            *x43    
      y_x2_fp     =y_x2_fp     
     9  +coeff( 18)    *x22*x31        
     1  +coeff( 19)            *x41*x52
     2  +coeff( 20)    *x23    *x41    
     3  +coeff( 21)    *x22*x31    *x51
     4  +coeff( 22)*x11    *x31        
     5  +coeff( 23)*x11        *x41    
     6  +coeff( 24)    *x21*x31    *x51
     7  +coeff( 25)*x11*x21    *x41    
     8  +coeff( 26)        *x31    *x53
      y_x2_fp     =y_x2_fp     
     9  +coeff( 27)            *x41*x53
     1  +coeff( 28)        *x33        
     2  +coeff( 29)    *x21*x33        
     3  +coeff( 30)        *x31*x42*x51
     4  +coeff( 31)    *x21*x32*x41*x51
     5  +coeff( 32)    *x23    *x41*x51
     6  +coeff( 33)    *x22    *x41*x52
     7  +coeff( 34)        *x31    *x54
     8  +coeff( 35)        *x35    *x51
      y_x2_fp     =y_x2_fp     
     9  +coeff( 36)    *x21*x32*x41*x52
     1  +coeff( 37)                *x51
     2  +coeff( 38)    *x21    *x43    
     3  +coeff( 39)            *x43*x51
     4  +coeff( 40)*x11        *x41*x51
     5  +coeff( 41)    *x21    *x41*x52
     6  +coeff( 42)    *x22*x32*x41    
     7  +coeff( 43)    *x22*x31*x42    
     8  +coeff( 44)    *x22    *x43    
      y_x2_fp     =y_x2_fp     
     9  +coeff( 45)        *x33    *x52
     1  +coeff( 46)        *x31*x42*x52
     2  +coeff( 47)    *x24    *x41    
     3  +coeff( 48)    *x21*x31    *x53
     4  +coeff( 49)    *x21    *x41*x53
     5  +coeff( 50)            *x41*x54
     6  +coeff( 51)        *x34*x41*x51
     7  +coeff( 52)    *x21*x31*x42*x52
     8  +coeff( 53)        *x33    *x53
      y_x2_fp     =y_x2_fp     
     9  +coeff( 54)*x11*x22*x31    *x51
     1  +coeff( 55)    *x22*x33    *x52
     2  +coeff( 56)    *x23*x35        
     3  +coeff( 57)        *x32*x41*x51
     4  +coeff( 58)        *x32*x43    
     5  +coeff( 59)    *x21*x31*x42*x51
     6  +coeff( 60)    *x21    *x43*x51
     7  +coeff( 61)*x11*x22    *x41    
     8  +coeff( 62)    *x23*x31    *x51
      y_x2_fp     =y_x2_fp     
     9  +coeff( 63)    *x22*x31    *x52
     1  +coeff( 64)    *x21*x34*x41    
     2  +coeff( 65)    *x23*x32*x41    
     3  +coeff( 66)    *x22*x33    *x51
     4  +coeff( 67)    *x23*x33    *x51
     5  +coeff( 68)    *x21*x32*x43*x52
     6  +coeff( 69)*x11*x24*x31    *x51
     7  +coeff( 70)    *x22            
     8  +coeff( 71)    *x21        *x51
      y_x2_fp     =y_x2_fp     
     9  +coeff( 72)    *x22        *x51
     1  +coeff( 73)    *x21*x31*x42    
     2  +coeff( 74)*x11*x21*x31        
     3  +coeff( 75)    *x23*x31        
     4  +coeff( 76)*x11    *x31    *x51
     5  +coeff( 77)        *x35        
     6  +coeff( 78)        *x31*x44    
     7  +coeff( 79)            *x43*x52
     8  +coeff( 80)*x11*x21*x31    *x51
      y_x2_fp     =y_x2_fp     
     9  +coeff( 81)    *x21*x35        
     1  +coeff( 82)    *x23*x33        
     2  +coeff( 83)    *x22*x31*x42*x51
     3  +coeff( 84)    *x21*x33    *x52
     4  +coeff( 85)*x12        *x42    
     5  +coeff( 86)*x11*x22    *x41*x51
     6  +coeff( 87)    *x22*x31    *x53
     7  +coeff( 88)    *x21*x31    *x54
     8  +coeff( 89)    *x21*x34*x41*x51
      y_x2_fp     =y_x2_fp     
     9  +coeff( 90)        *x34*x41*x52
     1  +coeff( 91)*x11*x21*x33    *x51
     2  +coeff( 92)    *x22    *x43*x52
     3  +coeff( 93)    *x21    *x43*x53
     4  +coeff( 94)    *x22    *x41*x54
     5  +coeff( 95)    *x22*x34*x41*x51
c
      return
      end
      function p_x2_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 91)
      data ncoeff/ 90/
      data avdat/  0.2196344E-03/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     + -0.19247606E-03,-0.30165680E-04, 0.16618835E-01,-0.33575770E-01,
     +  0.59792213E-02, 0.16139746E-01,-0.11178758E-01,-0.10163613E-01,
     +  0.61485847E-02, 0.17383876E-02, 0.84913295E-03, 0.39883023E-02,
     +  0.19369408E-02, 0.32752035E-02, 0.14064120E-02, 0.29731516E-03,
     + -0.32807780E-04, 0.63954369E-03,-0.28680597E-03, 0.23024317E-02,
     + -0.68208243E-03, 0.31436022E-04,-0.10014443E-03, 0.47689694E-03,
     + -0.10665732E-02,-0.66188094E-03, 0.19328165E-02, 0.10972802E-02,
     +  0.74816262E-03,-0.58737048E-03,-0.79633479E-04,-0.15594561E-03,
     + -0.15043474E-02,-0.61121839E-03,-0.20683093E-03, 0.42261272E-02,
     +  0.50478464E-03, 0.78611322E-06,-0.77213437E-04, 0.15843919E-03,
     +  0.11704338E-03, 0.58864127E-03,-0.13246690E-03, 0.52982970E-03,
     + -0.14078873E-03,-0.99488789E-04,-0.11321964E-02,-0.22779974E-03,
     +  0.52329469E-04,-0.89309883E-03,-0.52607898E-03,-0.68552117E-03,
     + -0.37006255E-04, 0.98339422E-03, 0.65866212E-03, 0.24396670E-02,
     +  0.36673259E-03,-0.44283070E-03, 0.54640364E-03, 0.14600269E-04,
     +  0.85503457E-03,-0.63520286E-03,-0.45840099E-03, 0.20366987E-03,
     + -0.40659201E-03,-0.17443770E-04,-0.18554163E-04,-0.69404073E-03,
     + -0.98915087E-04, 0.47473633E-03, 0.19318609E-03, 0.49849419E-03,
     +  0.77905883E-04,-0.16408801E-03, 0.40580024E-03, 0.11818596E-03,
     + -0.90256834E-03, 0.15120785E-02, 0.33168474E-03,-0.82049601E-05,
     + -0.62369321E-06,-0.91499360E-05,-0.65904278E-04,-0.27207885E-04,
     + -0.15975360E-04, 0.23506873E-03, 0.18569137E-03, 0.49771792E-04,
     + -0.10858798E-03,-0.57668908E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      p_x2_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)        *x31        
     4  +coeff(  4)            *x41    
     5  +coeff(  5)    *x21*x31        
     6  +coeff(  6)    *x21    *x41    
     7  +coeff(  7)        *x31    *x51
     8  +coeff(  8)            *x41*x51
      p_x2_fp     =p_x2_fp     
     9  +coeff(  9)    *x22    *x41    
     1  +coeff( 10)        *x31    *x52
     2  +coeff( 11)    *x24*x31        
     3  +coeff( 12)    *x22*x31        
     4  +coeff( 13)            *x41*x52
     5  +coeff( 14)        *x31*x42    
     6  +coeff( 15)            *x43    
     7  +coeff( 16)    *x21*x31    *x51
     8  +coeff( 17)        *x34*x41    
      p_x2_fp     =p_x2_fp     
     9  +coeff( 18)        *x33        
     1  +coeff( 19)*x11        *x41    
     2  +coeff( 20)        *x32*x41    
     3  +coeff( 21)        *x31    *x53
     4  +coeff( 22)                *x51
     5  +coeff( 23)*x11    *x31        
     6  +coeff( 24)    *x21    *x41*x51
     7  +coeff( 25)    *x23    *x41    
     8  +coeff( 26)            *x41*x53
      p_x2_fp     =p_x2_fp     
     9  +coeff( 27)    *x22    *x43    
     1  +coeff( 28)        *x31    *x54
     2  +coeff( 29)    *x22    *x43*x52
     3  +coeff( 30)    *x23*x31        
     4  +coeff( 31)    *x21*x33        
     5  +coeff( 32)*x11*x21    *x41    
     6  +coeff( 33)    *x21*x31*x42    
     7  +coeff( 34)    *x21    *x43    
     8  +coeff( 35)        *x33    *x51
      p_x2_fp     =p_x2_fp     
     9  +coeff( 36)    *x22*x31*x42    
     1  +coeff( 37)    *x21*x31    *x53
     2  +coeff( 38)*x11*x23*x31        
     3  +coeff( 39)    *x23*x32*x41    
     4  +coeff( 40)    *x21*x34*x41    
     5  +coeff( 41)        *x34*x41*x51
     6  +coeff( 42)    *x22    *x43*x51
     7  +coeff( 43)        *x31*x44*x51
     8  +coeff( 44)    *x24*x32*x41    
      p_x2_fp     =p_x2_fp     
     9  +coeff( 45)    *x22*x33    *x52
     1  +coeff( 46)*x11*x21*x31        
     2  +coeff( 47)    *x21*x32*x41    
     3  +coeff( 48)    *x22*x31    *x51
     4  +coeff( 49)*x11        *x41*x51
     5  +coeff( 50)        *x32*x41*x51
     6  +coeff( 51)        *x31*x42*x51
     7  +coeff( 52)            *x43*x51
     8  +coeff( 53)    *x21*x31    *x52
      p_x2_fp     =p_x2_fp     
     9  +coeff( 54)    *x22*x33        
     1  +coeff( 55)    *x24    *x41    
     2  +coeff( 56)    *x22*x32*x41    
     3  +coeff( 57)    *x23*x31    *x51
     4  +coeff( 58)    *x21*x32*x41*x51
     5  +coeff( 59)        *x33    *x52
     6  +coeff( 60)    *x22    *x41*x52
     7  +coeff( 61)        *x31*x42*x52
     8  +coeff( 62)    *x21    *x41*x53
      p_x2_fp     =p_x2_fp     
     9  +coeff( 63)            *x41*x54
     1  +coeff( 64)    *x24    *x41*x51
     2  +coeff( 65)        *x33    *x53
     3  +coeff( 66)    *x22            
     4  +coeff( 67)*x11    *x31    *x51
     5  +coeff( 68)    *x22    *x41*x51
     6  +coeff( 69)    *x21    *x41*x52
     7  +coeff( 70)    *x21*x31*x42*x51
     8  +coeff( 71)    *x21    *x43*x51
      p_x2_fp     =p_x2_fp     
     9  +coeff( 72)    *x22*x31    *x52
     1  +coeff( 73)    *x23*x33        
     2  +coeff( 74)    *x21*x32*x41*x52
     3  +coeff( 75)    *x21*x31    *x54
     4  +coeff( 76)    *x21    *x41*x54
     5  +coeff( 77)    *x23    *x43*x52
     6  +coeff( 78)    *x21*x32*x43*x52
     7  +coeff( 79)    *x24*x31    *x53
     8  +coeff( 80)            *x42    
      p_x2_fp     =p_x2_fp     
     9  +coeff( 81)                *x52
     1  +coeff( 82)*x11    *x31*x41    
     2  +coeff( 83)    *x22*x31*x41    
     3  +coeff( 84)    *x22    *x42    
     4  +coeff( 85)        *x32    *x52
     5  +coeff( 86)        *x33*x42    
     6  +coeff( 87)        *x31*x44    
     7  +coeff( 88)    *x21*x31*x41*x52
     8  +coeff( 89)            *x43*x52
      p_x2_fp     =p_x2_fp     
     9  +coeff( 90)    *x23*x31*x42    
c
      return
      end
      function l_x2_fp     (x,m)
      dimension x(m)
      dimension xmin(10),xmax(10),scale(10),xmean(10)
      dimension coeff( 51)
      data ncoeff/ 50/
      data avdat/ -0.2804277E-01/
      data xmin/
     1 -0.19989E-02,-0.60671E-01,-0.74953E-01,-0.32916E-01,-0.48949E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data xmax/
     1  0.19982E-02, 0.64590E-01, 0.74761E-01, 0.32818E-01, 0.49987E-01,
     2  0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00, 0.00000E+00/
      data scale /10*0./
      data coeff/
     +  0.17539877E-01,-0.32732999E+00,-0.31988200E-01, 0.61839488E-02,
     + -0.39297767E-01,-0.16090857E-01,-0.84665399E-02,-0.45203017E-02,
     + -0.26476302E-03, 0.51656528E-02,-0.52165622E-02, 0.12420099E-02,
     +  0.94227716E-02, 0.29172925E-02, 0.71176854E-02,-0.53125736E-02,
     + -0.15457046E-02, 0.37932212E-02, 0.83981352E-02,-0.60622342E-03,
     + -0.22106476E-02, 0.23280368E-02, 0.30429251E-02,-0.11501996E-02,
     +  0.40463492E-03, 0.35656593E-02, 0.30201916E-02,-0.11274476E-02,
     + -0.62607118E-03, 0.12246178E-02, 0.11895770E-01, 0.10518112E-01,
     +  0.15079837E-02,-0.18096036E-03, 0.91663009E-03,-0.10121992E-02,
     + -0.40625266E-03, 0.21294005E-01, 0.17504724E-01, 0.17675704E-02,
     + -0.22992280E-02,-0.32825879E-03, 0.21171642E-01,-0.32848585E-02,
     +  0.26577590E-02, 0.15666885E-01,-0.24000383E-02, 0.17043669E-01,
     + -0.18115180E-01,-0.29672720E-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
      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
      x51 = x5
      x52 = x51*x5
      x53 = x52*x5
      x54 = x53*x5
c
c                  function
c
      l_x2_fp     =avdat
     1  +coeff(  1)                    
     2  +coeff(  2)    *x21            
     3  +coeff(  3)                *x51
     4  +coeff(  4)*x11                
     5  +coeff(  5)    *x22            
     6  +coeff(  6)                *x52
     7  +coeff(  7)            *x42    
     8  +coeff(  8)    *x21        *x51
      l_x2_fp     =l_x2_fp     
     9  +coeff(  9)        *x34        
     1  +coeff( 10)    *x21    *x44    
     2  +coeff( 11)        *x31*x41    
     3  +coeff( 12)*x11*x21            
     4  +coeff( 13)    *x21*x31*x41    
     5  +coeff( 14)                *x53
     6  +coeff( 15)    *x23*x32        
     7  +coeff( 16)        *x32        
     8  +coeff( 17)    *x23            
      l_x2_fp     =l_x2_fp     
     9  +coeff( 18)    *x21*x32        
     1  +coeff( 19)    *x21    *x42    
     2  +coeff( 20)*x11            *x51
     3  +coeff( 21)    *x22        *x51
     4  +coeff( 22)        *x32    *x51
     5  +coeff( 23)        *x31*x41*x51
     6  +coeff( 24)    *x21        *x52
     7  +coeff( 25)*x11*x22            
     8  +coeff( 26)    *x22*x31*x41    
      l_x2_fp     =l_x2_fp     
     9  +coeff( 27)    *x22    *x42    
     1  +coeff( 28)    *x23        *x51
     2  +coeff( 29)    *x21*x32    *x51
     3  +coeff( 30)    *x21    *x42*x51
     4  +coeff( 31)    *x23*x31*x41    
     5  +coeff( 32)    *x23    *x42    
     6  +coeff( 33)    *x24*x32        
     7  +coeff( 34)            *x42*x51
     8  +coeff( 35)    *x24            
      l_x2_fp     =l_x2_fp     
     9  +coeff( 36)        *x31*x41*x52
     1  +coeff( 37)                *x54
     2  +coeff( 38)    *x21*x32*x42    
     3  +coeff( 39)    *x21*x31*x43    
     4  +coeff( 40)    *x22*x31*x43    
     5  +coeff( 41)    *x23*x31*x41*x51
     6  +coeff( 42)    *x22*x32    *x52
     7  +coeff( 43)    *x23*x33*x41    
     8  +coeff( 44)    *x23*x32*x42    
      l_x2_fp     =l_x2_fp     
     9  +coeff( 45)    *x21*x34    *x52
     1  +coeff( 46)    *x21*x33*x41*x52
     2  +coeff( 47)    *x21*x31*x41*x54
     3  +coeff( 48)    *x23*x34*x42    
     4  +coeff( 49)    *x23*x33*x41*x52
     5  +coeff( 50)    *x23*x32    *x54
c
      return
      end
