c
c
      program qOPairy
c
c
      external qquad
      real*16 x0,x1,a1(320),b1(320),a2(320),b2(320),
     *e(320),z(83),w(83),q1mach,qepsma,dendl(3),dendr(3),
     *qeps,dxp(1),dyp(1),dxfer(320),dwfer(320),dalpha(40),
     *dbeta(40),dbe(40),dx(320),dw(320),dxm(960),dwm(960),
     *dp0(960),dp1(960),dp2(960)
      logical finld,finrd
      common/q/x0,x1,a1,b1,a2,b2,e
      common/qa/z,w
c
c This routine generates in quadruple precision the first 40
c recurrence coefficients of the orthogonal polynomials for 
c the Airy weight function. Printed are the quadruple-precision
c values of the alpha- and beta-coefficients. The routine calls
c upon the quadruple-precision analogues of the ORTHPOL routines
c recur, r1mach, mcdis.
c
      nmax=40
      x0=2.q0
      x1=6.q0
      call qrecur(320,1,0.q0,0.q0,a1,b1,ier1)
      open(unit=5,file='qablee320',status='old')
      read(*,1) (a2(k),b2(k),k=1,320)
    1 format(1x,2d37.29)
      open(unit=5,file='qzwlag',status='old')
      read(*,2) (z(k),w(k),k=1,83)
    2 format(1x,2d38.30)
      finld=.true.
      finrd=.false.
      qepsma=q1mach(3)
c
c qepsma is the machine quadruple precision.
c
      iq=1
      idelta=1
      irout=1
      mcd=3
      mp=0
      ncpmd=320
c
c Set up the partition for the discretization of the inner product.
c
      dendl(1)=0.q0
      dendr(1)=2.q0
      dendl(2)=2.q0
      dendr(2)=6.q0
      dendl(3)=6.q0
      qeps=1.q5*qepsma
c
c Compute the desired recursion coefficients by the multiple-component
c discretization procedure. 
c
      call qmcdis(nmax,ncpmd,mcd,mp,dxp,dyp,qquad,qeps,iq,idelta,irout,
     *  finld,finrd,dendl,dendr,dxfer,dwfer,dalpha,dbeta,ncapd,kountd,
     *  ierrd,ied,dbe,dx,dw,dxm,dwm,dp0,dp1,dp2)
c
c Print the results.
c
      write(*,3)
    3 format(/5x,'k',7x,'dalpha(k)',14x,'dbeta(k)'/)
      do 10 k=1,nmax
        km1=k-1
        write(*,4) km1,dalpha(k),dbeta(k)
    4   format(1x,i5,2d36.28)
   10 continue
      stop
      end 

      subroutine qquad(n,dx,dw,i,ierr)
      real*16 dx(n),dw(n),x0,x1,a1(320),b1(320),
     *a2(320),b2(320),e(320),z(83),w(83),pi,c,c1,
     *qgamma,eps,q1mach,x,qairy
      common/q/x0,x1,a1,b1,a2,b2,e
      common/qa/z,w
      pi=4.q0*qatan(1.q0)
      c=.5q0**(2.q0/3.q0)
      c1=(2.q0**(2.q0/3.q0))*pi/((3.q0**(1.q0/6.q0))*
     *  qgamma(2.q0/3.q0,ie))
      eps=q1mach(3)
      if(i.eq.1) then
        call qgauss(n,a1,b1,eps,dx,dw,ierg,e)
        do 10 k=1,n
          x=.5q0*x0*(1.q0+dx(k))
          dw(k)=.5q0*x0*c1*qairy(c*(x**2))*qexp(
     *      -(x**3)/3.q0)*dw(k)
          dx(k)=x**3/3.q0
   10   continue
      else if(i.eq.2) then
        call qgauss(n,a1,b1,eps,dx,dw,ierg,e)
        do 20 k=1,n
          x=.5q0*((x1-x0)*dx(k)+x1+x0)
          dw(k)=.5q0*(x1-x0)*c1*qairy(c*(x**2))*qexp(
     *      -(x**3)/3.q0)*dw(k)
          dx(k)=x**3/3.q0
   20   continue
      else
        call qgauss(n,a2,b2,eps,dx,dw,ierg,e)
        do 30 k=1,n
          x=x1+dx(k)
          dw(k)=c1*qexp(-(x1**3)/3.q0)*qairy(c*(x**2))*
     *      qexp(-x1*(x1*dx(k)+dx(k)**2))*dw(k)
          dx(k)=(x**3)/3.q0
   30 continue
      end if
      return
      end

      real*16 function qairy(x)
      real*16 z(83),w(83),x,x1,pi,c1,c2,c3,
     *series,qgamma,sum
      common/qa/z,w
      if(x.le.0.q0) then
        write(*,1)
    1   format(1x,'invalid x-argument in qairy')
        stop
      end if
      pi=4.q0*qatan(1.q0)
      c1=1.q0/((3.q0**(2.q0/3.q0))*qgamma(2.q0/3.q0,ie))
      c2=1.q0/((3.q0**(4.q0/3.q0))*qgamma(4.q0/3.q0,ie))
      if(x.le.3.q0**(2.q0/3.q0)) then
        qairy=c1*series(x,-1.q0/3.q0)-c2*x*series(x,1.q0/3.q0)
        return
      else
        n=83
        c3=1.q0/(qsqrt(pi)*(48.q0**(1.q0/6.q0))*
     *    qgamma(5.q0/6.q0,ie))
        sum=0.q0
        x1=2.q0*((qsqrt(x))**3)/3.q0
        do 10 k=1,n
          sum=sum+w(k)*(2.q0+z(k)/x1)**(-1.q0/6.q0)
   10   continue
        qairy=c3*qexp(-x1)*(x1**(-1.q0/6.q0))*sum
      end if
      return
      end

      real*16 function series(x,qnu)
      real*16 x,qnu,x3,sum1,t1,qk,sum,t
      k=0
      x3=x**3
      sum1=1.q0
      t1=1.q0
   10 k=k+1
      qk=qreal(k)
      sum=sum1
      t=t1
      t1=x3*t/(9.q0*qk*(qnu+qk))
      sum1=sum+t1
      if(sum1.ne.sum) goto 10
      series=sum1
      return
      end
  
      real*16 function qwf(x,i)
      real*16 x
      return
      end
c
c
                              file qablee320
                 (containing the recursion coefficients
               for the weight function w(x)=exp(-x**3/3)
c
c
   0.72901113294722698141863626470D+00  0.12878993168540690872006831600D+01
   0.10422198256747440911363056128D+01  0.24500097941742094234842574778D+00
   0.12537306422019648129872014375D+01  0.35307351727990706406124026858D+00
   0.14061820889340038794916026341D+01  0.45380654475472013456007795220D+00
   0.15304717088698266166121572948D+01  0.54670915163293606036098161021D+00
   0.16371146876931010219364478718D+01  0.63279143126565640324575431662D+00
   0.17313265280009313616975788254D+01  0.71359155024155923743765162023D+00
   0.18162157284093989632512284557D+01  0.79017160081817896753892393209D+00
   0.18938033162945060603236080352D+01  0.86327669550039949557728853867D+00
   0.19654868263312374176204257310D+01  0.93345297398370756848555274697D+00
   0.20322783394582393893632309784D+01  0.10011143019264015566512732609D+01
   0.20949374105669606086639947880D+01  0.10665830763064051841449652052D+01
   0.21540505128026897945531915917D+01  0.11301163055170059928542543211D+01
   0.22100811161203423821005766708D+01  0.11919228879829774855536020259D+01
   0.22634026387069417553930361631D+01  0.12521754391488299058313985041D+01
   0.23143209101304160436902296711D+01  0.13110186225966271500022586854D+01
   0.23630899469916816496199493530D+01  0.13685751676146236967777835370D+01
   0.24099233167562483428250753554D+01  0.14249503121802225111773503609D+01
   0.24550025049603012817942699756D+01  0.14802351473683328458639907371D+01
   0.24984831949819342537544238777D+01  0.15345091782276951177983174137D+01
   0.25405000612283119886343859922D+01  0.15878423147510053923208124784D+01
   0.25811704827847374855997795188D+01  0.16402964409812249078988145264D+01
   0.26205974593929454159753483323D+01  0.16919266668748729763938073675D+01
   0.26588719287966830624731312932D+01  0.17427823381753443922811013891D+01
   0.26960746284907256780896244087D+01  0.17929078592990917007414275534D+01
   0.27322776063048205311543683532D+01  0.18423433700265341465386347559D+01
   0.27675454571717529825350311567D+01  0.18911253066553886888308721438D+01
   0.28019363441248307021832069394D+01  0.19392868709402412502509847931D+01
   0.28355028476091342383228373014D+01  0.19868584247626864617226799548D+01
   0.28682926769586859715160422938D+01  0.20338678244810629338780896435D+01
   0.29003492702998949393655050712D+01  0.20803407059072388143606695077D+01
   0.29317123034448058398525621051D+01  0.21263007285787534405587402412D+01
   0.29624181240177932712046692203D+01  0.21717697862468153545210753620D+01
   0.29925001237515478508316567445D+01  0.22167681891479147815272905429D+01
   0.30219890593323144507492897823D+01  0.22613148225707880517364009382D+01
   0.30509133301826576729319198050D+01  0.23054272853994264791569129249D+01
   0.30792992200056565662026822862D+01  0.23491220116538714297358240852D+01
   0.31071711076765261931000936492D+01  0.23924143775243274223168909095D+01
   0.31345516520812322083528718800D+01  0.24353187959711047733893880535D+01
   0.31614619547104344737770880348D+01  0.24778488006206924407801133331D+01
   0.31879217031784896668469767058D+01  0.25200171204097750101451189831D+01
   0.32139492983187900444063269085D+01  0.25618357462011029676474818209D+01
   0.32395619670834674350353415975D+01  0.26033159904076207984155255131D+01
   0.32647758631281435653843488033D+01  0.26444685405062113226316180429D+01
   0.32896061566759224289324702383D+01  0.26853035071935925689466141477D+01
   0.33140671150174094039733706216D+01  0.27258304678293812119728773869D+01
   0.33381721748059084524232658575D+01  0.27660585057212022564344507508D+01
   0.33619340071417189318168347544D+01  0.28059962457308532772267956978D+01
   0.33853645763007463348135812536D+01  0.28456518866164136221120016917D+01
   0.34084751928457419094411058459D+01  0.28850332304707991396648148310D+01
   0.34312765617595960376041941071D+01  0.29241477095709586338158508017D+01
   0.34537788261561544431768330992D+01  0.29630024109123511349908536057D+01
   0.34759916070525041412493462378D+01  0.30016040986694381331530614433D+01
   0.34979240396255451337072365690D+01  0.30399592347937722454236587831D+01
   0.35195848063232490969243473417D+01  0.30780739979361209012368856335D+01
   0.35409821671559283152598001499D+01  0.31159543008573155256570004599D+01
   0.35621239874539566537727945795D+01  0.31536058064736520924642133963D+01
   0.35830177633447533951473732170D+01  0.31910339426662611998490758405D+01
   0.36036706451726725775312274834D+01  0.32282439159695578250450222160D+01
   0.36240894590600769341402307676D+01  0.32652407242413726126893945141D+01
   0.36442807267857649769128023220D+01  0.33020291684064051108580312619D+01
   0.36642506841375982734441561159D+01  0.33386138633550110670348501655D+01
   0.36840052978792522322969286686D+01  0.33749992480708595369493872675D+01
   0.37035502814561566951076850089D+01  0.34111895950535170793750245271D+01
   0.37228911095526214764804371042D+01  0.34471890190954043541709414041D+01
   0.37420330316006177540628824937D+01  0.34830014854667127140132933819D+01
   0.37609810843305050698293014723D+01  0.35186308175566686331790121413D+01
   0.37797401034449813458563829941D+01  0.35540807040149093464079616204D+01
   0.37983147344895402934771362924D+01  0.35893547054326126414826570336D+01
   0.38167094429856184599796970890D+01  0.36244562605993458579220530902D+01
   0.38349285238862922632976205839D+01  0.36593886923683105704849373538D+01
   0.38529761104087481765375391123D+01  0.36941552131597140246024573812D+01
   0.38708561822927140715147774030D+01  0.37287589301293560275112665007D+01
   0.38885725735295349252171175067D+01  0.37632028500271457480886523326D+01
   0.39061289796025393397789679483D+01  0.37974898837681262529448422773D+01
   0.39235289642757203826403084601D+01  0.38316228507366589487839642055D+01
   0.39407759659644977960027091361D+01  0.38656044828426820583204465028D+01
   0.39578733037193972884408295863D+01  0.38994374283473863329246854821D+01
   0.39748241828508401824577094002D+01  0.39331242554742293784653419571D+01
   0.39916317002208513705669331421D+01  0.39666674558199213686993360922D+01
   0.40082988492253374003300086911D+01  0.40000694475788455386635810099D+01
   0.40248285244886349789028024118D+01  0.40333325785933143124796793901D+01
   0.40412235262902615794070290823D+01  0.40664591292410952705084454306D+01
   0.40574865647421949958482378802D+01  0.40994513151707606872443400964D+01
   0.40736202637335506906816180474D+01  0.41323112898946114495650073236D+01
   0.40896271646581996011693246908D+01  0.41650411472481931258417576874D+01
   0.41055097299396614002625581596D+01  0.41976429237247519702085102240D+01
   0.41212703463665072110875996604D+01  0.42301186006923656220439145292D+01
   0.41369113282505009205080260076D+01  0.42624701065009217588737063834D+01
   0.41524349204187901424905469280D+01  0.42946993184856031170045878911D+01
   0.41678433010506181722466485273D+01  0.43268080648730647523791806788D+01
   0.41831385843682594650891910268D+01  0.43587981265960552688739027230D+01
   0.41983228231911765753530940597D+01  0.43906712390218344864044532151D+01
   0.42133980113617501037031731113D+01  0.44224290935993725052689091234D+01
   0.42283660860503396468925787233D+01  0.44540733394299765100255640389D+01
   0.42432289299468881922710002829D+01  0.44856055847656793905525085402D+01
   0.42579883733457805072449119094D+01  0.45170273984394360347965368168D+01
   0.42726461961302039323580779280D+01  0.45483403112309068855953545354D+01
   0.42872041296618340729975580636D+01  0.45795458171713621710917781131D+01
   0.43016638585812750209366629307D+01  0.46106453747910124116294116296D+01
   0.43160270225243210528448621246D+01  0.46416404083118598327344156830D+01
   0.43302952177587716536340216100D+01  0.46725323087889697745840768213D+01
   0.43444699987462218495989232233D+01  0.47033224352028798129494590022D+01
   0.43585528796329630821247269162D+01  0.47340121155056959406273093583D+01
   0.43725453356738642778034909409D+01  0.47646026476232687518362566125D+01
   0.43864488045928566236365014500D+01  0.47950953004156971689340653531D+01
   0.44002646878834172451457843169D+01  0.48254913145982719798780223649D+01
   0.44139943520522350632750540747D+01  0.48557919036248455213896588274D+01
   0.44276391298090452553225034466D+01  0.48859982545354965212442356457D+01
   0.44412003212054357652809444846D+01  0.49161115287702497404863411281D+01
   0.44546791947252591050275378360D+01  0.49461328629505080263163755231D+01
   0.44680769883291242604893010421D+01  0.49760633696297591438705619574D+01
   0.44813949104552959533023105207D+01  0.50059041380150307916650814110D+01
   0.44946341409791909739081763281D+01  0.50356562346604840549062027378D+01
   0.45077958321335330327192584121D+01  0.50653207041344577854163769310D+01
   0.45208811093911078724764205490D+01  0.50948985696612036237702835097D+01
   0.45338910723119486061516366785D+01  0.51243908337384832373663838782D+01
   0.45468267953566768027011826211D+01  0.51537984787321355055551657812D+01
   0.45596893286676271978947057084D+01  0.51831224674486615340119970744D+01
   0.45724796988192925634923608767D+01  0.52123637436868192438083493193D+01
   0.45851989095395397693869543792D+01  0.52415232327691665966589753450D+01
   0.45978479424029680006175642941D+01  0.52706018420544430473010517392D+01
   0.46104277574977050583242219133D+01  0.52996004614316323360125140114D+01
   0.46229392940668673250696477323D+01  0.53285199637965060448834762255D+01
   0.46353834711258429825958019898D+01  0.53573612055114062520499433694D+01
   0.46477611880564961315008052231D+01  0.53861250268489869542817858710D+01
   0.46600733251793312982512983225D+01  0.54148122524205975286828345297D+01
   0.46723207443046031674118115151D+01  0.54434236915899572193387562987D+01
   0.46845042892633050075141197146D+01  0.54719601388727373260309275497D+01
   0.46966247864189209470404892648D+01  0.55004223743226373111953112067D+01
   0.47086830451607817984157538475D+01  0.55288111639045123089925300162D+01
   0.47206798583798213336413229370D+01  0.55571272598550824060559264062D+01
   0.47326160029274896100720378708D+01  0.55853714010317284644519876199D+01
   0.47444922400585419663985555019D+01  0.56135443132498550781275525629D+01
   0.47563093158583865064531768365D+01  0.56416467096092784058256456389D+01
   0.47680679616556391221118824148D+01  0.56696792908100750234961043301D+01
   0.47797688944205032463837172049D+01  0.56976427454583075106806691306D+01
   0.47914128171495614530058813536D+01  0.57255377503620231565574294792D+01
   0.48030004192375376172148982302D+01  0.57533649708179038755024507843D+01
   0.48145323768365615193767335871D+01  0.57811250608889280968921848096D+01
   0.48260093532034424115905000965D+01  0.58088186636733889813285488007D+01
   0.48374319990354340866660608026D+01  0.58364464115655977612859356738D+01
   0.48488009527949513046238200306D+01  0.58640089265085862577064914621D+01
   0.48601168410236759653962072598D+01  0.58915068202391086379810138910D+01
   0.48713802786464710943346006350D+01  0.59189406945252292108059068263D+01
   0.48825918692655014609517275814D+01  0.59463111413967704582228291043D+01
   0.48937522054449414171077705228D+01  0.59736187433688835460073099888D+01
   0.49048618689866332588392369087D+01  0.60008640736589921942273142744D+01
   0.49159214311970430304008657008D+01  0.60280476963973499962920708701D+01
   0.49269314531458451476371141209D+01  0.60551701668314410153426845742D+01
   0.49378924859164524716844589909D+01  0.60822320315244437315760790496D+01
   0.49488050708487944675226267138D+01  0.61092338285479691350724556246D+01
   0.49596697397746327922461544972D+01  0.61361760876692749296697233181D+01
   0.49704870152456910350288853988D+01  0.61630593305331494097613521643D+01
   0.49812574107548633370341414244D+01  0.61898840708386505704524740635D+01
   0.49919814309507552197593550407D+01  0.62166508145108783905584754858D+01
   0.50026595718457991114581400068D+01  0.62433600598679509670432944422D+01
   0.50132923210181767523540786726D+01  0.62700122977833482594640428396D+01
   0.50238801578077708512468219750D+01  0.62966080118437806057475049017D+01
   0.50344235535063590314818116533D+01  0.63231476785027328791713201961D+01
   0.50449229715422542174326176192D+01  0.63496317672298291547536196737D+01
   0.50553788676595871494956563891D+01  0.63760607406561570262972158046D+01
   0.50657916900924186534288278239D+01  0.64024350547156852488883696300D+01
   0.50761618797338616073336130244D+01  0.64287551587829031623328524053D+01
   0.50864898703003852265992446010D+01  0.64550214958068053662059205920D+01
   0.50967760884914673047905892247D+01  0.64812345024413403549974625087D+01
   0.51070209541447533889635506242D+01  0.65073946091724372710207371018D+01
   0.51172248803868755144600833984D+01  0.65335022404417205827250563345D+01
   0.51273882737800770610624876859D+01  0.65595578147670183368032134589D+01
   0.51375115344647845045718713678D+01  0.65855617448597656545623521309D+01
   0.51475950562982613113651571574D+01  0.66115144377394013375024995684D+01
   0.51576392269894739450207365097D+01  0.66374162948448518054799727340D+01
   0.51676444282302949111727359775D+01  0.66632677121431931052432115203D+01
   0.51776110358231629475456711457D+01  0.66890690802355783899686742537D+01
   0.51875394198053158594797736436D+01  0.67148207844605150745533027251D+01
   0.51974299445697070966441049074D+01  0.67405232049945728100793375718D+01
   0.52072829689827129540907689117D+01  0.67661767169506004876585481228D+01
   0.52170988464987332509146365435D+01  0.67917816904735276707267291962D+01
   0.52268779252717844836463440931D+01  0.68173384908338231600552342666D+01
   0.52366205482641807607005590278D+01  0.68428474785186808118346084879D+01
   0.52463270533523942907579825966D+01  0.68683090093210002510095970015D+01
   0.52559977734301838143356704226D+01  0.68937234344262277447155237394D+01
   0.52656330365090761268544186735D+01  0.69190911004971202195456754977D+01
   0.52752331658162827364802348240D+01  0.69444123497564932170678392050D+01
   0.52847984798901307244910548732D+01  0.69696875200680114803296257950D+01
   0.52943292926730840238268350034D+01  0.69949169450150788460823429128D+01
   0.53038259136024285970638305324D+01  0.70201009539778821793473906250D+01
   0.53132886476986923728530777083D+01  0.70452398722086422251715728618D+01
   0.53227177956518682847005976374D+01  0.70703340209051224635707355600D+01
   0.53321136539055063429295399189D+01  0.70953837172824453345156360561D+01
   0.53414765147387383550894818245D+01  0.71203892746432635473000745196D+01
   0.53508066663462966875387272618D+01  0.71453510024463325998298634614D+01
   0.53601043929165863272181265486D+01  0.71702692063735291055054510827D+01
   0.53693699747078674537668279593D+01  0.71951441883953580557990531221D+01
   0.53786036881226037643077875355D+01  0.72199762468349907328347091120D+01
   0.53878058057800299028483608821D+01  0.72447656764308736258733144228D+01
   0.53969765965869895298714248601D+01  0.72695127683979473963053738739D+01
   0.54061163258070938220761161628D+01  0.72942178104875136753913673593D+01
   0.54152252551282485142645806374D+01  0.73188810870457862654953307279D+01
   0.54243036427285959821129749525D+01  0.73435028790711621469610252223D+01
   0.54333517433409173132048979803D+01  0.73680834642702465672037059361D+01
   0.54423698083155378215720958462D+01  0.73926231171126654042433697283D+01
   0.54513580856817780255372799725D+01  0.74171221088846969520798911372D+01
   0.54603168202079907274650451099D+01  0.74415807077417542683778815485D+01
   0.54692462534602235047918789279D+01  0.74659991787597482543354801773D+01
   0.54781466238595446422267188851D+01  0.74903777839853607008736755744D+01
   0.54870181667380693031942927044D+01  0.75147167824852556329856796127D+01
   0.54958611143937215524373613746D+01  0.75390164303942564138779148747D+01
   0.55046756961437666992966559522D+01  0.75632769809625152311247851498D+01
   0.55134621383771473307326057770D+01  0.75874986846017007772158820462D+01
   0.55222206646056553429084150193D+01  0.76116817889302291554189364602D+01
   0.55309514955139712584666166738D+01  0.76358265388175622876892397697D+01
   0.55396548490086011319233185343D+01  0.76599331764275973733507558651D+01
   0.55483309402657403963700138219D+01  0.76840019412611702444273514401D+01
   0.55569799817780930894742037241D+01  0.77080330701976947848311356133D+01
   0.55656021834006740142340173564D+01  0.77320267975359599251783139725D+01
   0.55741977523956205387568692704D+01  0.77559833550341050919015081152D+01
   0.55827668934760399182447978845D+01  0.77799029719487943777002904563D+01
   0.55913098088489172301819281005D+01  0.78037858750736091093948477487D+01
   0.55998266982571082492878504961D+01  0.78276322887766779181325674386D+01
   0.56083177590204408510302977936D+01  0.78514424350375628648888693524D+01
   0.56167831860759478203349224974D+01  0.78752165334834196405787149237D+01
   0.56252231720172532545884353827D+01  0.78989548014244493441613589284D+01
   0.56336379071331340861464695973D+01  0.79226574538886588432145955488D+01
   0.56420275794452776084132558349D+01  0.79463247036559462389402652879D+01
   0.56503923747452552702798844671D+01  0.79699567612915274908308696394D+01
   0.56587324766307324054524720050D+01  0.79935538351787198046937742064D+01
   0.56670480665409329851675855111D+01  0.80171161315510969508347522793D+01
   0.56753393237913779242105412491D+01  0.80406438545240312564096144378D+01
   0.56836064256079149302859134278D+01  0.80641372061256366067464293221D+01
   0.56918495471600573649329862997D+01  0.80875963863271263943273792399D+01
   0.57000688615936490796557373203D+01  0.81110215930725999706246525552D+01
   0.57082645400628717030992279547D+01  0.81344130223082707846541264640D+01
   0.57164367517616103833308856772D+01  0.81577708680111490325073937885D+01
   0.57245856639541935329806481937D+01  0.81810953222171912938278622088D+01
   0.57327114420055216835873836627D+01  0.82043865750489292938078307501D+01
   0.57408142494106001284428769913D+01  0.82276448147425896024142200914D+01
   0.57488942478234896199938098074D+01  0.82508702276747157658298623743D+01
   0.57569515970856889879527351074D+01  0.82740629983883040581684079290D+01
   0.57649864552539631571976249731D+01  0.82972233096184637440414156657D+01
   0.57729989786276296698422057237D+01  0.83203513423176124541968620563D+01
   0.57809893217753164530906672572D+01  0.83434472756802169968927231627D+01
   0.57889576375612032232229218207D+01  0.83665112871670896566132908698D+01
   0.57969040771707585758798914691D+01  0.83895435525292497688870351305D+01
   0.58048287901359844833380976834D+01  0.84125442458313601050418983548D+01
   0.58127319243601796003004517361D+01  0.84355135394747473534664338955D+01
   0.58206136261422324705218317681D+01  0.84584516042200157440730008532D+01
   0.58284740402004554269842606070D+01  0.84813586092092626299319881512D+01
   0.58363133096959696880013794117D+01  0.85042347219879046142228987973D+01
   0.58441315762556518702426079466D+01  0.85270801085261225914972745644D+01
   0.58519289799946518669135465407D+01  0.85498949332399338595467578466D+01
   0.58597056595384917749124223748D+01  0.85726793590118992517022702489D+01
   0.58674617520447552984157827410D+01  0.85954335472114730389505164716D+01
   0.58751973932243768077542304190D+01  0.86181576577150031566426309986D+01
   0.58829127173625389913553313445D+01  0.86408518489253891215949492506D+01
   0.58906078573391878047005223421D+01  0.86635162777914048218588225352D+01
   0.58982829446491731934201769995D+01  0.86861510998266931831870835286D+01
   0.59059381094220238475994702332D+01  0.87087564691284395430776676903D+01
   0.59135734804413640308597028204D+01  0.87313325383957303950646922698D+01
   0.59211891851639803205961934420D+01  0.87538794589476040024946503960D+01
   0.59287853497385458946837571028D+01  0.87763973807407992222166925656D+01
   0.59363620990240098048010234448D+01  0.87988864523872087242831507026D+01
   0.59439195566076584870797650066D+01  0.88213468211710426437567114908D+01
   0.59514578448228565768665582707D+01  0.88437786330657085549162496478D+01
   0.59589770847664739158099487014D+01  0.88661820327504135164114554898D+01
   0.59664773963160054660819209504D+01  0.88885571636264937981088943568D+01
   0.59739588981463906781393339409D+01  0.89109041678334777663753877345D+01
   0.59814217077465386948666211309D+01  0.89332231862648872742393113844D+01
   0.59888659414355656160588957450D+01  0.89555143585837827761414356329D+01
   0.59962917143787498928536867179D+01  0.89777778232380572637231648098D+01
   0.60036991406032117717543191510D+01  0.90000137174754839991941969725D+01
   0.60110883330133225621680939962D+01  0.90222221773585229061701571477D+01
   0.60184594034058493597725574142D+01  0.90444033377788903643736571794D+01
   0.60258124624848407203927188543D+01  0.90665573324718970441529315525D+01
   0.60331476198762586452951256628D+01  0.90886842940305583092974191036D+01
   0.60404649841423621087597085775D+01  0.91107843539194816120293080015D+01
   0.60477646627958472323600134596D+01  0.91328576424885352022370940960D+01
   0.60550467623137490874536615816D+01  0.91549042889863023739075202959D+01
   0.60623113881511099878484029535D+01  0.91769244215733253752245983529D+01
   0.60695586447544190183595007748D+01  0.91989181673351430148602220552D+01
   0.60767886355748274319096090660D+01  0.92208856522951259055042457533D+01
   0.60840014630811444378444830827D+01  0.92428270014271131965994364573D+01
   0.60911972287726177971518681447D+01  0.92647423386678545614874594247D+01
   0.60983760331915035361850664256D+01  0.92866317869292611196674207790D+01
   0.61055379759354289892184250030D+01  0.93084954681104688925521218921D+01
   0.61126831556695532816137696765D+01  0.93303335031097183109149134270D+01
   0.61198116701385292694719687515D+01  0.93521460118360532140898063759D+01
   0.61269236161782708583024760797D+01  0.93739331132208427048592614677D+01
   0.61340190897275295323886826315D+01  0.93956949252291291497797465303D+01
   0.61410981858392838380835937998D+01  0.94174315648708055423985146010D+01
   0.61481609986919454781666269739D+01  0.94391431482116253763517186701D+01
   0.61552076216003855905584693274D+01  0.94608297903840481066512962425D+01
   0.61622381470267847030595374542D+01  0.94824916055979232105150701162D+01
   0.61692526665913097762834520690D+01  0.95041287071510157938218903265D+01
   0.61762512710826216695370407340D+01  0.95257412074393766256336258133D+01
   0.61832340504682162889917347652D+01  0.95473292179675594211721579845D+01
   0.61902010939046026040388508996D+01  0.95688928493586881331274490262D+01
   0.61971524897473206461660826667D+01  0.95904322113643769521588837576D+01
   0.62040883255608025349793614904D+01  0.96119474128745056598944053195D+01
   0.62110086881280795080696601263D+01  0.96334385619268529215897871649D+01
   0.62179136634603378652366003493D+01  0.96549057657165900508442812934D+01
   0.62248033368063266730798534228D+01  0.96763491306056377253406558082D+01
   0.62316777926616200131068490832D+01  0.96977687621318880804502677517D+01
   0.62385371147777364952343459870D+01  0.97191647650182945566814371191D+01
   0.62453813861711186988365651898D+01  0.97405372431818318273172265341D+01
   0.62522106891319751452698893272D+01  0.97618862997423280841530866591D+01
   0.62590251052329873490410143310D+01  0.97832120370311719119730289931D+01
   0.62658247153378844394406795063D+01  0.98045145565998959362633601416D+01
   0.62726095996098877904987621049D+01  0.98257939592286393836248421523D+01
   0.62793798375200280444899206225D+01  0.98470503449344916503776494050D+01
   0.62861355078553368628946303677D+01  0.98682838129797189319297901777D+01
   0.62928766887269156886620643646D+01  0.98894944618798759235707357669D+01
   0.62996034575778837547936470067D+01  0.99106823894118045624306775432D+01
   0.63063158911912075266351489882D+01  0.99318476926215217403857505764D+01
   0.63130140656974137187978521755D+01  0.99529904678319978786651450772D+01
   0.63196980565821879822935592124D+01  0.99741108106508282168024572449D+01
   0.63263679386938613132330022629D+01  0.99952088159777986313468273920D+01
   0.63330237862507861912724132877D+01  0.10016284578012347763386004504D+02
   0.63396656728486044138694652132D+01  0.10037338190260927198410779414D+02
   0.63462936714674085512991748844D+01  0.10058369745544261407346225609D+02
c
c
                                  file qzwlag
                     (containing the recursion coefficients
          for the generalized Laguerre polynomials with parameter -1/6)
c
c
   0.137546237035606481016626358259D-01  0.797263205789783349012395912137D-01
   0.828242819612625076589717747103D-01  0.137509571840463036234285241382D+00
   0.211045266878820067433757899195D+00  0.165624629277404627483798099224D+00
   0.398495883697246106215870196186D+00  0.169889385184249933348754420087D+00
   0.645246449413364192790056305710D+00  0.155979710986940400443018376349D+00
   0.951385724271680753484465397335D+00  0.130796100640281569497382999057D+00
   0.131702333539251889691557844124D+01  0.101195207699248819155416108185D+00
   0.174229039557434737528689876878D+01  0.726578417587654081289565811728D-01
   0.222733986799705719106134454584D+01  0.485868826476249722732414580795D-01
   0.277234690600003415615924289811D+01  0.303307104631899399813923717524D-01
   0.337750921867510044213142555742D+01  0.177032794850510052563288407594D-01
   0.404304747760954736833747106973D+01  0.967156114455361607595781378153D-02
   0.476920577124043653443101103548D+01  0.494910619781186603083484746505D-02
   0.555625211076505238359064027062D+01  0.237330404386327067480928921695D-02
   0.640447899087500233247019669842D+01  0.106685632091733307729458416886D-02
   0.731420400852675468739052130926D+01  0.449627009561703289099775291088D-03
   0.828577054315434627885977150222D+01  0.177668897044701158278497052572D-03
   0.931954850204318353083209913858D+01  0.658204611789929317057160248270D-04
   0.104159351349719001116590618545D+02  0.228583177116227029870073572656D-04
   0.115753559226768156120308674906D+02  0.744002738337168885093518101880D-05
   0.127982655441976107509841188907D+02  0.226901703190553122504707139869D-05
   0.140851489287260062254247699931D+02  0.648179317077785444771156745007D-06
   0.154365223982072337029903782879D+02  0.173374350391629422565860207969D-06
   0.168529349076447943295975436053D+02  0.434033071931337973251399342238D-07
   0.183349693908421258472298414543D+02  0.101649166902660479947132813417D-07
   0.198832442201918964152204611595D+02  0.222587973662653269368927846641D-08
   0.214984147901082737711203462908D+02  0.455480274668852349666205928872D-09
   0.231811752348055400413574663882D+02  0.870441993368912701957937182930D-10
   0.249322602923763714853936899471D+02  0.155247232505992515917844817607D-10
   0.267524473285366763289675037486D+02  0.258232958449522171238153558175D-11
   0.286425585350065368627866899483D+02  0.400286834186069782879872701788D-12
   0.306034633193184737637221713345D+02  0.577760961243634423791907788250D-13
   0.326360809049202573987031376541D+02  0.775825902274054523741659235002D-14
   0.347413831628116175240772736421D+02  0.968315825957348296265047292700D-15
   0.369203976986720441451715025151D+02  0.112222088751142422381790694889D-15
   0.391742112225594869113645755498D+02  0.120640534139625623786906828214D-16
   0.415039732318577168732425329215D+02  0.120165129567887359460562589892D-17
   0.439109000423080337295588679129D+02  0.110769648715769618173706021069D-18
   0.463962792067807294830331018748D+02  0.943789878500328636365175630221D-20
   0.489614743670463016078379152673D+02  0.742272692485229651999270764612D-21
   0.516079305903450972575107611951D+02  0.538109662993620555668774988116D-22
   0.543371802502086556629014574405D+02  0.359041280578153298065891265661D-23
   0.571508495199788819177592023398D+02  0.220135619307550878039271691198D-24
   0.600506655580756017373182580786D+02  0.123814220804856862367003999666D-25
   0.630384644766164020763062685364D+02  0.637675072786613776291097782778D-27
   0.661162001999135699395076038284D+02  0.300151290553555160313395799834D-28
   0.692859543371834871289812608312D+02  0.128854781356925667816580899399D-29
   0.725499472151586993304909197158D+02  0.503417039940762192462787883299D-31
   0.759105502420176879726828903691D+02  0.178568301155242434486452880392D-32
   0.793702998051896251298950709287D+02  0.573641442444169903709503390520D-34
   0.829319129434884668207156165314D+02  0.166444272509331931898542571984D-35
   0.865983050804015380146409935328D+02  0.434946794230483177409677542063D-37
   0.903726101624271247020764292365D+02  0.102045710427243504101765021931D-38
   0.942582036170247130721839820192D+02  0.214236223408145003516338376722D-40
   0.982587286328235051665608413176D+02  0.401017474137338970771521188961D-42
   0.102378126375283428945390999867D+03  0.666668913957876994400407470676D-44
   0.106620670890775924372995337396D+03  0.980154382381083332882641731528D-46
   0.110991009630172067651738754487D+03  0.126857516554975830826433092719D-47
   0.115494210751933587302323855079D+03  0.143811741447024234565344597829D-49
   0.120135818661547724937258921499D+03  0.142017869459526959683205241018D-51
   0.124921919632847734813611211887D+03  0.121434812626150933068394008129D-53
   0.129859219871094635968162995239D+03  0.893119478583897537667093867055D-56
   0.134955139066012235999553635266D+03  0.560854313044256699167304859348D-58
   0.140217923415493268362340097732D+03  0.298275616129906752106302388632D-60
   0.145656783381064149817949031613D+03  0.133122576316827967154259115778D-62
   0.151282063220011736775471598572D+03  0.493519765761388142050959783909D-65
   0.157105451864431884047737138328D+03  0.150228448116207251553416605374D-67
   0.163140248356336676851891466858D+03  0.370578636149056518630857554463D-70
   0.169401700395024578833681757223D+03  0.729692956908954916674429601704D-73
   0.175907442583447866375165292583D+03  0.112709538396289297639581375623D-75
   0.182678073320112747826259068435D+03  0.133814497292966958468592484789D-78
   0.189737928840241865645706653138D+03  0.119215045733700130665897289239D-81
   0.197116144848558851613910232716D+03  0.774386887778651369750370248306D-85
   0.204848150276338368407383895840D+03  0.354170552600900452690923626171D-88
   0.212977833265267549890905711157D+03  0.109223956430185599161239404714D-91
   0.221560797004253251076271759414D+03  0.215036403807807412441528603035D-95
   0.230669473231034355966830510998D+03  0.251699205160681803175493652384D-99
   0.240401604862943635610583231380D+03  0.159140640969695084408123457406-103
   0.250895342301616596497086665149D+03  0.474756323951450620940889550645-108
   0.262358759120349589116996153078D+03  0.545418947814451885889605946511-113
   0.275135789300494512317254606969D+03  0.172709205033590742874063800203-118
   0.289887449630615829882626700778D+03  0.795821688082107118597523913369-125
   0.308324624796163197455423024909D+03  0.104189157452535659404726837861-132
