function [ xtab, weight ] = laguerre_set ( norder )

%% LAGUERRE_SET sets abscissas and weights for Laguerre quadrature.
%
%  Discussion:
%
%    The integration interval is [ 0, +Infinity ).
%
%    The weight function w(x) = EXP ( - X );
%
%    The abscissas are the zeroes of the Laguerre polynomial L(NORDER)(X).
%
%
%    If the integral to approximate is:
%
%      Integral ( 0 <= X < +Infinity ) exp ( -X ) * F(X) dX
%
%    then the quadrature rule is:
%
%      Sum ( 1 <= I <= ORDER ) WEIGHT(I) * f ( XTAB(I) )
%
%    If the integral to approximate is:
%
%      Integral ( 0 <= X < +Infinity ) F(X) dX
%
%    then the quadrature rule is:
%
%      Sum ( 1 <= I <= ORDER ) WEIGHT(I) * exp ( XTAB(I) ) * f ( XTAB(I) )
%
%  Licensing:
%
%    This code is distributed under the GNU LGPL license.
%
%  Modified:
%
%    05 October 2007
%
%  Author:
%
%    John Burkardt
%
%  Reference:
%
%    Milton Abramowitz, Irene Stegun,
%    Handbook of Mathematical Functions,
%    National Bureau of Standards, 1964,
%    ISBN: 0-486-61272-4,
%    LC: QA47.A34.
%
%    Vladimir Krylov,
%    Approximate Calculation of Integrals,
%    Dover, 2006,
%    ISBN: 0486445798,
%    LC: QA311.K713.
%
%    Arthur Stroud, Don Secrest,
%    Gaussian Quadrature Formulas,
%    Prentice Hall, 1966,
%    LC: QA299.4G3S7.
%
%    Daniel Zwillinger, editor,
%    CRC Standard Mathematical Tables and Formulae,
%    30th Edition,
%    CRC Press, 1996,
%    ISBN: 0-8493-2479-3,
%    LC: QA47.M315.
%
%  Parameters:
%
%    Input, integer NORDER, the order of the rule.
%    NORDER must be between 1 and 20, 31, 32, 63, 64 or 127.
%
%    Output, real XTAB(NORDER), the abscissas of the rule.
%
%    Output, real WEIGHT(NORDER), the weights of the rule.
%    The weights are positive, and should add to 1.
%
  if ( norder == 1 )

    xtab(1) =    1.0;
    weight(1) =  1.0;

  elseif ( norder == 2 )

    xtab(1) =    0.585786437626904951198311275790;
    xtab(2) =    0.341421356237309504880168872421E+01;

    weight(1) =  0.853553390593273762200422181052;
    weight(2) =  0.146446609406726237799577818948;

  elseif ( norder == 3 )

    xtab(1) =    0.415774556783479083311533873128;
    xtab(2) =    0.229428036027904171982205036136E+01;
    xtab(3) =    0.628994508293747919686641576551E+01;

    weight(1) =  0.711093009929173015449590191143;
    weight(2) =  0.278517733569240848801444888457;
    weight(3) =  0.103892565015861357489649204007E-01;

  elseif ( norder == 4 )

    xtab(1) =    0.322547689619392311800361943361;
    xtab(2) =    0.174576110115834657568681671252E+01;
    xtab(3) =    0.453662029692112798327928538496E+01;
    xtab(4) =    0.939507091230113312923353644342E+01;

    weight(1) =  0.603154104341633601635966023818;
    weight(2) =  0.357418692437799686641492017458;
    weight(3) =  0.388879085150053842724381681562E-01;
    weight(4) =  0.539294705561327450103790567621E-03;

  elseif ( norder == 5 )

    xtab(1) =    0.263560319718140910203061943361;
    xtab(2) =    0.141340305910651679221840798019E+01;
    xtab(3) =    0.359642577104072208122318658878E+01;
    xtab(4) =    0.708581000585883755692212418111E+01;
    xtab(5) =    0.126408008442757826594332193066E+02;

    weight(1) =  0.521755610582808652475860928792;
    weight(2) =  0.398666811083175927454133348144;
    weight(3) =  0.759424496817075953876533114055E-01;
    weight(4) =  0.361175867992204845446126257304E-02;
    weight(5) =  0.233699723857762278911490845516E-04;

  elseif ( norder == 6 )

    xtab(1) =    0.222846604179260689464354826787;
    xtab(2) =    0.118893210167262303074315092194E+01;
    xtab(3) =    0.299273632605931407769132528451E+01;
    xtab(4) =    0.577514356910451050183983036943E+01;
    xtab(5) =    0.983746741838258991771554702994E+01;
    xtab(6) =    0.159828739806017017825457915674E+02;

    weight(1) =  0.458964673949963593568284877709;
    weight(2) =  0.417000830772120994113377566193;
    weight(3) =  0.113373382074044975738706185098;
    weight(4) =  0.103991974531490748989133028469E-01;
    weight(5) =  0.261017202814932059479242860001E-03;
    weight(6) =  0.898547906429621238825292052825E-06;

  elseif ( norder == 7 )

    xtab(1) =    0.193043676560362413838247885004;
    xtab(2) =    0.102666489533919195034519944317E+01;
    xtab(3) =    0.256787674495074620690778622666E+01;
    xtab(4) =    0.490035308452648456810171437810E+01;
    xtab(5) =    0.818215344456286079108182755123E+01;
    xtab(6) =    0.127341802917978137580126424582E+02;
    xtab(7) =    0.193957278622625403117125820576E+02;

    weight(1) =  0.409318951701273902130432880018;
    weight(2) =  0.421831277861719779929281005417;
    weight(3) =  0.147126348657505278395374184637;
    weight(4) =  0.206335144687169398657056149642E-01;
    weight(5) =  0.107401014328074552213195962843E-02;
    weight(6) =  0.158654643485642012687326223234E-04;
    weight(7) =  0.317031547899558056227132215385E-07;

  elseif ( norder == 8 )

    xtab(1) =    0.170279632305100999788861856608;
    xtab(2) =    0.903701776799379912186020223555;
    xtab(3) =    0.225108662986613068930711836697E+01;
    xtab(4) =    0.426670017028765879364942182690E+01;
    xtab(5) =    0.704590540239346569727932548212E+01;
    xtab(6) =    0.107585160101809952240599567880E+02;
    xtab(7) =    0.157406786412780045780287611584E+02;
    xtab(8) =    0.228631317368892641057005342974E+02;

    weight(1) =  0.369188589341637529920582839376;
    weight(2) =  0.418786780814342956076978581333;
    weight(3) =  0.175794986637171805699659866777;
    weight(4) =  0.333434922612156515221325349344E-01;
    weight(5) =  0.279453623522567252493892414793E-02;
    weight(6) =  0.907650877335821310423850149336E-04;
    weight(7) =  0.848574671627253154486801830893E-06;
    weight(8) =  0.104800117487151038161508853552E-08;

  elseif ( norder == 9 )

    xtab(1) =    0.152322227731808247428107073127;
    xtab(2) =    0.807220022742255847741419210952;
    xtab(3) =    0.200513515561934712298303324701E+01;
    xtab(4) =    0.378347397333123299167540609364E+01;
    xtab(5) =    0.620495677787661260697353521006E+01;
    xtab(6) =    0.937298525168757620180971073215E+01;
    xtab(7) =    0.134662369110920935710978818397E+02;
    xtab(8) =    0.188335977889916966141498992996E+02;
    xtab(9) =    0.263740718909273767961410072937E+02;

    weight(1) =  0.336126421797962519673467717606;
    weight(2) =  0.411213980423984387309146942793;
    weight(3) =  0.199287525370885580860575607212;
    weight(4) =  0.474605627656515992621163600479E-01;
    weight(5) =  0.559962661079458317700419900556E-02;
    weight(6) =  0.305249767093210566305412824291E-03;
    weight(7) =  0.659212302607535239225572284875E-05;
    weight(8) =  0.411076933034954844290241040330E-07;
    weight(9) =  0.329087403035070757646681380323E-10;

  elseif ( norder == 10 )

    xtab(1) =    0.137793470540492430830772505653;
    xtab(2) =    0.729454549503170498160373121676;
    xtab(3) =    0.180834290174031604823292007575E+01;
    xtab(4) =    0.340143369785489951448253222141E+01;
    xtab(5) =    0.555249614006380363241755848687E+01;
    xtab(6) =    0.833015274676449670023876719727E+01;
    xtab(7) =    0.118437858379000655649185389191E+02;
    xtab(8) =    0.162792578313781020995326539358E+02;
    xtab(9) =    0.219965858119807619512770901956E+02;
    xtab(10) =   0.299206970122738915599087933408E+02;

    weight(1) =  0.308441115765020141547470834678;
    weight(2) =  0.401119929155273551515780309913;
    weight(3) =  0.218068287611809421588648523475;
    weight(4) =  0.620874560986777473929021293135E-01;
    weight(5) =  0.950151697518110055383907219417E-02;
    weight(6) =  0.753008388587538775455964353676E-03;
    weight(7) =  0.282592334959956556742256382685E-04;
    weight(8) =  0.424931398496268637258657665975E-06;
    weight(9) =  0.183956482397963078092153522436E-08;
    weight(10) = 0.991182721960900855837754728324E-12;

  elseif ( norder == 11 )

    xtab(1) =    0.125796442187967522675794577516;
    xtab(2) =    0.665418255839227841678127839420;
    xtab(3) =    0.164715054587216930958700321365E+01;
    xtab(4) =    0.309113814303525495330195934259E+01;
    xtab(5) =    0.502928440157983321236999508366E+01;
    xtab(6) =    0.750988786380661681941099714450E+01;
    xtab(7) =    0.106059509995469677805559216457E+02;
    xtab(8) =    0.144316137580641855353200450349E+02;
    xtab(9) =    0.191788574032146786478174853989E+02;
    xtab(10) =   0.252177093396775611040909447797E+02;
    xtab(11) =   0.334971928471755372731917259395E+02;

    weight(1) =  0.284933212894200605056051024724;
    weight(2) =  0.389720889527849377937553508048;
    weight(3) =  0.232781831848991333940223795543;
    weight(4) =  0.765644535461966864008541790132E-01;
    weight(5) =  0.143932827673506950918639187409E-01;
    weight(6) =  0.151888084648487306984777640042E-02;
    weight(7) =  0.851312243547192259720424170600E-04;
    weight(8) =  0.229240387957450407857683270709E-05;
    weight(9) =  0.248635370276779587373391491114E-07;
    weight(10) = 0.771262693369132047028152590222E-10;
    weight(11) = 0.288377586832362386159777761217E-13;

  elseif ( norder == 12 )

    xtab(1) =    0.115722117358020675267196428240;
    xtab(2) =    0.611757484515130665391630053042;
    xtab(3) =    0.151261026977641878678173792687E+01;
    xtab(4) =    0.283375133774350722862747177657E+01;
    xtab(5) =    0.459922763941834848460572922485E+01;
    xtab(6) =    0.684452545311517734775433041849E+01;
    xtab(7) =    0.962131684245686704391238234923E+01;
    xtab(8) =    0.130060549933063477203460524294E+02;
    xtab(9) =    0.171168551874622557281840528008E+02;
    xtab(10) =   0.221510903793970056699218950837E+02;
    xtab(11) =   0.284879672509840003125686072325E+02;
    xtab(12) =   0.370991210444669203366389142764E+02;

    weight(1) =  0.264731371055443190349738892056;
    weight(2) =  0.377759275873137982024490556707;
    weight(3) =  0.244082011319877564254870818274;
    weight(4) =  0.904492222116809307275054934667E-01;
    weight(5) =  0.201023811546340965226612867827E-01;
    weight(6) =  0.266397354186531588105415760678E-02;
    weight(7) =  0.203231592662999392121432860438E-03;
    weight(8) =  0.836505585681979874533632766396E-05;
    weight(9) =  0.166849387654091026116989532619E-06;
    weight(10) = 0.134239103051500414552392025055E-08;
    weight(11) = 0.306160163503502078142407718971E-11;
    weight(12) = 0.814807746742624168247311868103E-15;

  elseif ( norder == 13 )

    xtab(1) =    0.107142388472252310648493376977;
    xtab(2) =    0.566131899040401853406036347177;
    xtab(3) =    0.139856433645101971792750259921E+01;
    xtab(4) =    0.261659710840641129808364008472E+01;
    xtab(5) =    0.423884592901703327937303389926E+01;
    xtab(6) =    0.629225627114007378039376523025E+01;
    xtab(7) =    0.881500194118697804733348868036E+01;
    xtab(8) =    0.118614035888112425762212021880E+02;
    xtab(9) =    0.155107620377037527818478532958E+02;
    xtab(10) =   0.198846356638802283332036594634E+02;
    xtab(11) =   0.251852638646777580842970297823E+02;
    xtab(12) =   0.318003863019472683713663283526E+02;
    xtab(13) =   0.407230086692655795658979667001E+02;

    weight(1) =  0.247188708429962621346249185964;
    weight(2) =  0.365688822900521945306717530893;
    weight(3) =  0.252562420057658502356824288815;
    weight(4) =  0.103470758024183705114218631672;
    weight(5) =  0.264327544155616157781587735702E-01;
    weight(6) =  0.422039604025475276555209292644E-02;
    weight(7) =  0.411881770472734774892472527082E-03;
    weight(8) =  0.235154739815532386882897300772E-04;
    weight(9) =  0.731731162024909910401047197761E-06;
    weight(10) = 0.110884162570398067979150974759E-07;
    weight(11) = 0.677082669220589884064621459082E-10;
    weight(12) = 0.115997995990507606094507145382E-12;
    weight(13) = 0.224509320389275841599187226865E-16;

  elseif ( norder == 14 )

    xtab(1) =    0.997475070325975745736829452514E-01;
    xtab(2) =    0.526857648851902896404583451502;
    xtab(3) =    0.130062912125149648170842022116E+01;
    xtab(4) =    0.243080107873084463616999751038E+01;
    xtab(5) =    0.393210282229321888213134366778E+01;
    xtab(6) =    0.582553621830170841933899983898E+01;
    xtab(7) =    0.814024014156514503005978046052E+01;
    xtab(8) =    0.109164995073660188408130510904E+02;
    xtab(9) =    0.142108050111612886831059780825E+02;
    xtab(10) =   0.181048922202180984125546272083E+02;
    xtab(11) =   0.227233816282696248232280886985E+02;
    xtab(12) =   0.282729817232482056954158923218E+02;
    xtab(13) =   0.351494436605924265828643121364E+02;
    xtab(14) =   0.443660817111174230416312423666E+02;

    weight(1) =  0.231815577144864977840774861104;
    weight(2) =  0.353784691597543151802331301273;
    weight(3) =  0.258734610245428085987320561144;
    weight(4) =  0.115482893556923210087304988673;
    weight(5) =  0.331920921593373600387499587137E-01;
    weight(6) =  0.619286943700661021678785967675E-02;
    weight(7) =  0.739890377867385942425890907080E-03;
    weight(8) =  0.549071946684169837857331777667E-04;
    weight(9) =  0.240958576408537749675775256553E-05;
    weight(10) = 0.580154398167649518088619303904E-07;
    weight(11) = 0.681931469248497411961562387084E-09;
    weight(12) = 0.322120775189484793980885399656E-11;
    weight(13) = 0.422135244051658735159797335643E-14;
    weight(14) = 0.605237502228918880839870806281E-18;

  elseif ( norder == 15 )

    xtab(1) =    0.933078120172818047629030383672E-01;
    xtab(2) =    0.492691740301883908960101791412;
    xtab(3) =    0.121559541207094946372992716488E+01;
    xtab(4) =    0.226994952620374320247421741375E+01;
    xtab(5) =    0.366762272175143727724905959436E+01;
    xtab(6) =    0.542533662741355316534358132596E+01;
    xtab(7) =    0.756591622661306786049739555812E+01;
    xtab(8) =    0.101202285680191127347927394568E+02;
    xtab(9) =    0.131302824821757235640991204176E+02;
    xtab(10) =   0.166544077083299578225202408430E+02;
    xtab(11) =   0.207764788994487667729157175676E+02;
    xtab(12) =   0.256238942267287801445868285977E+02;
    xtab(13) =   0.314075191697539385152432196202E+02;
    xtab(14) =   0.385306833064860094162515167595E+02;
    xtab(15) =   0.480260855726857943465734308508E+02;

    weight(1) =  0.218234885940086889856413236448;
    weight(2) =  0.342210177922883329638948956807;
    weight(3) =  0.263027577941680097414812275022;
    weight(4) =  0.126425818105930535843030549378;
    weight(5) =  0.402068649210009148415854789871E-01;
    weight(6) =  0.856387780361183836391575987649E-02;
    weight(7) =  0.121243614721425207621920522467E-02;
    weight(8) =  0.111674392344251941992578595518E-03;
    weight(9) =  0.645992676202290092465319025312E-05;
    weight(10) = 0.222631690709627263033182809179E-06;
    weight(11) = 0.422743038497936500735127949331E-08;
    weight(12) = 0.392189726704108929038460981949E-10;
    weight(13) = 0.145651526407312640633273963455E-12;
    weight(14) = 0.148302705111330133546164737187E-15;
    weight(15) = 0.160059490621113323104997812370E-19;

  elseif ( norder == 16 )

    xtab(1) =    0.876494104789278403601980973401E-01;
    xtab(2) =    0.462696328915080831880838260664;
    xtab(3) =    0.114105777483122685687794501811E+01;
    xtab(4) =    0.212928364509838061632615907066E+01;
    xtab(5) =    0.343708663389320664523510701675E+01;
    xtab(6) =    0.507801861454976791292305830814E+01;
    xtab(7) =    0.707033853504823413039598947080E+01;
    xtab(8) =    0.943831433639193878394724672911E+01;
    xtab(9) =    0.122142233688661587369391246088E+02;
    xtab(10) =   0.154415273687816170767647741622E+02;
    xtab(11) =   0.191801568567531348546631409497E+02;
    xtab(12) =   0.235159056939919085318231872752E+02;
    xtab(13) =   0.285787297428821403675206137099E+02;
    xtab(14) =   0.345833987022866258145276871778E+02;
    xtab(15) =   0.419404526476883326354722330252E+02;
    xtab(16) =   0.517011603395433183643426971197E+02;

    weight(1) =  0.206151714957800994334273636741;
    weight(2) =  0.331057854950884165992983098710;
    weight(3) =  0.265795777644214152599502020650;
    weight(4) =  0.136296934296377539975547513526;
    weight(5) =  0.473289286941252189780623392781E-01;
    weight(6) =  0.112999000803394532312490459701E-01;
    weight(7) =  0.184907094352631086429176783252E-02;
    weight(8) =  0.204271915308278460126018338421E-03;
    weight(9) =  0.148445868739812987713515067551E-04;
    weight(10) = 0.682831933087119956439559590327E-06;
    weight(11) = 0.188102484107967321388159920418E-07;
    weight(12) = 0.286235024297388161963062629156E-09;
    weight(13) = 0.212707903322410296739033610978E-11;
    weight(14) = 0.629796700251786778717446214552E-14;
    weight(15) = 0.505047370003551282040213233303E-17;
    weight(16) = 0.416146237037285519042648356116E-21;

  elseif ( norder == 17 )

    xtab(1) =    0.826382147089476690543986151980E-01;
    xtab(2) =    0.436150323558710436375959029847;
    xtab(3) =    0.107517657751142857732980316755E+01;
    xtab(4) =    0.200519353164923224070293371933E+01;
    xtab(5) =    0.323425612404744376157380120696E+01;
    xtab(6) =    0.477351351370019726480932076262E+01;
    xtab(7) =    0.663782920536495266541643929703E+01;
    xtab(8) =    0.884668551116980005369470571184E+01;
    xtab(9) =    0.114255293193733525869726151469E+02;
    xtab(10) =   0.144078230374813180021982874959E+02;
    xtab(11) =   0.178382847307011409290658752412E+02;
    xtab(12) =   0.217782682577222653261749080522E+02;
    xtab(13) =   0.263153178112487997766149598369E+02;
    xtab(14) =   0.315817716804567331343908517497E+02;
    xtab(15) =   0.377960938374771007286092846663E+02;
    xtab(16) =   0.453757165339889661829258363215E+02;
    xtab(17) =   0.553897517898396106640900199790E+02;

    weight(1) =  0.195332205251770832145927297697;
    weight(2) =  0.320375357274540281336625631970;
    weight(3) =  0.267329726357171097238809604160;
    weight(4) =  0.145129854358758625407426447473;
    weight(5) =  0.544369432453384577793805803066E-01;
    weight(6) =  0.143572977660618672917767247431E-01;
    weight(7) =  0.266282473557277256843236250006E-02;
    weight(8) =  0.343679727156299920611775097985E-03;
    weight(9) =  0.302755178378287010943703641131E-04;
    weight(10) = 0.176851505323167689538081156159E-05;
    weight(11) = 0.657627288681043332199222748162E-07;
    weight(12) = 0.146973093215954679034375821888E-08;
    weight(13) = 0.181691036255544979555476861323E-10;
    weight(14) = 0.109540138892868740297645078918E-12;
    weight(15) = 0.261737388222337042155132062413E-15;
    weight(16) = 0.167293569314615469085022374652E-18;
    weight(17) = 0.106562631627404278815253271162E-22;

  elseif ( norder == 18 )

    xtab(1) =    0.781691666697054712986747615334E-01;
    xtab(2) =    0.412490085259129291039101536536;
    xtab(3) =    0.101652017962353968919093686187E+01;
    xtab(4) =    0.189488850996976091426727831954E+01;
    xtab(5) =    0.305435311320265975115241130719E+01;
    xtab(6) =    0.450420553888989282633795571455E+01;
    xtab(7) =    0.625672507394911145274209116326E+01;
    xtab(8) =    0.832782515660563002170470261564E+01;
    xtab(9) =    0.107379900477576093352179033397E+02;
    xtab(10) =   0.135136562075550898190863812108E+02;
    xtab(11) =   0.166893062819301059378183984163E+02;
    xtab(12) =   0.203107676262677428561313764553E+02;
    xtab(13) =   0.244406813592837027656442257980E+02;
    xtab(14) =   0.291682086625796161312980677805E+02;
    xtab(15) =   0.346279270656601721454012429438E+02;
    xtab(16) =   0.410418167728087581392948614284E+02;
    xtab(17) =   0.488339227160865227486586093290E+02;
    xtab(18) =   0.590905464359012507037157810181E+02;

    weight(1) =  0.185588603146918805623337752284;
    weight(2) =  0.310181766370225293649597595713;
    weight(3) =  0.267866567148536354820854394783;
    weight(4) =  0.152979747468074906553843082053;
    weight(5) =  0.614349178609616527076780103487E-01;
    weight(6) =  0.176872130807729312772600233761E-01;
    weight(7) =  0.366017976775991779802657207890E-02;
    weight(8) =  0.540622787007735323128416319257E-03;
    weight(9) =  0.561696505121423113817929049294E-04;
    weight(10) = 0.401530788370115755858883625279E-05;
    weight(11) = 0.191466985667567497969210011321E-06;
    weight(12) = 0.583609526863159412918086289717E-08;
    weight(13) = 0.107171126695539012772851317562E-09;
    weight(14) = 0.108909871388883385562011298291E-11;
    weight(15) = 0.538666474837830887608094323164E-14;
    weight(16) = 0.104986597803570340877859934846E-16;
    weight(17) = 0.540539845163105364356554467358E-20;
    weight(18) = 0.269165326920102862708377715980E-24;

  elseif ( norder == 19 )

    xtab(1) =    0.741587837572050877131369916024E-01;
    xtab(2) =    0.391268613319994607337648350299;
    xtab(3) =    0.963957343997958058624879377130;
    xtab(4) =    0.179617558206832812557725825252E+01;
    xtab(5) =    0.289365138187378399116494713237E+01;
    xtab(6) =    0.426421553962776647436040018167E+01;
    xtab(7) =    0.591814156164404855815360191408E+01;
    xtab(8) =    0.786861891533473373105668358176E+01;
    xtab(9) =    0.101324237168152659251627415800E+02;
    xtab(10) =   0.127308814638423980045092979656E+02;
    xtab(11) =   0.156912783398358885454136069861E+02;
    xtab(12) =   0.190489932098235501532136429732E+02;
    xtab(13) =   0.228508497608294829323930586693E+02;
    xtab(14) =   0.271606693274114488789963947149E+02;
    xtab(15) =   0.320691222518622423224362865906E+02;
    xtab(16) =   0.377129058012196494770647508283E+02;
    xtab(17) =   0.443173627958314961196067736013E+02;
    xtab(18) =   0.523129024574043831658644222420E+02;
    xtab(19) =   0.628024231535003758413504690673E+02;

    weight(1) =  0.176768474915912502251035479815;
    weight(2) =  0.300478143607254379482156807712;
    weight(3) =  0.267599547038175030772695440648;
    weight(4) =  0.159913372135580216785512147895;
    weight(5) =  0.682493799761491134552355368344E-01;
    weight(6) =  0.212393076065443249244062193091E-01;
    weight(7) =  0.484162735114839596725013121019E-02;
    weight(8) =  0.804912747381366766594647138204E-03;
    weight(9) =  0.965247209315350170843161738801E-04;
    weight(10) = 0.820730525805103054408982992869E-05;
    weight(11) = 0.483056672473077253944806671560E-06;
    weight(12) = 0.190499136112328569993615674552E-07;
    weight(13) = 0.481668463092806155766936380273E-09;
    weight(14) = 0.734825883955114437684376840171E-11;
    weight(15) = 0.620227538757261639893719012423E-13;
    weight(16) = 0.254143084301542272371866857954E-15;
    weight(17) = 0.407886129682571235007187465134E-18;
    weight(18) = 0.170775018759383706100412325084E-21;
    weight(19) = 0.671506464990818995998969111749E-26;

  elseif ( norder == 20 )

    xtab(1) =    0.705398896919887533666890045842E-01;
    xtab(2) =    0.372126818001611443794241388761;
    xtab(3) =    0.916582102483273564667716277074;
    xtab(4) =    0.170730653102834388068768966741E+01;
    xtab(5) =    0.274919925530943212964503046049E+01;
    xtab(6) =    0.404892531385088692237495336913E+01;
    xtab(7) =    0.561517497086161651410453988565E+01;
    xtab(8) =    0.745901745367106330976886021837E+01;
    xtab(9) =    0.959439286958109677247367273428E+01;
    xtab(10) =   0.120388025469643163096234092989E+02;
    xtab(11) =   0.148142934426307399785126797100E+02;
    xtab(12) =   0.179488955205193760173657909926E+02;
    xtab(13) =   0.214787882402850109757351703696E+02;
    xtab(14) =   0.254517027931869055035186774846E+02;
    xtab(15) =   0.299325546317006120067136561352E+02;
    xtab(16) =   0.350134342404790000062849359067E+02;
    xtab(17) =   0.408330570567285710620295677078E+02;
    xtab(18) =   0.476199940473465021399416271529E+02;
    xtab(19) =   0.558107957500638988907507734445E+02;
    xtab(20) =   0.665244165256157538186403187915E+02;

    weight(1) =  0.168746801851113862149223899689;
    weight(2) =  0.291254362006068281716795323812;
    weight(3) =  0.266686102867001288549520868998;
    weight(4) =  0.166002453269506840031469127816;
    weight(5) =  0.748260646687923705400624639615E-01;
    weight(6) =  0.249644173092832210728227383234E-01;
    weight(7) =  0.620255084457223684744754785395E-02;
    weight(8) =  0.114496238647690824203955356969E-02;
    weight(9) =  0.155741773027811974779809513214E-03;
    weight(10) = 0.154014408652249156893806714048E-04;
    weight(11) = 0.108648636651798235147970004439E-05;
    weight(12) = 0.533012090955671475092780244305E-07;
    weight(13) = 0.175798117905058200357787637840E-08;
    weight(14) = 0.372550240251232087262924585338E-10;
    weight(15) = 0.476752925157819052449488071613E-12;
    weight(16) = 0.337284424336243841236506064991E-14;
    weight(17) = 0.115501433950039883096396247181E-16;
    weight(18) = 0.153952214058234355346383319667E-19;
    weight(19) = 0.528644272556915782880273587683E-23;
    weight(20) = 0.165645661249902329590781908529E-27;

  elseif ( order == 31 )

    xtab(  1) =   0.45901947621108290743496080275224E-01;
    xtab(  2) =   0.24198016382477204890408974151714E+00;
    xtab(  3) =   0.59525389422235073707330165005414E+00;
    xtab(  4) =    1.1066894995329987162111308789792E+00;
    xtab(  5) =    1.7775956928747727211593727482675E+00;
    xtab(  6) =    2.6097034152566806503893375925315E+00;
    xtab(  7) =    3.6051968023400442698805817554243E+00;
    xtab(  8) =    4.7667470844717611313629127271123E+00;
    xtab(  9) =    6.0975545671817409269925429328463E+00;
    xtab( 10) =    7.6014009492331374229360106942867E+00;
    xtab( 11) =    9.2827143134708894182536695297710E+00;
    xtab( 12) =    11.146649755619291358993815629587E+00;
    xtab( 13) =    13.199189576244998522464925028637E+00;
    xtab( 14) =    15.447268315549310075809325891801E+00;
    xtab( 15) =    17.898929826644757646725793817752E+00;
    xtab( 16) =    20.563526336715822170743048968779E+00;
    xtab( 17) =    23.451973482011858591050255575933E+00;
    xtab( 18) =    26.577081352118260459975876986478E+00;
    xtab( 19) =    29.953990872346445506951917840024E+00;
    xtab( 20) =    33.600759532902202735410313885784E+00;
    xtab( 21) =    37.539164407330440882887902558001E+00;
    xtab( 22) =    41.795830870182219981347945853330E+00;
    xtab( 23) =    46.403866806411123136029227604386E+00;
    xtab( 24) =    51.405314476797755161861461088395E+00;
    xtab( 25) =    56.854992868715843620511922055660E+00;
    xtab( 26) =    62.826855908786321453677523304806E+00;
    xtab( 27) =    69.425277191080345623322251656443E+00;
    xtab( 28) =    76.807047763862732837609972285484E+00;
    xtab( 29) =    85.230358607545669169387065607043E+00;
    xtab( 30) =    95.188939891525629981308606853957E+00;
    xtab( 31) =    107.95224382757871475002440117666E+00;
 
    weight(  1) =   0.11252789550372583820847728082801E+00;
    weight(  2) =   0.21552760818089123795222505285045E+00;
    weight(  3) =   0.23830825164569654731905788089234E+00;
    weight(  4) =   0.19538830929790229249915303390711E+00;
    weight(  5) =   0.12698283289306190143635272904602E+00;
    weight(  6) =   0.67186168923899300670929441993508E-01;
    weight(  7) =   0.29303224993879487404888669311974E-01;
    weight(  8) =   0.10597569915295736089529380314433E-01;
    weight(  9) =   0.31851272582386980320974842433019E-02;
    weight( 10) =   0.79549548307940382922092149012477E-03;
    weight( 11) =   0.16480052126636687317862967116412E-03;
    weight( 12) =   0.28229237864310816393860971468993E-04;
    weight( 13) =   0.39802902551008580387116174900106E-05;
    weight( 14) =   0.45931839841801061673729694510289E-06;
    weight( 15) =   0.43075545187731100930131457465897E-07;
    weight( 16) =   0.32551249938271570855175749257884E-08;
    weight( 17) =   0.19620246675410594996247151593142E-09;
    weight( 18) =   0.93190499086617587129534716431331E-11;
    weight( 19) =   0.34377541819411620520312597898311E-12;
    weight( 20) =   0.96795247130446716997405035776206E-14;
    weight( 21) =   0.20368066110115247398010624219291E-15;
    weight( 22) =   0.31212687280713526831765358632585E-17;
    weight( 23) =   0.33729581704161052453395678308350E-19;
    weight( 24) =   0.24672796386616696011038363242541E-21;
    weight( 25) =   0.11582201904525643634834564576593E-23;
    weight( 26) =   0.32472922591425422434798022809020E-26;
    weight( 27) =   0.49143017308057432740820076259666E-29;
    weight( 28) =   0.34500071104808394132223135953806E-32;
    weight( 29) =   0.87663710117162041472932760732881E-36;
    weight( 30) =   0.50363643921161490411297172316582E-40;
    weight( 31) =   0.19909984582531456482439549080330E-45;

  elseif ( order == 32 )

    xtab( 1) =   0.04448936583326720E+00;
    xtab( 2) =   0.2345261095196173E+00;
    xtab( 3) =   0.5768846293018861E+00;
    xtab( 4) =    1.072448753817818E+00;
    xtab( 5) =    1.722408776444646E+00;
    xtab( 6) =    2.528336706425794E+00;
    xtab( 7) =    3.492213273021993E+00;
    xtab( 8) =    4.616456769749767E+00;
    xtab( 9) =    5.903958504174245E+00;
    xtab(10) =    7.358126733186242E+00;
    xtab(11) =    8.982940924212595E+00;
    xtab(12) =    10.78301863253997E+00;
    xtab(13) =    12.76369798674272E+00;
    xtab(14) =    14.93113975552256E+00;
    xtab(15) =    17.29245433671532E+00;
    xtab(16) =    19.85586094033605E+00;
    xtab(17) =    22.63088901319678E+00;
    xtab(18) =    25.62863602245925E+00;
    xtab(19) =    28.86210181632347E+00;
    xtab(20) =    32.34662915396473E+00;
    xtab(21) =    36.10049480575197E+00;
    xtab(22) =    40.14571977153944E+00;
    xtab(23) =    44.50920799575494E+00;
    xtab(24) =    49.22439498730864E+00;
    xtab(25) =    54.33372133339691E+00;
    xtab(26) =    59.89250916213402E+00;
    xtab(27) =    65.97537728793505E+00;
    xtab(28) =    72.68762809066271E+00;
    xtab(29) =    80.18744697791352E+00;
    xtab(30) =    88.73534041789240E+00;
    xtab(31) =    98.82954286828397E+00;
    xtab(32) =    111.7513980979377E+00;
 
    weight( 1) =   0.1092183419523677E+00;
    weight( 2) =   0.2104431079388168E+00;
    weight( 3) =   0.2352132296698471E+00;
    weight( 4) =   0.1959033359728783E+00;
    weight( 5) =   0.1299837862860714E+00;
    weight( 6) =   0.7057862386571766E-01;
    weight( 7) =   0.3176091250917535E-01;
    weight( 8) =   0.1191821483483859E-01;
    weight( 9) =   0.3738816294611523E-02;
    weight(10) =   0.9808033066149506E-03;
    weight(11) =   0.2148649188013644E-03;
    weight(12) =   0.3920341967987944E-04;
    weight(13) =   0.5934541612868650E-05;
    weight(14) =   0.7416404578667550E-06;
    weight(15) =   0.7604567879120788E-07;
    weight(16) =   0.6350602226625806E-08;
    weight(17) =   0.4281382971040924E-09;
    weight(18) =   0.2305899491891338E-10;
    weight(19) =   0.9799379288727102E-12;
    weight(20) =   0.3237801657729275E-13;
    weight(21) =   0.8171823443420754E-15;
    weight(22) =   0.1542133833393825E-16;
    weight(23) =   0.2119792290163629E-18;
    weight(24) =   0.2054429673788038E-20;
    weight(25) =   0.1346982586637393E-22;
    weight(26) =   0.5661294130397363E-25;
    weight(27) =   0.1418560545463052E-27;
    weight(28) =   0.1913375494454211E-30;
    weight(29) =   0.1192248760098218E-33;
    weight(30) =   0.2671511219240120E-37;
    weight(31) =   0.1338616942106269E-41;
    weight(32) =   0.4510536193898970E-47;

  elseif ( order == 63 )

    xtab(  1) =   0.22768893732576153785994330248562E-01;
    xtab(  2) =   0.11998325242727824715771416426383E+00;
    xtab(  3) =   0.29494185444770149577427738517405E+00;
    xtab(  4) =   0.54779087896237725363865073775856E+00;
    xtab(  5) =   0.87869061179931901673895567052285E+00;
    xtab(  6) =    1.2878464335919706302309207788611E+00;
    xtab(  7) =    1.7755123815388553763979463268728E+00;
    xtab(  8) =    2.3419925567085989256055628337716E+00;
    xtab(  9) =    2.9876423223246473939976731053629E+00;
    xtab( 10) =    3.7128695992018000346299637413422E+00;
    xtab( 11) =    4.5181363349503584391105568561550E+00;
    xtab( 12) =    5.4039601781825946286902599782736E+00;
    xtab( 13) =    6.3709163787865330220392250891777E+00;
    xtab( 14) =    7.4196399339311711154888493199004E+00;
    xtab( 15) =    8.5508280008403328312589048722235E+00;
    xtab( 16) =    9.7652425999245366807004592977996E+00;
    xtab( 17) =    11.063713635140661736220550410604E+00;
    xtab( 18) =    12.447142262356492749798687569289E+00;
    xtab( 19) =    13.916504641057818562912967008183E+00;
    xtab( 20) =    15.472856110036296424777143607779E+00;
    xtab( 21) =    17.117335833863588753116900303886E+00;
    xtab( 22) =    18.851171974154856850873483787506E+00;
    xtab( 23) =    20.675687448056515660377265667433E+00;
    xtab( 24) =    22.592306346311528381292277759986E+00;
    xtab( 25) =    24.602561094972638883700642760037E+00;
    xtab( 26) =    26.708100458737343969779087998829E+00;
    xtab( 27) =    28.910698500451382640177718103234E+00;
    xtab( 28) =    31.212264631175912885477773820802E+00;
    xtab( 29) =    33.614854909101154836598842888345E+00;
    xtab( 30) =    36.120684774484823056306328740825E+00;
    xtab( 31) =    38.732143442933582145626041607663E+00;
    xtab( 32) =    41.451810222318741191114726181363E+00;
    xtab( 33) =    44.282473071479233839358857134636E+00;
    xtab( 34) =    47.227149784295686898935095231536E+00;
    xtab( 35) =    50.289112264240695761749021839419E+00;
    xtab( 36) =    53.471914456788652808348280619542E+00;
    xtab( 37) =    56.779424636342062213099781057119E+00;
    xtab( 38) =    60.215862909019862886417550114424E+00;
    xtab( 39) =    63.785845004235974631701139601836E+00;
    xtab( 40) =    67.494433702293885830374325695045E+00;
    xtab( 41) =    71.347199604295266286654803376075E+00;
    xtab( 42) =    75.350293425653234254290504744279E+00;
    xtab( 43) =    79.510532629986309149555391354778E+00;
    xtab( 44) =    83.835506080872257843339817658508E+00;
    xtab( 45) =    88.333701570354369086112766326498E+00;
    xtab( 46) =    93.014662728558547405303399037100E+00;
    xtab( 47) =    97.889184147578140043386727677112E+00;
    xtab( 48) =    102.96955690741381650783952746778E+00;
    xtab( 49) =    108.26988161961595392226350967206E+00;
    xtab( 50) =    113.80647350287462738934485955901E+00;
    xtab( 51) =    119.59839538830458666962452963285E+00;
    xtab( 52) =    125.66817255856119431291196303280E+00;
    xtab( 53) =    132.04277272091165746585590583045E+00;
    xtab( 54) =    138.75498418103789078167590567526E+00;
    xtab( 55) =    145.84541318313540358283994248439E+00;
    xtab( 56) =    153.36548459497863623710815962660E+00;
    xtab( 57) =    161.38215194813761243562172669592E+00;
    xtab( 58) =    169.98570600665839438795175301156E+00;
    xtab( 59) =    179.30366247401580910251827858515E+00;
    xtab( 60) =    189.52789596532475473668721332981E+00;
    xtab( 61) =    200.97521159924656741628671841018E+00;
    xtab( 62) =    214.25368536638788642698056296400E+00;
    xtab( 63) =    230.93465747089703971246562985079E+00;
 
    weight(  1) =   0.57118633213868979811587283390476E-01;
    weight(  2) =   0.12067476090640395283319932036351E+00;
    weight(  3) =   0.15925001096581873723870561096472E+00;
    weight(  4) =   0.16875178327560799234596192963585E+00;
    weight(  5) =   0.15366641977668956696193711310131E+00;
    weight(  6) =   0.12368770614716481641086652261948E+00;
    weight(  7) =   0.89275098854848671545279150057422E-01;
    weight(  8) =   0.58258485446105944957571825725160E-01;
    weight(  9) =   0.34546657545992580874717085812508E-01;
    weight( 10) =   0.18675685985714656798286552591203E-01;
    weight( 11) =   0.92233449044093536528490075241649E-02;
    weight( 12) =   0.41671250684839592762582663470209E-02;
    weight( 13) =   0.17238120299900582715386728541955E-02;
    weight( 14) =   0.65320845029716311169340559359043E-03;
    weight( 15) =   0.22677644670909586952405173207471E-03;
    weight( 16) =   0.72127674154810668410750270234861E-04;
    weight( 17) =   0.21011261180466484598811536851241E-04;
    weight( 18) =   0.56035500893357212749181536071292E-05;
    weight( 19) =   0.13673642785604888017836641282292E-05;
    weight( 20) =   0.30507263930195817240736097189550E-06;
    weight( 21) =   0.62180061839309763559981775409241E-07;
    weight( 22) =   0.11566529551931711260022448996296E-07;
    weight( 23) =   0.19614588267565478081534781863335E-08;
    weight( 24) =   0.30286171195709411244334756404054E-09;
    weight( 25) =   0.42521344539400686769012963452599E-10;
    weight( 26) =   0.54202220578073819334698791381873E-11;
    weight( 27) =   0.62627306838597672554166850420603E-12;
    weight( 28) =   0.65474443156573322992307089591924E-13;
    weight( 29) =   0.61815575808729181846302500000047E-14;
    weight( 30) =   0.52592721363507381404263991342633E-15;
    weight( 31) =   0.40230920092646484015391506025408E-16;
    weight( 32) =   0.27600740511819536505013824207729E-17;
    weight( 33) =   0.16936946756968296053322009855265E-18;
    weight( 34) =   0.92689146872177087314963772462726E-20;
    weight( 35) =   0.45093739060365632939780140603959E-21;
    weight( 36) =   0.19435162876132376573629962695374E-22;
    weight( 37) =   0.73926270895169207037999639194513E-24;
    weight( 38) =   0.24714364154434632615980126000066E-25;
    weight( 39) =   0.72288649446741597655145390616476E-27;
    weight( 40) =   0.18407617292614039362985209905608E-28;
    weight( 41) =   0.40583498566841960105759537058880E-30;
    weight( 42) =   0.77000496416438368114463925286343E-32;
    weight( 43) =   0.12488505764999334328843314866038E-33;
    weight( 44) =   0.17185000226767010697663950619912E-35;
    weight( 45) =   0.19896372636672396938013975755522E-37;
    weight( 46) =   0.19199671378804058267713164416870E-39;
    weight( 47) =   0.15278588285522166920459714708240E-41;
    weight( 48) =   0.99054752688842142955854138884590E-44;
    weight( 49) =   0.51597523673029211884228858692990E-46;
    weight( 50) =   0.21249846664084111245693912887783E-48;
    weight( 51) =   0.67903852766852910591172042494884E-51;
    weight( 52) =   0.16466654148296177467908300517887E-53;
    weight( 53) =   0.29509065402691055027053659375033E-56;
    weight( 54) =   0.37838420647571051984882241014675E-59;
    weight( 55) =   0.33358130068542431878174667995217E-62;
    weight( 56) =   0.19223461022273880981363303073329E-65;
    weight( 57) =   0.67812696961083016872779388922288E-69;
    weight( 58) =   0.13404752802440604607620468935693E-72;
    weight( 59) =   0.13109745101805029757648048223928E-76;
    weight( 60) =   0.52624863881401787388694579143866E-81;
    weight( 61) =   0.63780013856587414257760666006511E-86;
    weight( 62) =   0.12997078942372924566347473916943E-91;
    weight( 63) =   0.10008511496968754063443740168421E-98;
 
  elseif ( order == 64 )

    xtab( 1) =   0.02241587414670647E+00;
    xtab( 2) =   0.1181225120967683E+00;
    xtab( 3) =   0.2903657440180379E+00;
    xtab( 4) =   0.5392862212279776E+00;
    xtab( 5) =   0.8650370046481133E+00;
    xtab( 6) =    1.267814040775243E+00;
    xtab( 7) =    1.747859626059437E+00;
    xtab( 8) =    2.305463739307508E+00;
    xtab( 9) =    2.940965156725251E+00;
    xtab(10) =    3.654752650207290E+00;
    xtab(11) =    4.447266343313095E+00;
    xtab(12) =    5.318999254496392E+00;
    xtab(13) =    6.270499046923656E+00;
    xtab(14) =    7.302370002587396E+00;
    xtab(15) =    8.415275239483025E+00;
    xtab(16) =    9.609939192796110E+00;
    xtab(17) =    10.88715038388637E+00;
    xtab(18) =    12.24776450424430E+00;
    xtab(19) =    13.69270784554751E+00;
    xtab(20) =    15.22298111152473E+00;
    xtab(21) =    16.83966365264874E+00;
    xtab(22) =    18.54391817085919E+00;
    xtab(23) =    20.33699594873024E+00;
    xtab(24) =    22.22024266595088E+00;
    xtab(25) =    24.19510487593326E+00;
    xtab(26) =    26.26313722711849E+00;
    xtab(27) =    28.42601052750102E+00;
    xtab(28) =    30.68552076752597E+00;
    xtab(29) =    33.04359923643783E+00;
    xtab(30) =    35.50232389114121E+00;
    xtab(31) =    38.06393216564647E+00;
    xtab(32) =    40.73083544445863E+00;
    xtab(33) =    43.50563546642153E+00;
    xtab(34) =    46.39114297861619E+00;
    xtab(35) =    49.39039902562469E+00;
    xtab(36) =    52.50669934134631E+00;
    xtab(37) =    55.74362241327838E+00;
    xtab(38) =    59.10506191901710E+00;
    xtab(39) =    62.59526440015139E+00;
    xtab(40) =    66.21887325124756E+00;
    xtab(41) =    69.98098037714684E+00;
    xtab(42) =    73.88718723248296E+00;
    xtab(43) =    77.94367743446313E+00;
    xtab(44) =    82.15730377831930E+00;
    xtab(45) =    86.53569334945652E+00;
    xtab(46) =    91.08737561313309E+00;
    xtab(47) =    95.82194001552074E+00;
    xtab(48) =    100.7502319695140E+00;
    xtab(49) =    105.8845994687999E+00;
    xtab(50) =    111.2392075244396E+00;
    xtab(51) =    116.8304450513065E+00;
    xtab(52) =    122.6774602685386E+00;
    xtab(53) =    128.8028787692377E+00;
    xtab(54) =    135.2337879495258E+00;
    xtab(55) =    142.0031214899315E+00;
    xtab(56) =    149.1516659000494E+00;
    xtab(57) =    156.7310751326712E+00;
    xtab(58) =    164.8086026551505E+00;
    xtab(59) =    173.4749468364243E+00;
    xtab(60) =    182.8582046914315E+00;
    xtab(61) =    193.1511360370729E+00;
    xtab(62) =    204.6720284850595E+00;
    xtab(63) =    218.0318519353285E+00;
    xtab(64) =    234.8095791713262E+00;
 
    weight( 1) =   0.05625284233891167E+00;
    weight( 2) =   0.1190239873124324E+00;
    weight( 3) =   0.1574964038621455E+00;
    weight( 4) =   0.1675470504157663E+00;
    weight( 5) =   0.1533528557792350E+00;
    weight( 6) =   0.1242210536093002E+00;
    weight( 7) =   0.9034230098648614E-01;
    weight( 8) =   0.5947775576835436E-01;
    weight( 9) =   0.3562751890403622E-01;
    weight(10) =   0.1948041043116633E-01;
    weight(11) =   0.9743594899381976E-02;
    weight(12) =   0.4464310364166273E-02;
    weight(13) =   0.1875359581323072E-02;
    weight(14) =   0.7226469815750038E-03;
    weight(15) =   0.2554875328334950E-03;
    weight(16) =   0.8287143534396861E-04;
    weight(17) =   0.2465686396788580E-04;
    weight(18) =   0.6726713878829671E-05;
    weight(19) =   0.1681785369964099E-05;
    weight(20) =   0.3850812981546702E-06;
    weight(21) =   0.8068728040990521E-07;
    weight(22) =   0.1545723706757701E-07;
    weight(23) =   0.2704480147617480E-08;
    weight(24) =   0.4316775475427202E-09;
    weight(25) =   0.6277752541761450E-10;
    weight(26) =   0.8306317376288941E-11;
    weight(27) =   0.9984031787220194E-12;
    weight(28) =   0.1088353887116663E-12;
    weight(29) =   0.1074017403441591E-13;
    weight(30) =   0.9575737231574444E-15;
    weight(31) =   0.7697028023648578E-16;
    weight(32) =   0.5564881137454025E-17;
    weight(33) =   0.3609756409010444E-18;
    weight(34) =   0.2095095369548946E-19;
    weight(35) =   0.1084793301097549E-20;
    weight(36) =   0.4994699486363812E-22;
    weight(37) =   0.2037836974598821E-23;
    weight(38) =   0.7339537564278870E-25;
    weight(39) =   0.2323783082198704E-26;
    weight(40) =   0.6438234706908816E-28;
    weight(41) =   0.1553121095788267E-29;
    weight(42) =   0.3244250092019555E-31;
    weight(43) =   0.5832386267836164E-33;
    weight(44) =   0.8963254833102871E-35;
    weight(45) =   0.1168703989550731E-36;
    weight(46) =   0.1282055984359980E-38;
    weight(47) =   0.1172094937404993E-40;
    weight(48) =   0.8835339672328639E-43;
    weight(49) =   0.5424955590306203E-45;
    weight(50) =   0.2675542666678911E-47;
    weight(51) =   0.1042917031411369E-49;
    weight(52) =   0.3152902351957753E-52;
    weight(53) =   0.7229541910647446E-55;
    weight(54) =   0.1224235301230072E-57;
    weight(55) =   0.1482168504901918E-60;
    weight(56) =   0.1232519348814539E-63;
    weight(57) =   0.6691499004571088E-67;
    weight(58) =   0.2220465941850503E-70;
    weight(59) =   0.4120946094738935E-74;
    weight(60) =   0.3774399061896465E-78;
    weight(61) =   0.1414115052917608E-82;
    weight(62) =   0.1591833064041367E-87;
    weight(63) =   0.2989484348860626E-93;
    weight(64) =   0.02089063508436960E-99;

  elseif ( order == 127 )

    xtab(  1) =   0.11339635298518611691893169631306E-01;
    xtab(  2) =   0.59749753435726620281348237057387E-01;
    xtab(  3) =   0.14685098690746167612388223687431E+00;
    xtab(  4) =   0.27267590735859553131378008278900E+00;
    xtab(  5) =   0.43724600644192665554577035869932E+00;
    xtab(  6) =   0.64058688222566929533576416399983E+00;
    xtab(  7) =   0.88272968639058364481487653650042E+00;
    xtab(  8) =    1.1637114160166537661560584700951E+00;
    xtab(  9) =    1.4835750152834613891313584861012E+00;
    xtab( 10) =    1.8423694351613565380686320809853E+00;
    xtab( 11) =    2.2401496839579024244513315656522E+00;
    xtab( 12) =    2.6769768780141303692167869961238E+00;
    xtab( 13) =    3.1529182957082825565771508308846E+00;
    xtab( 14) =    3.6680474360304752540226339926515E+00;
    xtab( 15) =    4.2224440823301888455977876667425E+00;
    xtab( 16) =    4.8161943715870502475665535087286E+00;
    xtab( 17) =    5.4493908694559416755862178908416E+00;
    xtab( 18) =    6.1221326512997254193944584763155E+00;
    xtab( 19) =    6.8345253894122668112237994973336E+00;
    xtab( 20) =    7.5866814466367472174205986836847E+00;
    xtab( 21) =    8.3787199765932725254842120659452E+00;
    xtab( 22) =    9.2107670307426558777922506102445E+00;
    xtab( 23) =    10.082955672528643809166439353647E+00;
    xtab( 24) =    10.995426098858125429803147358780E+00;
    xtab( 25) =    11.948325769197725997610605127857E+00;
    xtab( 26) =    12.941809542585531053723381098192E+00;
    xtab( 27) =    13.976039822878506520014405668679E+00;
    xtab( 28) =    15.051186712579523631574796365435E+00;
    xtab( 29) =    16.167428175612852922977395051768E+00;
    xtab( 30) =    17.324950209443673446561163712616E+00;
    xtab( 31) =    18.523947026965688560811711309349E+00;
    xtab( 32) =    19.764621248611504104071669386884E+00;
    xtab( 33) =    21.047184105173183606877044020054E+00;
    xtab( 34) =    22.371855651855542817648123918101E+00;
    xtab( 35) =    23.738864994122497183652313788712E+00;
    xtab( 36) =    25.148450525937368234077278385644E+00;
    xtab( 37) =    26.600860181041749607253384279755E+00;
    xtab( 38) =    28.096351697964619201753961292129E+00;
    xtab( 39) =    29.635192899504178910610227138642E+00;
    xtab( 40) =    31.217661987479759144214467152615E+00;
    xtab( 41) =    32.844047853610430460522951341338E+00;
    xtab( 42) =    34.514650407441149149105635947422E+00;
    xtab( 43) =    36.229780922306804019615388508885E+00;
    xtab( 44) =    37.989762400399956435968780140278E+00;
    xtab( 45) =    39.794929958089961778396437141707E+00;
    xtab( 46) =    41.645631232730180705153990897484E+00;
    xtab( 47) =    43.542226812286859549950892993822E+00;
    xtab( 48) =    45.485090689228791137996151336673E+00;
    xtab( 49) =    47.474610740231964719468766599146E+00;
    xtab( 50) =    49.511189233379087716728884584381E+00;
    xtab( 51) =    51.595243364671244443182771266934E+00;
    xtab( 52) =    53.727205825819316758288140069145E+00;
    xtab( 53) =    55.907525405447553305830605991732E+00;
    xtab( 54) =    58.136667626022439197077526025660E+00;
    xtab( 55) =    60.415115419018590295707192053805E+00;
    xtab( 56) =    62.743369841051809700207126742685E+00;
    xtab( 57) =    65.121950833949996311956025417139E+00;
    xtab( 58) =    67.551398031997886314411872443149E+00;
    xtab( 59) =    70.032271619884584511229871192030E+00;
    xtab( 60) =    72.565153245206849090888669416801E+00;
    xtab( 61) =    75.150646989739935299354362325096E+00;
    xtab( 62) =    77.789380404085816000647405462136E+00;
    xtab( 63) =    80.482005610750729205803962926758E+00;
    xtab( 64) =    83.229200481195914886796120019048E+00;
    xtab( 65) =    86.031669892953582966798238732643E+00;
    xtab( 66) =    88.890147073512051099652518544282E+00;
    xtab( 67) =    91.805395038358177994971250170499E+00;
    xtab( 68) =    94.778208131331583205387031034825E+00;
    xtab( 69) =    97.809413676305116411054110115424E+00;
    xtab( 70) =    100.89987375017285940371939762172E+00;
    xtab( 71) =    104.05048708821598934704076845022E+00;
    xtab( 72) =    107.26219113414600428423116401414E+00;
    xtab( 73) =    110.53596424851500530602771351277E+00;
    xtab( 74) =    113.87282809075839485348376187652E+00;
    xtab( 75) =    117.27385019192517774095477886379E+00;
    xtab( 76) =    120.74014673718880106173978002719E+00;
    xtab( 77) =    124.27288557955698354259506446928E+00;
    xtab( 78) =    127.87328950885942645093841745425E+00;
    xtab( 79) =    131.54263980314366921809377742137E+00;
    xtab( 80) =    135.28228009311836970132738106369E+00;
    xtab( 81) =    139.09362057432970013964422086977E+00;
    xtab( 82) =    142.97814260643601776808227753574E+00;
    xtab( 83) =    146.93740374437366549441080969072E+00;
    xtab( 84) =    150.97304325252187127492511437460E+00;
    xtab( 85) =    155.08678816034612572229641420609E+00;
    xtab( 86) =    159.28045992663288235401956989889E+00;
    xtab( 87) =    163.55598178957571104015967182053E+00;
    xtab( 88) =    167.91538689194360134245547184721E+00;
    xtab( 89) =    172.36082728473812536838156191681E+00;
    xtab( 90) =    176.89458392960192176311674993508E+00;
    xtab( 91) =    181.51907784036813069227528834025E+00;
    xtab( 92) =    186.23688252828112373861202530357E+00;
    xtab( 93) =    191.05073794450929196790836610789E+00;
    xtab( 94) =    195.96356614879879837839002542988E+00;
    xtab( 95) =    200.97848897600025153696475526130E+00;
    xtab( 96) =    206.09884802468871112127283042753E+00;
    xtab( 97) =    211.32822735671655260572377256981E+00;
    xtab( 98) =    216.67047937658230323477089465777E+00;
    xtab( 99) =    222.12975445929687246267304963754E+00;
    xtab(100) =    227.71053502072232419089132431317E+00;
    xtab(101) =    233.41767488282602453367775322563E+00;
    xtab(102) =    239.25644498830308620018749667089E+00;
    xtab(103) =    245.23258677871567172531254018984E+00;
    xtab(104) =    251.35237488718128030005500991754E+00;
    xtab(105) =    257.62269123792061413076191882313E+00;
    xtab(106) =    264.05111322908240551754377241831E+00;
    xtab(107) =    270.64601945722796749299111718606E+00;
    xtab(108) =    277.41671750163651071798388218104E+00;
    xtab(109) =    284.37359974220870326674402873120E+00;
    xtab(110) =    291.52833521346495719581282021650E+00;
    xtab(111) =    298.89410837028248600878895615414E+00;
    xtab(112) =    306.48591978262611320418112423947E+00;
    xtab(113) =    314.32096986471177487400007507615E+00;
    xtab(114) =    322.41915589128679683349440361344E+00;
    xtab(115) =    330.80372663802405651933847334878E+00;
    xtab(116) =    339.50216127832433747735367595958E+00;
    xtab(117) =    348.54737559472697355480761787441E+00;
    xtab(118) =    357.97942028029845454049007443090E+00;
    xtab(119) =    367.84794520076004578858341422871E+00;
    xtab(120) =    378.21590623135532818332979188889E+00;
    xtab(121) =    389.16539141251004101579475325153E+00;
    xtab(122) =    400.80729331451702589996361286427E+00;
    xtab(123) =    413.29853681779384418008260081859E+00;
    xtab(124) =    426.87579153663675538288509017051E+00;
    xtab(125) =    441.93085485310841412460309271842E+00;
    xtab(126) =    459.21804639888429981971267313224E+00;
    xtab(127) =    480.69378263388373859704269229304E+00;
 
    weight(  1) =   0.28773246692000124355770010301506E-01;
    weight(  2) =   0.63817468175134649363480949265236E-01;
    weight(  3) =   0.91919669721570571389864194652717E-01;
    weight(  4) =   0.11054167914413766381245463002967E+00;
    weight(  5) =   0.11879771633375850188328329422643E+00;
    weight(  6) =   0.11737818530052695148804451630074E+00;
    weight(  7) =   0.10819305984180551488335145581193E+00;
    weight(  8) =   0.93827075290489628080377261401107E-01;
    weight(  9) =   0.76966450960588843995822485928431E-01;
    weight( 10) =   0.59934903912939714332570730063476E-01;
    weight( 11) =   0.44417742073889001371708316272923E-01;
    weight( 12) =   0.31385080966252320983009372215062E-01;
    weight( 13) =   0.21172316041924506411370709025015E-01;
    weight( 14) =   0.13650145364230541652171185564626E-01;
    weight( 15) =   0.84172852710599172279366657385445E-02;
    weight( 16) =   0.49674990059882760515912858620175E-02;
    weight( 17) =   0.28069903895001884631961957446400E-02;
    weight( 18) =   0.15192951003941952460445341057817E-02;
    weight( 19) =   0.78789028751796084086217287140548E-03;
    weight( 20) =   0.39156751064868450584507324648999E-03;
    weight( 21) =   0.18652434268825860550093566260060E-03;
    weight( 22) =   0.85173160415576621908809828160247E-04;
    weight( 23) =   0.37285639197853037712145321577724E-04;
    weight( 24) =   0.15648416791712993947447805296768E-04;
    weight( 25) =   0.62964340695224829035692735524979E-05;
    weight( 26) =   0.24288929711328724574541379938222E-05;
    weight( 27) =   0.89824607890051007201922871545035E-06;
    weight( 28) =   0.31844174740760353710742966328091E-06;
    weight( 29) =   0.10821272905566839211861807542741E-06;
    weight( 30) =   0.35245076750635536015902779085340E-07;
    weight( 31) =   0.11001224365719347407063839761738E-07;
    weight( 32) =   0.32904079616717932125329343003261E-08;
    weight( 33) =   0.94289145237889976419772700772988E-09;
    weight( 34) =   0.25882578904668318184050195309296E-09;
    weight( 35) =   0.68047437103370762630942259017560E-10;
    weight( 36) =   0.17131398805120837835399564475632E-10;
    weight( 37) =   0.41291744524052865469443922304935E-11;
    weight( 38) =   0.95264189718807273220707664873469E-12;
    weight( 39) =   0.21032604432442425932962942047474E-12;
    weight( 40) =   0.44427151938729352860940434285789E-13;
    weight( 41) =   0.89760500362833703323319846405449E-14;
    weight( 42) =   0.17341511407769287074627948346848E-14;
    weight( 43) =   0.32028099548988356631494379835210E-15;
    weight( 44) =   0.56531388950793682022660742095189E-16;
    weight( 45) =   0.95329672799026591234588044025896E-17;
    weight( 46) =   0.15353453477310142565288509437552E-17;
    weight( 47) =   0.23608962179467365686057842132176E-18;
    weight( 48) =   0.34648742794456611332193876653230E-19;
    weight( 49) =   0.48515241897086461320126957663545E-20;
    weight( 50) =   0.64786228633519813428137373790678E-21;
    weight( 51) =   0.82476020965403242936448553126316E-22;
    weight( 52) =   0.10005361880214719793491658282977E-22;
    weight( 53) =   0.11561395116207304954233181263632E-23;
    weight( 54) =   0.12719342731167922655612134264961E-24;
    weight( 55) =   0.13316584714165372967340004160814E-25;
    weight( 56) =   0.13261218454678944033646108509198E-26;
    weight( 57) =   0.12554995447643949807286074138324E-27;
    weight( 58) =   0.11294412178579462703240913107219E-28;
    weight( 59) =   0.96491020279562119228500608131696E-30;
    weight( 60) =   0.78241846768302099396733076955632E-31;
    weight( 61) =   0.60181503542219626658249939076636E-32;
    weight( 62) =   0.43882482704961741551510518054138E-33;
    weight( 63) =   0.30314137647517256304035802501863E-34;
    weight( 64) =   0.19826016543944539545224676057020E-35;
    weight( 65) =   0.12267623373665926559013654872402E-36;
    weight( 66) =   0.71763931692508888943812834967620E-38;
    weight( 67) =   0.39659378833836963584113716149270E-39;
    weight( 68) =   0.20688970553868040099581951696677E-40;
    weight( 69) =   0.10179587017979517245268418427523E-41;
    weight( 70) =   0.47200827745986374625714293679649E-43;
    weight( 71) =   0.20606828985553374825744353490744E-44;
    weight( 72) =   0.84627575907305987245899032156188E-46;
    weight( 73) =   0.32661123687088798658026998931647E-47;
    weight( 74) =   0.11833939207883162380564134612682E-48;
    weight( 75) =   0.40211209123895013807243250164050E-50;
    weight( 76) =   0.12799824394111125389430292847476E-51;
    weight( 77) =   0.38123877747548846504399051365162E-53;
    weight( 78) =   0.10612057542701156767898551949650E-54;
    weight( 79) =   0.27571446947200403594113572720812E-56;
    weight( 80) =   0.66772544240928492881306904862856E-58;
    weight( 81) =   0.15052438383868234954068178600268E-59;
    weight( 82) =   0.31538986800113758526689068500772E-61;
    weight( 83) =   0.61326614299483180785237418887960E-63;
    weight( 84) =   0.11048510030324810567549119229368E-64;
    weight( 85) =   0.18410563538091348076979665543900E-66;
    weight( 86) =   0.28323926570052832195543883237652E-68;
    weight( 87) =   0.40154409843763655508670978777418E-70;
    weight( 88) =   0.52351530215683708779772201956106E-72;
    weight( 89) =   0.62634476665005100555787696642851E-74;
    weight( 90) =   0.68612210535666530365348093803922E-76;
    weight( 91) =   0.68651298840956019297134099761855E-78;
    weight( 92) =   0.62581388433728084867318704240915E-80;
    weight( 93) =   0.51833271237514904046803469968027E-82;
    weight( 94) =   0.38893621571918443533108973497673E-84;
    weight( 95) =   0.26357711379476932781525533730623E-86;
    weight( 96) =   0.16078851293917979699005509638883E-88;
    weight( 97) =   0.87978042070968939637972577886624E-91;
    weight( 98) =   0.43013405077495109903408697802188E-93;
    weight( 99) =   0.18713435881342838527144321803729E-95;
    weight(100) =   0.72125744708060471675805761366523E-98;
    weight(101) =   0.24508746062177874383231742333023E-100;
    weight(102) =   0.73042094619470875777647865078327E-103;
    weight(103) =   0.18983290818383463537886818579820E-105;
    weight(104) =   0.42757400244246684123093264825902E-108;
    weight(105) =   0.82894681420515755691423485228897E-111;
    weight(106) =   0.13729432219324400013067050156048E-113;
    weight(107) =   0.19265464126404973222043166489406E-116;
    weight(108) =   0.22693344503301354826140809941334E-119;
    weight(109) =   0.22209290603717355061909071271535E-122;
    weight(110) =   0.17851087685544512662856555121755E-125;
    weight(111) =   0.11630931990387164467431190485525E-128;
    weight(112) =   0.60524443584652392290952805077893E-132;
    weight(113) =   0.24729569115063528647628375096400E-135;
    weight(114) =   0.77789065006489410364997205809045E-139;
    weight(115) =   0.18409738662712607039570678274636E-142;
    weight(116) =   0.31900921131079114970179071968597E-146;
    weight(117) =   0.39179487139174199737617666077555E-150;
    weight(118) =   0.32782158394188697053774429820559E-154;
    weight(119) =   0.17793590713138888062819640128739E-158;
    weight(120) =   0.58882353408932623157467835381214E-163;
    weight(121) =   0.10957236509071169877747203273886E-167;
    weight(122) =   0.10281621114867000898285076975760E-172;
    weight(123) =   0.41704725557697758145816510853967E-178;
    weight(124) =   0.58002877720316101774638319601971E-184;
    weight(125) =   0.18873507745825517106171619101120E-190;
    weight(126) =   0.69106601826730911682786705950895E-198;
    weight(127) =   0.43506813201105855628383313334402E-207;
 
  else

    fprintf ( 1, '\n' );
    fprintf ( 1, 'LAGUERRE_SET - Fatal error!\n' );
    fprintf ( 1, '  Illegal value of NORDER = %d\n', norder );
    fprintf ( 1, '  Legal values are 1 to 20, 31, 32, 63, 64 or 127\n' );
    error ( 'LAGUERRE_SET - Fatal error!' );

  end