1 SUBROUTINE cunhj(Z, FNU, IPMTR, TOL, PHI, ARG, ZETA1, ZETA2,
34 COMPLEX ARG, ASUM, BSUM, CFNU, CONE, CR, CZERO, DR, P, PHI,
35 * przth, ptfn, rfn13, rtzta, rzth, suma, sumb, tfn, t2, up, w, w2,
36 * z, za, zb, zc, zeta, zeta1, zeta2, zth
37 REAL ALFA, ANG, AP, AR, ATOL, AW2, AZTH, BETA, BR, BTOL, C, EX1,
38 * ex2, fnu, fn13, fn23, gama, hpi, pi, pp, rfnu, rfnu2, thpi, tol,
39 *
wi, wr, zci, zcr, zetai, zetar, zthi, zthr, asumr, asumi, bsumr,
40 * bsumi,
test, tstr, tsti, ac
41 INTEGER IAS, IBS, IPMTR, IS, J, JR, JU, K, KMAX, KP1, KS, L, LR,
43 dimension ar(14), br(14), c(105), alfa(180), beta(210), gama(30),
44 * ap(30), p(30), up(14), cr(14), dr(14)
45 DATA ar(1), ar(2), ar(3), ar(4), ar(5), ar(6), ar(7), ar(8),
46 1 ar(9), ar(10), ar(11), ar(12), ar(13), ar(14)/
47 2 1.00000000000000000
e+00, 1.04166666666666667
e-01,
48 3 8.35503472222222222
e-02, 1.28226574556327160
e-01,
49 4 2.91849026464140464
e-01, 8.81627267443757652
e-01,
50 5 3.32140828186276754
e+00, 1.49957629868625547
e+01,
51 6 7.89230130115865181
e+01, 4.74451538868264323
e+02,
52 7 3.20749009089066193
e+03, 2.40865496408740049
e+04,
53 8 1.98923119169509794
e+05, 1.79190200777534383
e+06/
54 DATA br(1), br(2), br(3), br(4), br(5), br(6), br(7), br(8),
55 1 br(9), br(10), br(11), br(12), br(13), br(14)/
56 2 1.00000000000000000
e+00, -1.45833333333333333
e-01,
57 3 -9.87413194444444444
e-02, -1.43312053915895062
e-01,
58 4 -3.17227202678413548
e-01, -9.42429147957120249
e-01,
59 5 -3.51120304082635426
e+00, -1.57272636203680451
e+01,
60 6 -8.22814390971859444
e+01, -4.92355370523670524
e+02,
61 7 -3.31621856854797251
e+03, -2.48276742452085896
e+04,
62 8 -2.04526587315129788
e+05, -1.83844491706820990
e+06/
63 DATA c(1), c(2), c(3), c(4), c(5), c(6), c(7), c(8), c(9), c(10),
64 1 c(11), c(12), c(13), c(14), c(15), c(16), c(17), c(18),
65 2 c(19), c(20), c(21), c(22), c(23), c(24)/
66 3 1.00000000000000000
e+00, -2.08333333333333333
e-01,
67 4 1.25000000000000000
e-01, 3.34201388888888889
e-01,
68 5 -4.01041666666666667
e-01, 7.03125000000000000
e-02,
69 6 -1.02581259645061728
e+00, 1.84646267361111111
e+00,
70 7 -8.91210937500000000
e-01, 7.32421875000000000
e-02,
71 8 4.66958442342624743
e+00, -1.12070026162229938
e+01,
72 9 8.78912353515625000
e+00, -2.36408691406250000
e+00,
73 a 1.12152099609375000
e-01, -2.82120725582002449
e+01,
74 b 8.46362176746007346
e+01, -9.18182415432400174
e+01,
75 c 4.25349987453884549
e+01, -7.36879435947963170
e+00,
76 d 2.27108001708984375
e-01, 2.12570130039217123
e+02,
77 e -7.65252468141181642
e+02, 1.05999045252799988
e+03/
78 DATA c(25), c(26), c(27), c(28), c(29), c(30), c(31), c(32),
79 1 c(33), c(34), c(35), c(36), c(37), c(38), c(39), c(40),
80 2 c(41), c(42), c(43), c(44), c(45), c(46), c(47), c(48)/
81 3 -6.99579627376132541
e+02, 2.18190511744211590
e+02,
82 4 -2.64914304869515555
e+01, 5.72501420974731445
e-01,
83 5 -1.91945766231840700
e+03, 8.06172218173730938
e+03,
84 6 -1.35865500064341374
e+04, 1.16553933368645332
e+04,
85 7 -5.30564697861340311
e+03, 1.20090291321635246
e+03,
86 8 -1.08090919788394656
e+02, 1.72772750258445740
e+00,
87 9 2.02042913309661486
e+04, -9.69805983886375135
e+04,
88 a 1.92547001232531532
e+05, -2.03400177280415534
e+05,
89 b 1.22200464983017460
e+05, -4.11926549688975513
e+04,
90 c 7.10951430248936372
e+03, -4.93915304773088012
e+02,
91 d 6.07404200127348304
e+00, -2.42919187900551333
e+05,
92 e 1.31176361466297720
e+06, -2.99801591853810675
e+06/
93 DATA c(49), c(50), c(51), c(52), c(53), c(54), c(55), c(56),
94 1 c(57), c(58), c(59), c(60), c(61), c(62), c(63), c(64),
95 2 c(65), c(66), c(67), c(68), c(69), c(70), c(71), c(72)/
96 3 3.76327129765640400
e+06, -2.81356322658653411
e+06,
97 4 1.26836527332162478
e+06, -3.31645172484563578
e+05,
98 5 4.52187689813627263
e+04, -2.49983048181120962
e+03,
99 6 2.43805296995560639
e+01, 3.28446985307203782
e+06,
100 7 -1.97068191184322269
e+07, 5.09526024926646422
e+07,
101 8 -7.41051482115326577
e+07, 6.63445122747290267
e+07,
102 9 -3.75671766607633513
e+07, 1.32887671664218183
e+07,
103 a -2.78561812808645469
e+06, 3.08186404612662398
e+05,
104 b -1.38860897537170405
e+04, 1.10017140269246738
e+02,
105 c -4.93292536645099620
e+07, 3.25573074185765749
e+08,
106 d -9.39462359681578403
e+08, 1.55359689957058006
e+09,
107 e -1.62108055210833708
e+09, 1.10684281682301447
e+09/
108 DATA c(73), c(74), c(75), c(76), c(77), c(78), c(79), c(80),
109 1 c(81), c(82), c(83), c(84), c(85), c(86), c(87), c(88),
110 2 c(89), c(90), c(91), c(92), c(93), c(94), c(95), c(96)/
111 3 -4.95889784275030309
e+08, 1.42062907797533095
e+08,
112 4 -2.44740627257387285
e+07, 2.24376817792244943
e+06,
113 5 -8.40054336030240853
e+04, 5.51335896122020586
e+02,
114 6 8.14789096118312115
e+08, -5.86648149205184723
e+09,
115 7 1.86882075092958249
e+10, -3.46320433881587779
e+10,
116 8 4.12801855797539740
e+10, -3.30265997498007231
e+10,
117 9 1.79542137311556001
e+10, -6.56329379261928433
e+09,
118 a 1.55927986487925751
e+09, -2.25105661889415278
e+08,
119 b 1.73951075539781645
e+07, -5.49842327572288687
e+05,
120 c 3.03809051092238427
e+03, -1.46792612476956167
e+10,
121 d 1.14498237732025810
e+11, -3.99096175224466498
e+11,
122 e 8.19218669548577329
e+11, -1.09837515608122331
e+12/
123 DATA c(97), c(98), c(99), c(100), c(101), c(102), c(103), c(104),
125 2 1.00815810686538209
e+12, -6.45364869245376503
e+11,
126 3 2.87900649906150589
e+11, -8.78670721780232657
e+10,
127 4 1.76347306068349694
e+10, -2.16716498322379509
e+09,
128 5 1.43157876718888981
e+08, -3.87183344257261262
e+06,
129 6 1.82577554742931747
e+04/
130 DATA alfa(1), alfa(2), alfa(3), alfa(4), alfa(5), alfa(6),
131 1 alfa(7), alfa(8), alfa(9), alfa(10), alfa(11), alfa(12),
132 2 alfa(13), alfa(14), alfa(15), alfa(16), alfa(17), alfa(18),
133 3 alfa(19), alfa(20), alfa(21), alfa(22)/
134 4 -4.44444444444444444
e-03, -9.22077922077922078
e-04,
135 5 -8.84892884892884893
e-05, 1.65927687832449737
e-04,
136 6 2.46691372741792910
e-04, 2.65995589346254780
e-04,
137 7 2.61824297061500945
e-04, 2.48730437344655609
e-04,
138 8 2.32721040083232098
e-04, 2.16362485712365082
e-04,
139 9 2.00738858762752355
e-04, 1.86267636637545172
e-04,
140 a 1.73060775917876493
e-04, 1.61091705929015752
e-04,
141 b 1.50274774160908134
e-04, 1.40503497391269794
e-04,
142 c 1.31668816545922806
e-04, 1.23667445598253261
e-04,
143 d 1.16405271474737902
e-04, 1.09798298372713369
e-04,
144 e 1.03772410422992823
e-04, 9.82626078369363448
e-05/
145 DATA alfa(23), alfa(24), alfa(25), alfa(26), alfa(27), alfa(28),
146 1 alfa(29), alfa(30), alfa(31), alfa(32), alfa(33), alfa(34),
147 2 alfa(35), alfa(36), alfa(37), alfa(38), alfa(39), alfa(40),
148 3 alfa(41), alfa(42), alfa(43), alfa(44)/
149 4 9.32120517249503256
e-05, 8.85710852478711718
e-05,
150 5 8.42963105715700223
e-05, 8.03497548407791151
e-05,
151 6 7.66981345359207388
e-05, 7.33122157481777809
e-05,
152 7 7.01662625163141333
e-05, 6.72375633790160292
e-05,
153 8 6.93735541354588974
e-04, 2.32241745182921654
e-04,
154 9 -1.41986273556691197
e-05, -1.16444931672048640
e-04,
155 a -1.50803558053048762
e-04, -1.55121924918096223
e-04,
156 b -1.46809756646465549
e-04, -1.33815503867491367
e-04,
157 c -1.19744975684254051
e-04, -1.06184319207974020
e-04,
158 d -9.37699549891194492
e-05, -8.26923045588193274
e-05,
159 e -7.29374348155221211
e-05, -6.44042357721016283
e-05/
160 DATA alfa(45), alfa(46), alfa(47), alfa(48), alfa(49), alfa(50),
161 1 alfa(51), alfa(52), alfa(53), alfa(54), alfa(55), alfa(56),
162 2 alfa(57), alfa(58), alfa(59), alfa(60), alfa(61), alfa(62),
163 3 alfa(63), alfa(64), alfa(65), alfa(66)/
164 4 -5.69611566009369048
e-05, -5.04731044303561628
e-05,
165 5 -4.48134868008882786
e-05, -3.98688727717598864
e-05,
166 6 -3.55400532972042498
e-05, -3.17414256609022480
e-05,
167 7 -2.83996793904174811
e-05, -2.54522720634870566
e-05,
168 8 -2.28459297164724555
e-05, -2.05352753106480604
e-05,
169 9 -1.84816217627666085
e-05, -1.66519330021393806
e-05,
170 a -1.50179412980119482
e-05, -1.35554031379040526
e-05,
171 b -1.22434746473858131
e-05, -1.10641884811308169
e-05,
172 c -3.54211971457743841
e-04, -1.56161263945159416
e-04,
173 d 3.04465503594936410
e-05, 1.30198655773242693
e-04,
174 e 1.67471106699712269
e-04, 1.70222587683592569
e-04/
175 DATA alfa(67), alfa(68), alfa(69), alfa(70), alfa(71), alfa(72),
176 1 alfa(73), alfa(74), alfa(75), alfa(76), alfa(77), alfa(78),
177 2 alfa(79), alfa(80), alfa(81), alfa(82), alfa(83), alfa(84),
178 3 alfa(85), alfa(86), alfa(87), alfa(88)/
179 4 1.56501427608594704
e-04, 1.36339170977445120
e-04,
180 5 1.14886692029825128
e-04, 9.45869093034688111
e-05,
181 6 7.64498419250898258
e-05, 6.07570334965197354
e-05,
182 7 4.74394299290508799
e-05, 3.62757512005344297
e-05,
183 8 2.69939714979224901
e-05, 1.93210938247939253
e-05,
184 9 1.30056674793963203
e-05, 7.82620866744496661
e-06,
185 a 3.59257485819351583
e-06, 1.44040049814251817
e-07,
186 b -2.65396769697939116
e-06, -4.91346867098485910
e-06,
187 c -6.72739296091248287
e-06, -8.17269379678657923
e-06,
188 d -9.31304715093561232
e-06, -1.02011418798016441
e-05,
189 e -1.08805962510592880
e-05, -1.13875481509603555
e-05/
190 DATA alfa(89), alfa(90), alfa(91), alfa(92), alfa(93), alfa(94),
191 1 alfa(95), alfa(96), alfa(97), alfa(98), alfa(99), alfa(100),
192 2 alfa(101), alfa(102), alfa(103), alfa(104), alfa(105),
193 3 alfa(106), alfa(107), alfa(108), alfa(109), alfa(110)/
194 4 -1.17519675674556414
e-05, -1.19987364870944141
e-05,
195 5 3.78194199201772914
e-04, 2.02471952761816167
e-04,
196 6 -6.37938506318862408
e-05, -2.38598230603005903
e-04,
197 7 -3.10916256027361568
e-04, -3.13680115247576316
e-04,
198 8 -2.78950273791323387
e-04, -2.28564082619141374
e-04,
199 9 -1.75245280340846749
e-04, -1.25544063060690348
e-04,
200 a -8.22982872820208365
e-05, -4.62860730588116458
e-05,
201 b -1.72334302366962267
e-05, 5.60690482304602267
e-06,
202 c 2.31395443148286800
e-05, 3.62642745856793957
e-05,
203 d 4.58006124490188752
e-05, 5.24595294959114050
e-05,
204 e 5.68396208545815266
e-05, 5.94349820393104052
e-05/
205 DATA alfa(111), alfa(112), alfa(113), alfa(114), alfa(115),
206 1 alfa(116), alfa(117), alfa(118), alfa(119), alfa(120),
207 2 alfa(121), alfa(122), alfa(123), alfa(124), alfa(125),
208 3 alfa(126), alfa(127), alfa(128), alfa(129), alfa(130)/
209 4 6.06478527578421742
e-05, 6.08023907788436497
e-05,
210 5 6.01577894539460388
e-05, 5.89199657344698500
e-05,
211 6 5.72515823777593053
e-05, 5.52804375585852577
e-05,
212 7 5.31063773802880170
e-05, 5.08069302012325706
e-05,
213 8 4.84418647620094842
e-05, 4.60568581607475370
e-05,
214 9 -6.91141397288294174
e-04, -4.29976633058871912
e-04,
215 a 1.83067735980039018
e-04, 6.60088147542014144
e-04,
216 b 8.75964969951185931
e-04, 8.77335235958235514
e-04,
217 c 7.49369585378990637
e-04, 5.63832329756980918
e-04,
218 d 3.68059319971443156
e-04, 1.88464535514455599
e-04/
219 DATA alfa(131), alfa(132), alfa(133), alfa(134), alfa(135),
220 1 alfa(136), alfa(137), alfa(138), alfa(139), alfa(140),
221 2 alfa(141), alfa(142), alfa(143), alfa(144), alfa(145),
222 3 alfa(146), alfa(147), alfa(148), alfa(149), alfa(150)/
223 4 3.70663057664904149
e-05, -8.28520220232137023
e-05,
224 5 -1.72751952869172998
e-04, -2.36314873605872983
e-04,
225 6 -2.77966150694906658
e-04, -3.02079514155456919
e-04,
226 7 -3.12594712643820127
e-04, -3.12872558758067163
e-04,
227 8 -3.05678038466324377
e-04, -2.93226470614557331
e-04,
228 9 -2.77255655582934777
e-04, -2.59103928467031709
e-04,
229 a -2.39784014396480342
e-04, -2.20048260045422848
e-04,
230 b -2.00443911094971498
e-04, -1.81358692210970687
e-04,
231 c -1.63057674478657464
e-04, -1.45712672175205844
e-04,
232 d -1.29425421983924587
e-04, -1.14245691942445952
e-04/
233 DATA alfa(151), alfa(152), alfa(153), alfa(154), alfa(155),
234 1 alfa(156), alfa(157), alfa(158), alfa(159), alfa(160),
235 2 alfa(161), alfa(162), alfa(163), alfa(164), alfa(165),
236 3 alfa(166), alfa(167), alfa(168), alfa(169), alfa(170)/
237 4 1.92821964248775885
e-03, 1.35592576302022234
e-03,
238 5 -7.17858090421302995
e-04, -2.58084802575270346
e-03,
239 6 -3.49271130826168475
e-03, -3.46986299340960628
e-03,
240 7 -2.82285233351310182
e-03, -1.88103076404891354
e-03,
241 8 -8.89531718383947600
e-04, 3.87912102631035228
e-06,
242 9 7.28688540119691412
e-04, 1.26566373053457758
e-03,
243 a 1.62518158372674427
e-03, 1.83203153216373172
e-03,
244 b 1.91588388990527909
e-03, 1.90588846755546138
e-03,
245 c 1.82798982421825727
e-03, 1.70389506421121530
e-03,
246 d 1.55097127171097686
e-03, 1.38261421852276159
e-03/
247 DATA alfa(171), alfa(172), alfa(173), alfa(174), alfa(175),
248 1 alfa(176), alfa(177), alfa(178), alfa(179), alfa(180)/
249 2 1.20881424230064774
e-03, 1.03676532638344962
e-03,
250 3 8.71437918068619115
e-04, 7.16080155297701002
e-04,
251 4 5.72637002558129372
e-04, 4.42089819465802277
e-04,
252 5 3.24724948503090564
e-04, 2.20342042730246599
e-04,
253 6 1.28412898401353882
e-04, 4.82005924552095464
e-05/
254 DATA beta(1), beta(2), beta(3), beta(4), beta(5), beta(6),
255 1 beta(7), beta(8), beta(9), beta(10), beta(11), beta(12),
256 2 beta(13), beta(14), beta(15), beta(16), beta(17), beta(18),
257 3 beta(19), beta(20), beta(21), beta(22)/
258 4 1.79988721413553309
e-02, 5.59964911064388073
e-03,
259 5 2.88501402231132779
e-03, 1.80096606761053941
e-03,
260 6 1.24753110589199202
e-03, 9.22878876572938311
e-04,
261 7 7.14430421727287357
e-04, 5.71787281789704872
e-04,
262 8 4.69431007606481533
e-04, 3.93232835462916638
e-04,
263 9 3.34818889318297664
e-04, 2.88952148495751517
e-04,
264 a 2.52211615549573284
e-04, 2.22280580798883327
e-04,
265 b 1.97541838033062524
e-04, 1.76836855019718004
e-04,
266 c 1.59316899661821081
e-04, 1.44347930197333986
e-04,
267 d 1.31448068119965379
e-04, 1.20245444949302884
e-04,
268 e 1.10449144504599392
e-04, 1.01828770740567258
e-04/
269 DATA beta(23), beta(24), beta(25), beta(26), beta(27), beta(28),
270 1 beta(29), beta(30), beta(31), beta(32), beta(33), beta(34),
271 2 beta(35), beta(36), beta(37), beta(38), beta(39), beta(40),
272 3 beta(41), beta(42), beta(43), beta(44)/
273 4 9.41998224204237509
e-05, 8.74130545753834437
e-05,
274 5 8.13466262162801467
e-05, 7.59002269646219339
e-05,
275 6 7.09906300634153481
e-05, 6.65482874842468183
e-05,
276 7 6.25146958969275078
e-05, 5.88403394426251749
e-05,
277 8 -1.49282953213429172
e-03, -8.78204709546389328
e-04,
278 9 -5.02916549572034614
e-04, -2.94822138512746025
e-04,
279 a -1.75463996970782828
e-04, -1.04008550460816434
e-04,
280 b -5.96141953046457895
e-05, -3.12038929076098340
e-05,
281 c -1.26089735980230047
e-05, -2.42892608575730389
e-07,
282 d 8.05996165414273571
e-06, 1.36507009262147391
e-05,
283 e 1.73964125472926261
e-05, 1.98672978842133780
e-05/
284 DATA beta(45), beta(46), beta(47), beta(48), beta(49), beta(50),
285 1 beta(51), beta(52), beta(53), beta(54), beta(55), beta(56),
286 2 beta(57), beta(58), beta(59), beta(60), beta(61), beta(62),
287 3 beta(63), beta(64), beta(65), beta(66)/
288 4 2.14463263790822639
e-05, 2.23954659232456514
e-05,
289 5 2.28967783814712629
e-05, 2.30785389811177817
e-05,
290 6 2.30321976080909144
e-05, 2.28236073720348722
e-05,
291 7 2.25005881105292418
e-05, 2.20981015361991429
e-05,
292 8 2.16418427448103905
e-05, 2.11507649256220843
e-05,
293 9 2.06388749782170737
e-05, 2.01165241997081666
e-05,
294 a 1.95913450141179244
e-05, 1.90689367910436740
e-05,
295 b 1.85533719641636667
e-05, 1.80475722259674218
e-05,
296 c 5.52213076721292790
e-04, 4.47932581552384646
e-04,
297 d 2.79520653992020589
e-04, 1.52468156198446602
e-04,
298 e 6.93271105657043598
e-05, 1.76258683069991397
e-05/
299 DATA beta(67), beta(68), beta(69), beta(70), beta(71), beta(72),
300 1 beta(73), beta(74), beta(75), beta(76), beta(77), beta(78),
301 2 beta(79), beta(80), beta(81), beta(82), beta(83), beta(84),
302 3 beta(85), beta(86), beta(87), beta(88)/
303 4 -1.35744996343269136
e-05, -3.17972413350427135
e-05,
304 5 -4.18861861696693365
e-05, -4.69004889379141029
e-05,
305 6 -4.87665447413787352
e-05, -4.87010031186735069
e-05,
306 7 -4.74755620890086638
e-05, -4.55813058138628452
e-05,
307 8 -4.33309644511266036
e-05, -4.09230193157750364
e-05,
308 9 -3.84822638603221274
e-05, -3.60857167535410501
e-05,
309 a -3.37793306123367417
e-05, -3.15888560772109621
e-05,
310 b -2.95269561750807315
e-05, -2.75978914828335759
e-05,
311 c -2.58006174666883713
e-05, -2.41308356761280200
e-05,
312 d -2.25823509518346033
e-05, -2.11479656768912971
e-05,
313 e -1.98200638885294927
e-05, -1.85909870801065077
e-05/
314 DATA beta(89), beta(90), beta(91), beta(92), beta(93), beta(94),
315 1 beta(95), beta(96), beta(97), beta(98), beta(99), beta(100),
316 2 beta(101), beta(102), beta(103), beta(104), beta(105),
317 3 beta(106), beta(107), beta(108), beta(109), beta(110)/
318 4 -1.74532699844210224
e-05, -1.63997823854497997
e-05,
319 5 -4.74617796559959808
e-04, -4.77864567147321487
e-04,
320 6 -3.20390228067037603
e-04, -1.61105016119962282
e-04,
321 7 -4.25778101285435204
e-05, 3.44571294294967503
e-05,
322 8 7.97092684075674924
e-05, 1.03138236708272200
e-04,
323 9 1.12466775262204158
e-04, 1.13103642108481389
e-04,
324 a 1.08651634848774268
e-04, 1.01437951597661973
e-04,
325 b 9.29298396593363896
e-05, 8.40293133016089978
e-05,
326 c 7.52727991349134062
e-05, 6.69632521975730872
e-05,
327 d 5.92564547323194704
e-05, 5.22169308826975567
e-05,
328 e 4.58539485165360646
e-05, 4.01445513891486808
e-05/
329 DATA beta(111), beta(112), beta(113), beta(114), beta(115),
330 1 beta(116), beta(117), beta(118), beta(119), beta(120),
331 2 beta(121), beta(122), beta(123), beta(124), beta(125),
332 3 beta(126), beta(127), beta(128), beta(129), beta(130)/
333 4 3.50481730031328081
e-05, 3.05157995034346659
e-05,
334 5 2.64956119950516039
e-05, 2.29363633690998152
e-05,
335 6 1.97893056664021636
e-05, 1.70091984636412623
e-05,
336 7 1.45547428261524004
e-05, 1.23886640995878413
e-05,
337 8 1.04775876076583236
e-05, 8.79179954978479373
e-06,
338 9 7.36465810572578444
e-04, 8.72790805146193976
e-04,
339 a 6.22614862573135066
e-04, 2.85998154194304147
e-04,
340 b 3.84737672879366102
e-06, -1.87906003636971558
e-04,
341 c -2.97603646594554535
e-04, -3.45998126832656348
e-04,
342 d -3.53382470916037712
e-04, -3.35715635775048757
e-04/
343 DATA beta(131), beta(132), beta(133), beta(134), beta(135),
344 1 beta(136), beta(137), beta(138), beta(139), beta(140),
345 2 beta(141), beta(142), beta(143), beta(144), beta(145),
346 3 beta(146), beta(147), beta(148), beta(149), beta(150)/
347 4 -3.04321124789039809
e-04, -2.66722723047612821
e-04,
348 5 -2.27654214122819527
e-04, -1.89922611854562356
e-04,
349 6 -1.55058918599093870
e-04, -1.23778240761873630
e-04,
350 7 -9.62926147717644187
e-05, -7.25178327714425337
e-05,
351 8 -5.22070028895633801
e-05, -3.50347750511900522
e-05,
352 9 -2.06489761035551757
e-05, -8.70106096849767054
e-06,
353 a 1.13698686675100290
e-06, 9.16426474122778849
e-06,
354 b 1.56477785428872620
e-05, 2.08223629482466847
e-05,
355 c 2.48923381004595156
e-05, 2.80340509574146325
e-05,
356 d 3.03987774629861915
e-05, 3.21156731406700616
e-05/
357 DATA beta(151), beta(152), beta(153), beta(154), beta(155),
358 1 beta(156), beta(157), beta(158), beta(159), beta(160),
359 2 beta(161), beta(162), beta(163), beta(164), beta(165),
360 3 beta(166), beta(167), beta(168), beta(169), beta(170)/
361 4 -1.80182191963885708
e-03, -2.43402962938042533
e-03,
362 5 -1.83422663549856802
e-03, -7.62204596354009765
e-04,
363 6 2.39079475256927218
e-04, 9.49266117176881141
e-04,
364 7 1.34467449701540359
e-03, 1.48457495259449178
e-03,
365 8 1.44732339830617591
e-03, 1.30268261285657186
e-03,
366 9 1.10351597375642682
e-03, 8.86047440419791759
e-04,
367 a 6.73073208165665473
e-04, 4.77603872856582378
e-04,
368 b 3.05991926358789362
e-04, 1.60315694594721630
e-04,
369 c 4.00749555270613286
e-05, -5.66607461635251611
e-05,
370 d -1.32506186772982638
e-04, -1.90296187989614057
e-04/
371 DATA beta(171), beta(172), beta(173), beta(174), beta(175),
372 1 beta(176), beta(177), beta(178), beta(179), beta(180),
373 2 beta(181), beta(182), beta(183), beta(184), beta(185),
374 3 beta(186), beta(187), beta(188), beta(189), beta(190)/
375 4 -2.32811450376937408
e-04, -2.62628811464668841
e-04,
376 5 -2.82050469867598672
e-04, -2.93081563192861167
e-04,
377 6 -2.97435962176316616
e-04, -2.96557334239348078
e-04,
378 7 -2.91647363312090861
e-04, -2.83696203837734166
e-04,
379 8 -2.73512317095673346
e-04, -2.61750155806768580
e-04,
380 9 6.38585891212050914
e-03, 9.62374215806377941
e-03,
381 a 7.61878061207001043
e-03, 2.83219055545628054
e-03,
382 b -2.09841352012720090
e-03, -5.73826764216626498
e-03,
383 c -7.70804244495414620
e-03, -8.21011692264844401
e-03,
384 d -7.65824520346905413
e-03, -6.47209729391045177
e-03/
385 DATA beta(191), beta(192), beta(193), beta(194), beta(195),
386 1 beta(196), beta(197), beta(198), beta(199), beta(200),
387 2 beta(201), beta(202), beta(203), beta(204), beta(205),
388 3 beta(206), beta(207), beta(208), beta(209), beta(210)/
389 4 -4.99132412004966473
e-03, -3.45612289713133280
e-03,
390 5 -2.01785580014170775
e-03, -7.59430686781961401
e-04,
391 6 2.84173631523859138
e-04, 1.10891667586337403
e-03,
392 7 1.72901493872728771
e-03, 2.16812590802684701
e-03,
393 8 2.45357710494539735
e-03, 2.61281821058334862
e-03,
394 9 2.67141039656276912
e-03, 2.65203073395980430
e-03,
395 a 2.57411652877287315
e-03, 2.45389126236094427
e-03,
396 b 2.30460058071795494
e-03, 2.13684837686712662
e-03,
397 c 1.95896528478870911
e-03, 1.77737008679454412
e-03,
398 d 1.59690280765839059
e-03, 1.42111975664438546
e-03/
399 DATA gama(1), gama(2), gama(3), gama(4), gama(5), gama(6),
400 1 gama(7), gama(8), gama(9), gama(10), gama(11), gama(12),
401 2 gama(13), gama(14), gama(15), gama(16), gama(17), gama(18),
402 3 gama(19), gama(20), gama(21), gama(22)/
403 4 6.29960524947436582
e-01, 2.51984209978974633
e-01,
404 5 1.54790300415655846
e-01, 1.10713062416159013
e-01,
405 6 8.57309395527394825
e-02, 6.97161316958684292
e-02,
406 7 5.86085671893713576
e-02, 5.04698873536310685
e-02,
407 8 4.42600580689154809
e-02, 3.93720661543509966
e-02,
408 9 3.54283195924455368
e-02, 3.21818857502098231
e-02,
409 a 2.94646240791157679
e-02, 2.71581677112934479
e-02,
410 b 2.51768272973861779
e-02, 2.34570755306078891
e-02,
411 c 2.19508390134907203
e-02, 2.06210828235646240
e-02,
412 d 1.94388240897880846
e-02, 1.83810633800683158
e-02,
413 e 1.74293213231963172
e-02, 1.65685837786612353
e-02/
414 DATA gama(23), gama(24), gama(25), gama(26), gama(27), gama(28),
415 1 gama(29), gama(30)/
416 2 1.57865285987918445
e-02, 1.50729501494095594
e-02,
417 3 1.44193250839954639
e-02, 1.38184805735341786
e-02,
418 4 1.32643378994276568
e-02, 1.27517121970498651
e-02,
419 5 1.22761545318762767
e-02, 1.18338262398482403
e-02/
420 DATA ex1, ex2, hpi, pi, thpi /
421 1 3.33333333333333333
e-01, 6.66666666666666667
e-01,
422 2 1.57079632679489662
e+00, 3.14159265358979324
e+00,
423 3 4.71238898038468986
e+00/
424 DATA czero, cone / (0.0e0,0.0e0), (1.0e0,0.0e0) /
433 test = r1mach(1)*1.0
e+3
435 IF (
abs(tstr).GT.ac .OR.
abs(tsti).GT.ac) go
to 15
437 zeta1 =
cmplx(ac,0.0e0)
438 zeta2 =
cmplx(fnu,0.0e0)
443 zb = z*
cmplx(rfnu,0.0e0)
450 rfn13 =
cmplx(1.0e0/fn13,0.0e0)
453 IF (aw2.GT.0.25e0) go
to 130
459 suma =
cmplx(gama(1),0.0e0)
461 IF (aw2.LT.tol) go
to 20
464 suma = suma + p(k)*
cmplx(gama(k),0.0e0)
466 IF (ap(k).LT.tol) go
to 20
472 arg = zeta*
cmplx(fn23,0.0e0)
474 zeta2 = csqrt(w2)*
cmplx(fnu,0.0e0)
475 zeta1 = zeta2*(cone+zeta*za*
cmplx(ex2,0.0e0))
477 phi = csqrt(za)*rfn13
478 IF (ipmtr.EQ.1) go
to 120
484 sumb = sumb + p(k)*
cmplx(beta(k),0.0e0)
490 btol = tol*cabs(bsum)
495 IF (rfnu2.LT.tol) go
to 110
499 IF (ias.EQ.1) go
to 60
503 suma = suma + p(k)*
cmplx(alfa(m),0.0e0)
504 IF (ap(k).LT.atol) go
to 50
507 asum = asum + suma*
cmplx(pp,0.0e0)
508 IF (pp.LT.tol) ias = 1
510 IF (ibs.EQ.1) go
to 90
514 sumb = sumb + p(k)*
cmplx(beta(m),0.0e0)
515 IF (ap(k).LT.atol) go
to 80
518 bsum = bsum + sumb*
cmplx(pp,0.0e0)
519 IF (pp.LT.btol) ibs = 1
521 IF (ias.EQ.1 .AND. ibs.EQ.1) go
to 110
527 pp = rfnu*
REAL(rfn13)
528 bsum = bsum*
cmplx(pp,0.0e0)
538 IF (wr.LT.0.0e0) wr = 0.0e0
539 IF (
wi.LT.0.0e0)
wi = 0.0e0
545 IF (zci.LT.0.0e0) zci = 0.0e0
546 IF (zci.GT.hpi) zci = hpi
547 IF (zcr.LT.0.0e0) zcr = 0.0e0
549 zth = (zc-w)*
cmplx(1.5e0,0.0e0)
550 cfnu =
cmplx(fnu,0.0e0)
557 IF (zthr.GE.0.0e0 .AND. zthi.LT.0.0e0) go
to 140
559 IF (zthr.EQ.0.0e0) go
to 140
560 ang = atan(zthi/zthr)
561 IF (zthr.LT.0.0e0) ang = ang + pi
567 IF (zetai.LT.0.0e0) zetai = 0.0e0
568 zeta =
cmplx(zetar,zetai)
569 arg = zeta*
cmplx(fn23,0.0e0)
572 phi = csqrt(za+za)*rfn13
573 IF (ipmtr.EQ.1) go
to 120
574 tfn =
cmplx(rfnu,0.0e0)/w
575 rzth =
cmplx(rfnu,0.0e0)/zth
576 zc = rzth*
cmplx(ar(2),0.0e0)
578 up(2) = (t2*
cmplx(c(2),0.0e0)+
cmplx(c(3),0.0e0))*tfn
581 IF (rfnu.LT.tol) go
to 220
588 btol = tol*(
abs(bsumr)+
abs(bsumi))
604 za =
cmplx(c(l),0.0e0)
607 za = za*t2 +
cmplx(c(l),0.0e0)
611 cr(ks) = przth*
cmplx(br(ks+1),0.0e0)
613 dr(ks) = przth*
cmplx(ar(ks+2),0.0e0)
616 IF (ias.EQ.1) go
to 180
621 suma = suma + cr(jr)*up(ju)
627 IF (pp.LT.tol .AND.
test.LT.tol) ias = 1
629 IF (ibs.EQ.1) go
to 200
630 sumb = up(lr+2) + up(lrp1)*zc
634 sumb = sumb + dr(jr)*up(ju)
640 IF (pp.LT.btol .AND.
test.LT.tol) ibs = 1
642 IF (ias.EQ.1 .AND. ibs.EQ.1) go
to 220
646 bsum = -bsum*rfn13/rtzta
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
octave_value sin(void) const
may be zero for pure relative error test tem the relative tolerance must be greater than or equal to
OCTAVE_EXPORT octave_value_list etc The functions then dimension(columns)
bool test(F fcn) const
Generic any/all test functionality with arbitrary predicate.
subroutine cunhj(Z, FNU, IPMTR, TOL, PHI, ARG, ZETA1, ZETA2, ASUM, BSUM)
octave_value cos(void) const
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup abs(local error in x(i))<