cbess0_cpp.txt

BESSEL PROGRAMS IN C/C

      ALFA[17]= 1.31668816545923e-0004;  ALFA[18]= 1.23667445598253e-0004;
      ALFA[19]= 1.16405271474738e-0004;  ALFA[20]= 1.09798298372713e-0004;
      ALFA[21]= 1.03772410422993e-0004;  ALFA[22]= 9.82626078369363e-0005;
      ALFA[23]= 9.32120517249503e-0005;  ALFA[24]= 8.85710852478712e-0005;
      ALFA[25]= 8.42963105715700e-0005;  ALFA[26]= 8.03497548407791e-0005;
      ALFA[27]= 7.66981345359207e-0005;  ALFA[28]= 7.33122157481778e-0005;
      ALFA[29]= 7.01662625163141e-0005;  ALFA[30]= 6.72375633790160e-0005;
      ALFA[31]= 6.93735541354589e-0004;  ALFA[32]= 2.32241745182922e-0004;
      ALFA[33]=-1.41986273556691e-0005;  ALFA[34]=-1.16444931672049e-0004;
      ALFA[35]=-1.50803558053049e-0004;  ALFA[36]=-1.55121924918096e-0004;
      ALFA[37]=-1.46809756646466e-0004;  ALFA[38]=-1.33815503867491e-0004;
      ALFA[39]=-1.19744975684254e-0004;  ALFA[40]=-1.06184319207974e-0004;
      ALFA[41]=-9.37699549891194e-0005;  ALFA[42]=-8.26923045588193e-0005;
      ALFA[43]=-7.29374348155221e-0005;  ALFA[44]=-6.44042357721016e-0005;
      ALFA[45]=-5.69611566009369e-0005;  ALFA[46]=-5.04731044303562e-0005;
      ALFA[47]=-4.48134868008883e-0005;  ALFA[48]=-3.98688727717599e-0005;
      ALFA[49]=-3.55400532972043e-0005;  ALFA[50]=-3.17414256609023e-0005;
      ALFA[51]=-2.83996793904175e-0005;  ALFA[52]=-2.54522720634871e-0005;
      ALFA[53]=-2.28459297164725e-0005;  ALFA[54]=-2.05352753106481e-0005;
      ALFA[55]=-1.84816217627666e-0005;  ALFA[56]=-1.66519330021394e-0005;
      ALFA[57]=-1.50179412980119e-0005;  ALFA[58]=-1.35554031379041e-0005;
      ALFA[59]=-1.22434746473858e-0005;  ALFA[60]=-1.10641884811308e-0005;
      ALFA[61]=-3.54211971457744e-0004;  ALFA[62]=-1.56161263945159e-0004;
      ALFA[63]= 3.04465503594936e-0005;  ALFA[64]= 1.30198655773243e-0004;
      ALFA[65]= 1.67471106699712e-0004;  ALFA[66]= 1.70222587683593e-0004;
      ALFA[67]= 1.56501427608595e-0004;  ALFA[68]= 1.36339170977445e-0004;
      ALFA[69]= 1.14886692029825e-0004;  ALFA[70]= 9.45869093034688e-0005;
      ALFA[71]= 7.64498419250898e-0005;  ALFA[72]= 6.07570334965197e-0005;
      ALFA[73]= 4.74394299290509e-0005;  ALFA[74]= 3.62757512005344e-0005;
      ALFA[75]= 2.69939714979225e-0005;  ALFA[76]= 1.93210938247939e-0005;
      ALFA[77]= 1.30056674793963e-0005;  ALFA[78]= 7.82620866744497e-0006;
      ALFA[79]= 3.59257485819352e-0006;  ALFA[80]= 1.44040049814252e-0007;
      ALFA[81]=-2.65396769697939e-0006;  ALFA[82]=-4.91346867098486e-0006;
      ALFA[83]=-6.72739296091248e-0006;  ALFA[84]=-8.17269379678658e-0006;
      ALFA[85]=-9.31304715093561e-0006;  ALFA[86]=-1.02011418798016e-0005;
      ALFA[87]=-1.08805962510593e-0005;  ALFA[88]=-1.13875481509604e-0005;
      ALFA[89]=-1.17519675674556e-0005;  ALFA[90]=-1.19987364870944e-0005;
      ALFA[91]= 3.78194199201773e-0004;  ALFA[92]= 2.02471952761816e-0004;
      ALFA[93]=-6.37938506318862e-0005;  ALFA[94]=-2.38598230603006e-0004;
      ALFA[95]=-3.10916256027362e-0004;  ALFA[96]=-3.13680115247576e-0004;
      ALFA[97]=-2.78950273791323e-0004;  ALFA[98]=-2.28564082619141e-0004;
      ALFA[99]=-1.75245280340847e-0004;  ALFA[100]=-1.25544063060690e-0004;
      ALFA[101]=-8.22982872820208e-0005;  ALFA[102]=-4.62860730588116e-0005;
      ALFA[103]=-1.72334302366962e-0005;  ALFA[104]= 5.60690482304602e-0006;
      ALFA[105]= 2.31395443148287e-0005;  ALFA[106]= 3.62642745856794e-0005;
      ALFA[107]= 4.58006124490189e-0005;  ALFA[108]= 5.24595294959114e-0005;
      ALFA[109]= 5.68396208545815e-0005;  ALFA[110]= 5.94349820393104e-0005;
      ALFA[111]= 6.06478527578422e-0005;  ALFA[112]= 6.08023907788436e-0005;
      ALFA[113]= 6.01577894539460e-0005;  ALFA[114]= 5.89199657344698e-0005;
      ALFA[115]= 5.72515823777593e-0005;  ALFA[116]= 5.52804375585853e-0005;
      ALFA[117]= 5.31063773802880e-0005;  ALFA[118]= 5.08069302012326e-0005;
      ALFA[119]= 4.84418647620095e-0005;  ALFA[120]= 4.60568581607475e-0005;
      ALFA[121]=-6.91141397288294e-0004;  ALFA[122]=-4.29976633058872e-0004;
      ALFA[123]= 1.83067735980039e-0004;  ALFA[124]= 6.60088147542014e-0004;
      ALFA[125]= 8.75964969951186e-0004;  ALFA[126]= 8.77335235958236e-0004;
      ALFA[127]= 7.49369585378991e-0004;  ALFA[128]= 5.63832329756981e-0004;
      ALFA[129]= 3.68059319971443e-0004;  ALFA[130]= 1.88464535514456e-0004;
      ALFA[131]= 3.70663057664904e-0005;  ALFA[132]=-8.28520220232137e-0005;
      ALFA[133]=-1.72751952869173e-0004;  ALFA[134]=-2.36314873605873e-0004;
      ALFA[135]=-2.77966150694907e-0004;  ALFA[136]=-3.02079514155457e-0004;
      ALFA[137]=-3.12594712643820e-0004;  ALFA[138]=-3.12872558758067e-0004;
      ALFA[139]=-3.05678038466324e-0004;  ALFA[140]=-2.93226470614557e-0004;
      ALFA[141]=-2.77255655582935e-0004;  ALFA[142]=-2.59103928467032e-0004;
      ALFA[143]=-2.39784014396480e-0004;  ALFA[144]=-2.20048260045423e-0004;
      ALFA[145]=-2.00443911094971e-0004;  ALFA[146]=-1.81358692210971e-0004;
      ALFA[147]=-1.63057674478657e-0004;  ALFA[148]=-1.45712672175206e-0004;
      ALFA[149]=-1.29425421983925e-0004;  ALFA[150]=-1.14245691942446e-0004;
      ALFA[151]= 1.92821964248776e-0003;  ALFA[152]= 1.35592576302022e-0003;
      ALFA[153]=-7.17858090421303e-0004;  ALFA[154]=-2.58084802575270e-0003;
      ALFA[155]=-3.49271130826168e-0003;  ALFA[156]=-3.46986299340961e-0003;
      ALFA[157]=-2.82285233351310e-0003;  ALFA[158]=-1.88103076404891e-0003;
      ALFA[159]=-8.89531718383948e-0004;  ALFA[160]= 3.87912102631035e-0006;
      ALFA[161]= 7.28688540119691e-0004;  ALFA[162]= 1.26566373053458e-0003;
      ALFA[163]= 1.62518158372674e-0003;  ALFA[164]= 1.83203153216373e-0003;
      ALFA[165]= 1.91588388990528e-0003;  ALFA[166]= 1.90588846755546e-0003;
      ALFA[167]= 1.82798982421826e-0003;  ALFA[168]= 1.70389506421122e-0003;
      ALFA[169]= 1.55097127171098e-0003;  ALFA[170]= 1.38261421852276e-0003;
      ALFA[171]= 1.20881424230065e-0003;  ALFA[172]= 1.03676532638345e-0003;
      ALFA[173]= 8.71437918068619e-0004;  ALFA[174]= 7.16080155297701e-0004;
      ALFA[175]= 5.72637002558129e-0004;  ALFA[176]= 4.42089819465802e-0004;
      ALFA[177]= 3.24724948503091e-0004;  ALFA[178]= 2.20342042730247e-0004;
      ALFA[179]= 1.28412898401354e-0004;  ALFA[180]= 4.82005924552095e-0005;

	  BETA[0]= 0.0;
      BETA[1]= 1.79988721413553e-0002;  BETA[2]= 5.59964911064388e-0003;
      BETA[3]= 2.88501402231133e-0003;  BETA[4]= 1.80096606761054e-0003;
      BETA[5]= 1.24753110589199e-0003;  BETA[6]= 9.22878876572938e-0004;
      BETA[7]= 7.14430421727287e-0004;  BETA[8]= 5.71787281789705e-0004;
      BETA[9]= 4.69431007606482e-0004;  BETA[10]= 3.93232835462917e-0004;
      BETA[11]= 3.34818889318298e-0004;  BETA[12]= 2.88952148495752e-0004;
      BETA[13]= 2.52211615549573e-0004;  BETA[14]= 2.22280580798883e-0004;
      BETA[15]= 1.97541838033063e-0004;  BETA[16]= 1.76836855019718e-0004;
      BETA[17]= 1.59316899661821e-0004;  BETA[18]= 1.44347930197334e-0004;
      BETA[19]= 1.31448068119965e-0004;  BETA[20]= 1.20245444949303e-0004;
      BETA[21]= 1.10449144504599e-0004;  BETA[22]= 1.01828770740567e-0004;
      BETA[23]= 9.41998224204238e-0005;  BETA[24]= 8.74130545753834e-0005;
      BETA[25]= 8.13466262162801e-0005;  BETA[26]= 7.59002269646219e-0005;
      BETA[27]= 7.09906300634154e-0005;  BETA[28]= 6.65482874842468e-0005;
      BETA[29]= 6.25146958969275e-0005;  BETA[30]= 5.88403394426252e-0005;
      BETA[31]=-1.49282953213429e-0003;  BETA[32]=-8.78204709546389e-0004;
      BETA[33]=-5.02916549572035e-0004;  BETA[34]=-2.94822138512746e-0004;
      BETA[35]=-1.75463996970783e-0004;  BETA[36]=-1.04008550460816e-0004;
      BETA[37]=-5.96141953046458e-0005;  BETA[38]=-3.12038929076098e-0005;
      BETA[39]=-1.26089735980230e-0005;  BETA[40]=-2.42892608575730e-0007;
      BETA[41]= 8.05996165414274e-0006;  BETA[42]= 1.36507009262147e-0005;
      BETA[43]= 1.73964125472926e-0005;  BETA[44]= 1.98672978842134e-0005;
      BETA[45]= 2.14463263790823e-0005;  BETA[46]= 2.23954659232457e-0005;
      BETA[47]= 2.28967783814713e-0005;  BETA[48]= 2.30785389811178e-0005;
      BETA[49]= 2.30321976080909e-0005;  BETA[50]= 2.28236073720349e-0005;
      BETA[51]= 2.25005881105292e-0005;  BETA[52]= 2.20981015361991e-0005;
      BETA[53]= 2.16418427448104e-0005;  BETA[54]= 2.11507649256221e-0005;
      BETA[55]= 2.06388749782171e-0005;  BETA[56]= 2.01165241997082e-0005;
      BETA[57]= 1.95913450141179e-0005;  BETA[58]= 1.90689367910437e-0005;
      BETA[59]= 1.85533719641637e-0005;  BETA[60]= 1.80475722259674e-0005;
      BETA[61]= 5.52213076721293e-0004;  BETA[62]= 4.47932581552385e-0004;
      BETA[63]= 2.79520653992021e-0004;  BETA[64]= 1.52468156198447e-0004;
      BETA[65]= 6.93271105657044e-0005;  BETA[66]= 1.76258683069991e-0005;
      BETA[67]=-1.35744996343269e-0005;  BETA[68]=-3.17972413350427e-0005;
      BETA[69]=-4.18861861696693e-0005;  BETA[70]=-4.69004889379141e-0005;
      BETA[71]=-4.87665447413787e-0005;  BETA[72]=-4.87010031186735e-0005;
      BETA[73]=-4.74755620890087e-0005;  BETA[74]=-4.55813058138628e-0005;
      BETA[75]=-4.33309644511266e-0005;  BETA[76]=-4.09230193157750e-0005;
      BETA[77]=-3.84822638603221e-0005;  BETA[78]=-3.60857167535411e-0005;
      BETA[79]=-3.37793306123367e-0005;  BETA[80]=-3.15888560772110e-0005;
      BETA[81]=-2.95269561750807e-0005;  BETA[82]=-2.75978914828336e-0005;
      BETA[83]=-2.58006174666884e-0005;  BETA[84]=-2.41308356761280e-0005;
      BETA[85]=-2.25823509518346e-0005;  BETA[86]=-2.11479656768913e-0005;
      BETA[87]=-1.98200638885295e-0005;  BETA[88]=-1.85909870801065e-0005;
      BETA[89]=-1.74532699844210e-0005;  BETA[90]=-1.63997823854498e-0005;
      BETA[91]=-4.74617796559960e-0004;  BETA[92]=-4.77864567147322e-0004;
      BETA[93]=-3.20390228067038e-0004;  BETA[94]=-1.61105016119962e-0004;
      BETA[95]=-4.25778101285435e-0005;  BETA[96]= 3.44571294294967e-0005;
      BETA[97]= 7.97092684075675e-0005;  BETA[98]= 1.03138236708272e-0004;
      BETA[99]= 1.12466775262204e-0004;  BETA[100]= 1.13103642108481e-0004;
      BETA[101]= 1.08651634848774e-0004;  BETA[102]= 1.01437951597662e-0004;
      BETA[103]= 9.29298396593364e-0005;  BETA[104]= 8.40293133016090e-0005;
      BETA[105]= 7.52727991349134e-0005;  BETA[106]= 6.69632521975731e-0005;
      BETA[107]= 5.92564547323195e-0005;  BETA[108]= 5.22169308826976e-0005;
      BETA[109]= 4.58539485165361e-0005;  BETA[110]= 4.01445513891487e-0005;
      BETA[111]= 3.50481730031328e-0005;  BETA[112]= 3.05157995034347e-0005;
      BETA[113]= 2.64956119950516e-0005;  BETA[114]= 2.29363633690998e-0005;
      BETA[115]= 1.97893056664022e-0005;  BETA[116]= 1.70091984636413e-0005;
      BETA[117]= 1.45547428261524e-0005;  BETA[118]= 1.23886640995878e-0005;
      BETA[119]= 1.04775876076583e-0005;  BETA[120]= 8.79179954978479e-0006;
      BETA[121]= 7.36465810572578e-0004;  BETA[122]= 8.72790805146194e-0004;
      BETA[123]= 6.22614862573135e-0004;  BETA[124]= 2.85998154194304e-0004;
      BETA[125]= 3.84737672879366e-0006;  BETA[126]=-1.87906003636972e-0004;
      BETA[127]=-2.97603646594555e-0004;  BETA[128]=-3.45998126832656e-0004;
      BETA[129]=-3.53382470916038e-0004;  BETA[130]=-3.35715635775049e-0004;
      BETA[131]=-3.04321124789040e-0004;  BETA[132]=-2.66722723047613e-0004;
      BETA[133]=-2.27654214122820e-0004;  BETA[134]=-1.89922611854562e-0004;
      BETA[135]=-1.55058918599094e-0004;  BETA[136]=-1.23778240761874e-0004;
      BETA[137]=-9.62926147717644e-0005;  BETA[138]=-7.25178327714425e-0005;
      BETA[139]=-5.22070028895634e-0005;  BETA[140]=-3.50347750511901e-0005;
      BETA[141]=-2.06489761035552e-0005;  BETA[142]=-8.70106096849767e-0006;
      BETA[143]= 1.13698686675100e-0006;  BETA[144]= 9.16426474122779e-0006;
      BETA[145]= 1.56477785428873e-0005;  BETA[146]= 2.08223629482467e-0005;
      BETA[147]= 2.48923381004595e-0005;  BETA[148]= 2.80340509574146e-0005;
      BETA[149]= 3.03987774629862e-0005;  BETA[150]= 3.21156731406701e-0005;
      BETA[151]=-1.80182191963886e-0003;  BETA[152]=-2.43402962938043e-0003;
      BETA[153]=-1.83422663549857e-0003;  BETA[154]=-7.62204596354010e-0004;
      BETA[155]= 2.39079475256927e-0004;  BETA[156]= 9.49266117176881e-0004;
      BETA[157]= 1.34467449701540e-0003;  BETA[158]= 1.48457495259449e-0003;
      BETA[159]= 1.44732339830618e-0003;  BETA[160]= 1.30268261285657e-0003;
      BETA[161]= 1.10351597375643e-0003;  BETA[162]= 8.86047440419792e-0004;
      BETA[163]= 6.73073208165665e-0004;  BETA[164]= 4.77603872856582e-0004;
      BETA[165]= 3.05991926358789e-0004;  BETA[166]= 1.60315694594722e-0004;
      BETA[167]= 4.00749555270613e-0005;  BETA[168]=-5.66607461635252e-0005;
      BETA[169]=-1.32506186772983e-0004;  BETA[170]=-1.90296187989614e-0004;
      BETA[171]=-2.32811450376937e-0004;  BETA[172]=-2.62628811464669e-0004;
      BETA[173]=-2.82050469867599e-0004;  BETA[174]=-2.93081563192861e-0004;
      BETA[175]=-2.97435962176317e-0004;  BETA[176]=-2.96557334239348e-0004;
      BETA[177]=-2.91647363312091e-0004;  BETA[178]=-2.83696203837734e-0004;
      BETA[179]=-2.73512317095673e-0004;  BETA[180]=-2.61750155806769e-0004;
      BETA[181]= 6.38585891212051e-0003;  BETA[182]= 9.62374215806378e-0003;
      BETA[183]= 7.61878061207001e-0003;  BETA[184]= 2.83219055545628e-0003;
      BETA[185]=-2.09841352012720e-0003;  BETA[186]=-5.73826764216626e-0003;
      BETA[187]=-7.70804244495415e-0003;  BETA[188]=-8.21011692264844e-0003;
      BETA[189]=-7.65824520346905e-0003;  BETA[190]=-6.47209729391045e-0003;
      BETA[191]=-4.99132412004966e-0003;  BETA[192]=-3.45612289713133e-0003;
      BETA[193]=-2.01785580014171e-0003;  BETA[194]=-7.59430686781961e-0004;
      BETA[195]= 2.84173631523859e-0004;  BETA[196]= 1.10891667586337e-0003;
      BETA[197]= 1.72901493872729e-0003;  BETA[198]= 2.16812590802685e-0003;
      BETA[199]= 2.45357710494540e-0003;  BETA[200]= 2.61281821058335e-0003;
      BETA[201]= 2.67141039656277e-0003;  BETA[202]= 2.65203073395980e-0003;
      BETA[203]= 2.57411652877287e-0003;  BETA[204]= 2.45389126236094e-0003;
      BETA[205]= 2.30460058071796e-0003;  BETA[206]= 2.13684837686713e-0003;
      BETA[207]= 1.95896528478871e-0003;  BETA[208]= 1.77737008679454e-0003;
      BETA[209]= 1.59690280765839e-0003;  BETA[210]= 1.42111975664439e-0003;

	  GAMA[0]= 0.0;
      GAMA[1]= 6.29960524947437e-0001;  GAMA[2]= 2.51984209978975e-0001;
      GAMA[3]= 1.54790300415656e-0001;  GAMA[4]= 1.10713062416159e-0001;
      GAMA[5]= 8.57309395527395e-0002;  GAMA[6]= 6.97161316958684e-0002;
      GAMA[7]= 5.86085671893714e-0002;  GAMA[8]= 5.04698873536311e-0002;
      GAMA[9]= 4.42600580689155e-0002;  GAMA[10]= 3.93720661543510e-0002;
      GAMA[11]= 3.54283195924455e-0002;  GAMA[12]= 3.21818857502098e-0002;
      GAMA[13]= 2.94646240791158e-0002;  GAMA[14]= 2.71581677112934e-0002;
      GAMA[15]= 2.51768272973862e-0002;  GAMA[16]= 2.34570755306079e-0002;
      GAMA[17]= 2.19508390134907e-0002;  GAMA[18]= 2.06210828235646e-0002;
      GAMA[19]= 1.94388240897881e-0002;  GAMA[20]= 1.83810633800683e-0002;
      GAMA[21]= 1.74293213231963e-0002;  GAMA[22]= 1.65685837786612e-0002;
      GAMA[23]= 1.57865285987918e-0002;  GAMA[24]= 1.50729501494096e-0002;
      GAMA[25]= 1.44193250839955e-0002;  GAMA[26]= 1.38184805735342e-0002;
      GAMA[27]= 1.32643378994277e-0002;  GAMA[28]= 1.27517121970499e-0002;
      GAMA[29]= 1.22761545318763e-0002;  GAMA[30]= 1.18338262398482e-0002;

      EX1=3.33333333333333333e-01; EX2=6.66666666666666667e-01;
      HPI=1.57079632679489662e+00; GXPI=3.14159265358979324e+00;
      THPI=4.71238898038468986;

      ZEROR=0.0; ZEROI=0.0; CONER=1.0; CONEI=0.0;

      RFNU = 1.0/FNU;
      ZBR = ZR*RFNU;
      ZBI = ZI*RFNU;
      RFNU2 = RFNU*RFNU;
/*----------------------------------------------------------------------
!     COMPUTE IN THE FOURTH QUADRANT
!---------------------------------------------------------------------*/
      FN13 = pow(FNU,EX1);
      FN23 = FN13*FN13;
      RFN13 = 1.0/FN13;
      W2R = CONER - ZBR*ZBR + ZBI*ZBI;
      W2I = CONEI - ZBR*ZBI - ZBR*ZBI;
      AW2 = ZABS(W2R,W2I);
      if (AW2 > 0.25) goto e130;
/*----------------------------------------------------------------------
!     POWER SERIES FOR CABS(W2) <= 0.25
!---------------------------------------------------------------------*/
      K = 1;
      PR[1] = CONER;
      XPI[1] = CONEI;
      SUMAR = GAMA[1];
      SUMAI = ZEROI;
      AP[1] = 1.0;
      if (AW2 < TOL) goto e20;
      for (K=2; K<=30; K++) {
        PR[K] = PR[K-1]*W2R - XPI[K-1]*W2I;
        XPI[K] = PR[K-1]*W2I + XPI[K-1]*W2R;
        SUMAR = SUMAR + PR[K]*GAMA[K];
        SUMAI = SUMAI + XPI[K]*GAMA[K];
        AP[K] = AP[K-1]*AW2;
        if (AP[K] < TOL) goto e20;
      }
      K = 30;
e20:  KMAX = K;
      ZETAR = W2R*SUMAR - W2I*SUMAI;
      ZETAI = W2R*SUMAI + W2I*SUMAR;
      *ARGR = ZETAR*FN23;
      *ARGI = ZETAI*FN23;
      ZSQRT(SUMAR, SUMAI, &ZAR, &ZAI);
      ZSQRT(W2R, W2I, &STR, &STI);
      *ZETA2R = STR*FNU;
      *ZETA2I = STI*FNU;
      STR = CONER + EX2*(ZETAR*ZAR-ZETAI*ZAI);
      STI = CONEI + EX2*(ZETAR*ZAI+ZETAI*ZAR);
      *ZETA1R = STR*(*ZETA2R) - STI*(*ZETA2I);
      *ZETA1I = STR*(*ZETA2I) + STI*(*ZETA2R);
      ZAR = ZAR + ZAR;
      ZAI = ZAI + ZAI;
      ZSQRT(ZAR, ZAI, &STR, &STI);
      *PHIR = STR*RFN13;
      *PHII = STI*RFN13;
      if (IPMTR == 1) goto e120;
/*----------------------------------------------------------------------
!     SUM SERIES FOR ASUM AND BSUM
!---------------------------------------------------------------------*/
      SUMBR = ZEROR;
      SUMBI = ZEROI;
      for (K=1; K<=KMAX; K++) {
        SUMBR = SUMBR + PR[K]*BETA[K];
        SUMBI = SUMBI + XPI[K]*BETA[K];
      }
      *ASUMR = ZEROR;
      *ASUMI = ZEROI;
      *BSUMR = SUMBR;
      *BSUMI = SUMBI;
      L1 = 0;
      L2 = 30;
      BTOL = TOL*(ABS(*BSUMR)+ABS(*BSUMI));
      ATOL = TOL;
      PP = 1.0;
      IAS = 0;
      IBS = 0;
      if (RFNU2 < TOL) goto e110;
      for (IS=2; IS<8; IS++) {
        ATOL = ATOL/RFNU2;
        PP = PP*RFNU2;
        if (IAS == 1) goto e60;
        SUMAR = ZEROR;
        SUMAI = ZEROI;
        for (K=1; K<=KMAX; K++) {
          M = L1 + K;
          SUMAR = SUMAR + PR[K]*ALFA[M];
          SUMAI = SUMAI + XPI[K]*ALFA[M];
          if (AP[K] < ATOL) goto e50;
        }
e50:    *ASUMR = *ASUMR + SUMAR*PP;
        *ASUMI = *ASUMI + SUMAI*PP;