version 1.1, 2003/12/16 13:39:53
|
version 1.3, 2005/04/01 16:51:48
|
|
|
subroutine hcf2r(qq,ww,f2,r) | subroutine hcf2r(qq,ww,f2,r) |
| |
c include 'model_liang.cmn' | c include 'model_liang.cmn' |
implicit real*8 (a-h,o-z) |
implicit none |
integer i,j,j1,i1,tmp | integer i,j,j1,i1,tmp |
real*8 qq,ww,w1,w2,ww1,ww2,w11,w22 |
real*8 qq,ww,w1,w2,ww1,ww2 !,w11,w22 |
real*8 r,f2,nu,mp,mp2,r1,f21,r2,f22 | real*8 r,f2,nu,mp,mp2,r1,f21,r2,f22 |
real*8 q4,q6,frac,xx,fac,rlog,q8 |
real*8 q4,q6,xx,fac,rlog,q8 !,frac |
real*8 rp(276), f2p(552),tmp1,tmp2 | real*8 rp(276), f2p(552),tmp1,tmp2 |
real*8 q,qth,dq,y,f2th,lin,binding | real*8 q,qth,dq,y,f2th,lin,binding |
real*8 scb,quad | real*8 scb,quad |
|
|
| |
| |
c write(*,*) qq,ww | c write(*,*) qq,ww |
if (ww.lt.1.15) then |
if (ww.lt.1.15d0) then |
w1 = 0.d0 | w1 = 0.d0 |
w2 = 0.d0 | w2 = 0.d0 |
else if (ww.gt.3.9d0) then | else if (ww.gt.3.9d0) then |
|
|
| |
| |
data rp/ | data rp/ |
> 0.33129, 1.32645, -0.79900, -3.11559, 2.5681, -1.82534, |
> 0.33129d0, 1.32645d0,-0.79900d0, -3.11559d0, 2.5681d0, -1.82534d0, |
> 0.28813, 0.90973, -0.54869, -3.76651, 3.4846, -2.02225, |
> 0.28813d0, 0.90973d0,-0.54869d0, -3.76651d0, 3.4846d0, -2.02225d0, |
> 0.56382, 0.43268, 0.08276, 23.37011, -25.4063, 13.71927, |
> 0.56382d0, 0.43268d0, 0.08276d0, 23.37011d0, -25.4063d0, 13.71927d0, |
> 0.67816, -0.09742, 1.91168, 51.09415, -61.0272, 17.29566, |
> 0.67816d0,-0.09742d0, 1.91168d0, 51.09415d0, -61.0272d0, 17.29566d0, |
> 0.64090, 0.51581, 0.28952, -69.99436, -98.0096, 133.50420, |
> 0.64090d0, 0.51581d0, 0.28952d0, -69.99436d0, -98.0096d0, 133.50420d0, |
> 0.46245, -0.02240, 0.31779, -3.26504, 77.9481, 14.49329, |
> 0.46245d0,-0.02240d0, 0.31779d0, -3.26504d0, 77.9481d0, 14.49329d0, |
> 0.39436, -0.02232, 0.34328, 24.69035, 50.0479, 13.91268, |
> 0.39436d0,-0.02232d0, 0.34328d0, 24.69035d0, 50.0479d0, 13.91268d0, |
> 0.39805, 0.03919, 0.07716, 3.75484, 10.3324, 8.08673, |
> 0.39805d0, 0.03919d0, 0.07716d0, 3.75484d0, 10.3324d0, 8.08673d0, |
> 0.38997, 0.00900, 0.15343, 16.36940, 4.9903, 7.44212, |
> 0.38997d0, 0.00900d0, 0.15343d0, 16.36940d0, 4.9903d0, 7.44212d0, |
> 0.29674, 0.02783, 0.06118, 27.36369, -17.4484, 4.32732, |
> 0.29674d0, 0.02783d0, 0.06118d0, 27.36369d0, -17.4484d0, 4.32732d0, |
> 0.35439, -0.01126, 0.23150, 53.46022, -18.2759, 7.18313, |
> 0.35439d0,-0.01126d0, 0.23150d0, 53.46022d0, -18.2759d0, 7.18313d0, |
> 0.31987, 0.01455, 0.23862, 81.55396, -74.7982, 5.96615, |
> 0.31987d0, 0.01455d0, 0.23862d0, 81.55396d0, -74.7982d0, 5.96615d0, |
> 0.31163, -0.00699, 0.12995, 51.23706, 18.1217, 4.82674, |
> 0.31163d0,-0.00699d0, 0.12995d0, 51.23706d0, 18.1217d0, 4.82674d0, |
> 0.32176, 0.12831, 0.00514, 68.39391, -73.3401, 3.86676, |
> 0.32176d0, 0.12831d0, 0.00514d0, 68.39391d0, -73.3401d0, 3.86676d0, |
> 0.32514, 0.20568, 0.02398, 31.46726, -34.1291, 3.42527, |
> 0.32514d0, 0.20568d0, 0.02398d0, 31.46726d0, -34.1291d0, 3.42527d0, |
> 0.48734, 0.32263, 0.07471, 42.71339, -49.2712, 7.19561, |
> 0.48734d0, 0.32263d0, 0.07471d0, 42.71339d0, -49.2712d0, 7.19561d0, |
> 0.38366, 0.16378, 0.00110, 26.58418, -29.0223, 3.80760, |
> 0.38366d0, 0.16378d0, 0.00110d0, 26.58418d0, -29.0223d0, 3.80760d0, |
> 0.52817, -0.37320, 10.22177, 44.87227, -53.3864, 15.03388, |
> 0.52817d0,-0.37320d0,10.22177d0, 44.87227d0, -53.3864d0, 15.03388d0, |
> 0.34971, 0.24432, 0.03942, 15.46841, -17.0274, 2.87509, |
> 0.34971d0, 0.24432d0, 0.03942d0, 15.46841d0, -17.0274d0, 2.87509d0, |
> 0.32309, 0.05971, 0.23399, 81.18824, -85.6957, 4.59427, |
> 0.32309d0, 0.05971d0, 0.23399d0, 81.18824d0, -85.6957d0, 4.59427d0, |
> 0.48408, 0.14606, -0.12127, 78.59016, -92.2077, 2.15488, |
> 0.48408d0, 0.14606d0,-0.12127d0, 78.59016d0, -92.2077d0, 2.15488d0, |
> 0.44745, 0.08725, -0.07341, 96.91505, -113.3328, 1.61255, |
> 0.44745d0, 0.08725d0,-0.07341d0, 96.91505d0, -113.3328d0, 1.61255d0, |
> 0.17621, 0.00267, 0.16446, 290.88310, -284.5832, 4.78794, |
> 0.17621d0, 0.00267d0, 0.16446d0, 290.88310d0, -284.5832d0, 4.78794d0, |
> 0.19235, 0.03305, 0.26250, 86.32166, -86.1215, 4.45001, |
> 0.19235d0, 0.03305d0, 0.26250d0, 86.32166d0, -86.1215d0, 4.45001d0, |
> 0.25966, 0.02309, 0.06332, 54.29687, -49.3171, 3.33319, |
> 0.25966d0, 0.02309d0, 0.06332d0, 54.29687d0, -49.3171d0, 3.33319d0, |
> 0.35243, 0.02105, -0.18868, 110.69520, -226.4427, 4.95000, |
> 0.35243d0, 0.02105d0,-0.18868d0, 110.69520d0, -226.4427d0, 4.95000d0, |
> 0.36534, -0.21314, 2.81021, -59.05300, -37.7424, 31.42299, |
> 0.36534d0,-0.21314d0, 2.81021d0, -59.05300d0, -37.7424d0, 31.42299d0, |
> 0.28897, -5.27934, 62.75448, 15.15152, 5.7736, 29.82406, |
> 0.28897d0,-5.27934d0,62.75448d0, 15.15152d0, 5.7736d0, 29.82406d0, |
> 0.31192, -0.03546, 0.09286, 19.06862, -33.4685, 1.56198, |
> 0.31192d0,-0.03546d0, 0.09286d0, 19.06862d0, -33.4685d0, 1.56198d0, |
> 0.35713, -0.20097, 2.05940, 436.08360, -517.4285, 14.75464, |
> 0.35713d0,-0.20097d0, 2.05940d0, 436.08360d0, -517.4285d0, 14.75464d0, |
> 0.32423, 0.05549, 0.04984, 74.70284, -85.5771, 3.61358, |
> 0.32423d0, 0.05549d0, 0.04984d0, 74.70284d0, -85.5771d0, 3.61358d0, |
> 0.26321, 0.02532, 0.04165, 291.81850, -330.3480, 3.42817, |
> 0.26321d0, 0.02532d0, 0.04165d0, 291.81850d0, -330.3480d0, 3.42817d0, |
> 0.19235, -0.00521, 0.30864, 262.23130, -271.3294, 4.52254, |
> 0.19235d0,-0.00521d0, 0.30864d0, 262.23130d0, -271.3294d0, 4.52254d0, |
> 0.17858, -0.00018, 0.28614, 195.28650, -195.6441, 3.92928, |
> 0.17858d0,-0.00018d0, 0.28614d0, 195.28650d0, -195.6441d0, 3.92928d0, |
> 0.14180, 0.00801, 0.21545, 131.05210, -121.0251, 3.20372, |
> 0.14180d0, 0.00801d0, 0.21545d0, 131.05210d0, -121.0251d0, 3.20372d0, |
> 0.19810, 0.00240, 0.29005, 122.32390, -121.0733, 3.60232, |
> 0.19810d0, 0.00240d0, 0.29005d0, 122.32390d0, -121.0733d0, 3.60232d0, |
> 0.24832, -0.00999, 0.46397, 81.22332, -85.5181, 3.40233, |
> 0.24832d0,-0.00999d0, 0.46397d0, 81.22332d0, -85.5181d0, 3.40233d0, |
> 0.31740, 0.12344, 0.07026, 53.01398, -64.9467, 3.08821, |
> 0.31740d0, 0.12344d0, 0.07026d0, 53.01398d0, -64.9467d0, 3.08821d0, |
> 0.32739, -0.06064, 0.79821, 104.36790, -124.3242, 4.74774, |
> 0.32739d0,-0.06064d0, 0.79821d0, 104.36790d0, -124.3242d0, 4.74774d0, |
> 0.34597, 0.11624, -0.04088, 50.29695, -64.8436, 2.35262, |
> 0.34597d0, 0.11624d0,-0.04088d0, 50.29695d0, -64.8436d0, 2.35262d0, |
> 0.33840, 0.11035, 0.03128, 39.27810, -51.3317, 2.38903, |
> 0.33840d0, 0.11035d0, 0.03128d0, 39.27810d0, -51.3317d0, 2.38903d0, |
> 0.38341, 0.12239, -0.12886, 59.12659, -80.0135, 1.94849, |
> 0.38341d0, 0.12239d0,-0.12886d0, 59.12659d0, -80.0135d0, 1.94849d0, |
> 0.34759, -0.06730, 0.99427, 87.77917, -116.0158, 4.31371, |
> 0.34759d0,-0.06730d0, 0.99427d0, 87.77917d0, -116.0158d0, 4.31371d0, |
> 0.35035, 0.13825, 0.05372, 37.00152, -51.0832, 2.35058, |
> 0.35035d0, 0.13825d0, 0.05372d0, 37.00152d0, -51.0832d0, 2.35058d0, |
> 0.31877, 0.09654, 0.12551, 70.29255, -93.9217, 2.72640, |
> 0.31877d0, 0.09654d0, 0.12551d0, 70.29255d0, -93.9217d0, 2.72640d0, |
> 0.32985, 0.25111, -0.24644, 11.54706, -17.2718, 1.00506/ |
> 0.32985d0, 0.25111d0,-0.24644d0, 11.54706d0, -17.2718d0, 1.00506d0/ |
| |
data f2p/ | data f2p/ |
> 0.00270, -0.02563, 0.06438, 0.13288, -0.18442, 59.57852, |
> 0.00270d0, -0.02563d0, 0.06438d0, 0.13288d0, -0.18442d0, 59.57852d0, |
> 126.68970,-156.52300, 259.6862, 45.71608, 51.94866,-337.30260, |
> 126.68970d0,-156.52300d0,259.6862d0, 45.71608d0, 51.94866d0,-337.30260d0, |
> -0.01266, 0.53825, -1.14310, -1.10755, 2.67523, -96.68052, |
> -0.01266d0, 0.53825d0, -1.14310d0, -1.10755d0, 2.67523d0, -96.68052d0, |
> -65.71029, 144.62170, -350.7377, 61.76398, 116.55460, 172.52990, |
> -65.71029d0, 144.62170d0,-350.7377d0, 61.76398d0, 116.55460d0, 172.52990d0, |
> -0.00820, 0.30865, -0.10750, -0.79825, 1.01081, -37.89112, |
> -0.00820d0, 0.30865d0, -0.10750d0, -0.79825d0, 1.01081d0, -37.89112d0, |
> -99.42540, 120.79090, -161.3490, -27.73087, -37.25240, 223.31690, |
> -99.42540d0, 120.79090d0,-161.3490d0, -27.73087d0, -37.25240d0, 223.31690d0, |
> -0.01626, 0.68263, -2.08330, 0.93416, 7.66115, -20.34949, |
> -0.01626d0, 0.68263d0, -2.08330d0, 0.93416d0, 7.66115d0, -20.34949d0, |
> -89.79055, 100.32260, -89.9365, -81.08472, 33.33122, 139.55320, |
> -89.79055d0, 100.32260d0, -89.9365d0, -81.08472d0, 33.33122d0, 139.55320d0, |
> 0.10085, -0.39180, 4.39375, 2.96299, -9.43805, 5.53277, |
> 0.10085d0, -0.39180d0, 4.39375d0, 2.96299d0, -9.43805d0, 5.53277d0, |
> -49.30489, 43.73740, 12.3122,-123.36040, 117.44570, -2.55460, |
> -49.30489d0, 43.73740d0, 12.3122d0,-123.36040d0, 117.44570d0, -2.55460d0, |
> 0.16801, 1.66346, -3.58292, -1.43295, 6.54949, 9.59234, |
> 0.16801d0, 1.66346d0, -3.58292d0, -1.43295d0, 6.54949d0, 9.59234d0, |
> -15.43190, 7.56539, 32.3811, -79.52365, 79.85619, -30.67446, |
> -15.43190d0, 7.56539d0, 32.3811d0, -79.52365d0, 79.85619d0, -30.67446d0, |
> 0.24036, 0.66839, -1.02483, 2.73454, -3.12791, 4.67220, |
> 0.24036d0, 0.66839d0, -1.02483d0, 2.73454d0, -3.12791d0, 4.67220d0, |
> -1.94161, -1.02171, 14.8272, -21.21147, 20.23464, -11.97409, |
> -1.94161d0, -1.02171d0, 14.8272d0, -21.21147d0, 20.23464d0, -11.97409d0, |
> 0.14395, -0.24303, 0.47600, 5.20955, 4.07727, -6.50083, |
> 0.14395d0, -0.24303d0, 0.47600d0, 5.20955d0, 4.07727d0, -6.50083d0, |
> -59.81225, 63.16642, -30.5782, -70.20249, 50.86131, 55.80522, |
> -59.81225d0, 63.16642d0, -30.5782d0, -70.20249d0, 50.86131d0, 55.80522d0, |
> 0.19595, 0.49078, 0.09058, 2.08639, 0.61890, 3.84467, |
> 0.19595d0, 0.49078d0, 0.09058d0, 2.08639d0, 0.61890d0, 3.84467d0, |
> 30.61177, -30.33699, 17.4414, 32.96692, -20.06578, -29.88399, |
> 30.61177d0, -30.33699d0, 17.4414d0, 32.96692d0, -20.06578d0, -29.88399d0, |
> 0.16363, 0.21846, 0.20264, 2.53441, -0.76028, 3.90782, |
> 0.16363d0, 0.21846d0, 0.20264d0, 2.53441d0, -0.76028d0, 3.90782d0, |
> 36.58707, -35.68809, 18.2261, 40.25367, -24.84555, -33.50865, |
> 36.58707d0, -35.68809d0, 18.2261d0, 40.25367d0, -24.84555d0, -33.50865d0, |
> 0.11976, -0.45730, 1.08627, 2.63121, -0.58440, -5.72117, |
> 0.11976d0, -0.45730d0, 1.08627d0, 2.63121d0, -0.58440d0, -5.72117d0, |
> -53.80934, 57.73899, -25.7978, -52.77444, 37.38223, 48.12931, |
> -53.80934d0, 57.73899d0, -25.7978d0, -52.77444d0, 37.38223d0, 48.12931d0, |
> 0.03129, 0.08877, 0.31201, 0.52214, 0.56099, 9.45694, |
> 0.03129d0, 0.08877d0, 0.31201d0, 0.52214d0, 0.56099d0, 9.45694d0, |
> 80.42830, -77.88046, 41.7567, 69.77966, -24.70381, -89.21936, |
> 80.42830d0, -77.88046d0, 41.7567d0, 69.77966d0, -24.70381d0, -89.21936d0, |
> 0.04149, 0.05422, 0.14645, 0.46049, 0.70058, 7.52857, |
> 0.04149d0, 0.05422d0, 0.14645d0, 0.46049d0, 0.70058d0, 7.52857d0, |
> 73.36507, -69.99093, 33.8845, 68.44820, -34.26609, -69.93432, |
> 73.36507d0, -69.99093d0, 33.8845d0, 68.44820d0, -34.26609d0, -69.93432d0, |
> 0.03743, 0.05538, 0.06989, 0.42875, 0.74432, 9.91687, |
> 0.03743d0, 0.05538d0, 0.06989d0, 0.42875d0, 0.74432d0, 9.91687d0, |
> 99.10915, -95.21824, 44.3231, 88.44396, -44.56128, -92.21440, |
> 99.10915d0, -95.21824d0, 44.3231d0, 88.44396d0, -44.56128d0, -92.21440d0, |
> 0.04214, 0.07361, -0.05129, 0.40464, 1.03509, 8.78332, |
> 0.04214d0, 0.07361d0, -0.05129d0, 0.40464d0, 1.03509d0, 8.78332d0, |
> 90.67270, -86.53600, 39.0777, 79.55918, -42.40906, -79.76051, |
> 90.67270d0, -86.53600d0, 39.0777d0, 79.55918d0, -42.40906d0, -79.76051d0, |
> 0.02977, 0.24989, -0.36480, 0.43282, 1.66379, 10.58702, |
> 0.02977d0, 0.24989d0, -0.36480d0, 0.43282d0, 1.66379d0, 10.58702d0, |
> 99.35179, -95.66632, 45.5028, 74.12631, -26.76438, -97.86623, |
> 99.35179d0, -95.66632d0, 45.5028d0, 74.12631d0, -26.76438d0, -97.86623d0, |
> 0.04183, 0.17483, -0.23191, 0.45452, 1.41034, 8.17936, |
> 0.04183d0, 0.17483d0, -0.23191d0, 0.45452d0, 1.41034d0, 8.17936d0, |
> 74.28743, -69.96535, 34.5675, 50.35936, -11.35428, -76.36961, |
> 74.28743d0, -69.96535d0, 34.5675d0, 50.35936d0, -11.35428d0, -76.36961d0, |
> 0.05173, 0.04448, 0.13424, 0.38647, 0.54959, 6.16662, |
> 0.05173d0, 0.04448d0, 0.13424d0, 0.38647d0, 0.54959d0, 6.16662d0, |
> 51.08549, -46.36961, 25.3451, 28.23540, 3.70411, -58.01252, |
> 51.08549d0, -46.36961d0, 25.3451d0, 28.23540d0, 3.70411d0, -58.01252d0, |
> 0.08047, 0.07636, -0.05934, 0.34274, 0.77572, 4.27257, |
> 0.08047d0, 0.07636d0, -0.05934d0, 0.34274d0, 0.77572d0, 4.27257d0, |
> 38.90837, -35.54622, 17.7386, 24.76801, -9.17321, -33.32258, |
> 38.90837d0, -35.54622d0, 17.7386d0, 24.76801d0, -9.17321d0, -33.32258d0, |
> 0.08916, 0.08720, -0.04772, 0.36316, 0.77081, 4.10228, |
> 0.08916d0, 0.08720d0, -0.04772d0, 0.36316d0, 0.77081d0, 4.10228d0, |
> 40.26231, -37.09061, 17.1693, 26.59951, -14.07698, -29.96852, |
> 40.26231d0, -37.09061d0, 17.1693d0, 26.59951d0, -14.07698d0, -29.96852d0, |
> 0.02688, 0.20865, 0.24902, 0.19521, 0.37927, -0.57453, |
> 0.02688d0, 0.20865d0, 0.24902d0, 0.19521d0, 0.37927d0, -0.57453d0, |
> -9.00651, 14.38374, -2.4497, -6.51968, 9.73428, 3.89960, |
> -9.00651d0, 14.38374d0, -2.4497d0, -6.51968d0, 9.73428d0, 3.89960d0, |
> 0.04139, 0.15301, 0.32248, 0.36027, 0.17856, 3.32167, |
> 0.04139d0, 0.15301d0, 0.32248d0, 0.36027d0, 0.17856d0, 3.32167d0, |
> 28.88752, -23.90850, 13.4160, 12.82948, 3.90580, -29.60664, |
> 28.88752d0, -23.90850d0, 13.4160d0, 12.82948d0, 3.90580d0, -29.60664d0, |
> 0.05169, -0.19792, 0.70477, 0.98245, 0.19970, 1.07366, |
> 0.05169d0, -0.19792d0, 0.70477d0, 0.98245d0, 0.19970d0, 1.07366d0, |
> -0.96632, 13.46862, 3.5916, -10.17772, 41.48280, -29.53078, |
> -0.96632d0, 13.46862d0, 3.5916d0, -10.17772d0, 41.48280d0, -29.53078d0, |
> 0.04810, -0.14599, 0.44319, 0.86999, 0.59196, 2.37207, |
> 0.04810d0, -0.14599d0, 0.44319d0, 0.86999d0, 0.59196d0, 2.37207d0, |
> 16.49698, -3.58537, 9.1814, 2.29473, 31.25257, -39.15542, |
> 16.49698d0, -3.58537d0, 9.1814d0, 2.29473d0, 31.25257d0, -39.15542d0, |
> 0.04731, -0.14069, 0.33899, 0.74002, 0.62827, -0.83045, |
> 0.04731d0, -0.14069d0, 0.33899d0, 0.74002d0, 0.62827d0, -0.83045d0, |
> -17.17104, 32.13834, -3.8414, -13.50958, 37.48258, -11.74448, |
> -17.17104d0, 32.13834d0, -3.8414d0, -13.50958d0, 37.48258d0, -11.74448d0, |
> 0.05221, -0.11177, 0.27687, 0.68277, 0.62925, 2.10930, |
> 0.05221d0, -0.11177d0, 0.27687d0, 0.68277d0, 0.62925d0, 2.10930d0, |
> 13.86561, -1.66573, 7.9977, 0.35896, 27.99104, -32.29613, |
> 13.86561d0, -1.66573d0, 7.9977d0, 0.35896d0, 27.99104d0, -32.29613d0, |
> 0.06442, -0.07257, 0.11952, 0.50396, 0.53317, 4.29581, |
> 0.06442d0, -0.07257d0, 0.11952d0, 0.50396d0, 0.53317d0, 4.29581d0, |
> 44.07500, -35.61928, 17.2312, 19.12265, 0.86142, -37.08731, |
> 44.07500d0, -35.61928d0, 17.2312d0, 19.12265d0, 0.86142d0, -37.08731d0, |
> 0.07351, -0.10911, 0.12156, 0.46602, 0.48331, 5.17804, |
> 0.07351d0, -0.10911d0, 0.12156d0, 0.46602d0, 0.48331d0, 5.17804d0, |
> 58.31429, -50.40657, 21.0203, 27.73335, -10.71966, -39.46624, |
> 58.31429d0, -50.40657d0, 21.0203d0, 27.73335d0, -10.71966d0, -39.46624d0, |
> 0.06943, -0.07336, 0.15002, 0.46675, 0.51645, 3.31844, |
> 0.06943d0, -0.07336d0, 0.15002d0, 0.46675d0, 0.51645d0, 3.31844d0, |
> 37.30048, -28.68561, 13.4842, 17.16954, -2.77701, -26.91567, |
> 37.30048d0, -28.68561d0, 13.4842d0, 17.16954d0, -2.77701d0, -26.91567d0, |
> 0.05361, -0.02633, 0.19376, 0.44691, 0.44814, 2.88220, |
> 0.05361d0, -0.02633d0, 0.19376d0, 0.44691d0, 0.44814d0, 2.88220d0, |
> 29.95785, -19.93872, 11.4738, 10.83239, 7.47675, -28.12408, |
> 29.95785d0, -19.93872d0, 11.4738d0, 10.83239d0, 7.47675d0, -28.12408d0, |
> 0.03795, -0.04353, 0.30496, 0.53973, 0.44511, 1.82131, |
> 0.03795d0, -0.04353d0, 0.30496d0, 0.53973d0, 0.44511d0, 1.82131d0, |
> 14.62905, -0.35599, 6.9752, 1.48548, 24.99437, -29.88984, |
> 14.62905d0, -0.35599d0, 6.9752d0, 1.48548d0, 24.99437d0, -29.88984d0, |
> 0.04790, 0.05632, 0.20660, 0.39304, 0.45673, 4.19462, |
> 0.04790d0, 0.05632d0, 0.20660d0, 0.39304d0, 0.45673d0, 4.19462d0, |
> 47.41709, -39.83233, 16.6487, 17.18184, -2.52659, -32.65720, |
> 47.41709d0, -39.83233d0, 16.6487d0, 17.18184d0, -2.52659d0, -32.65720d0, |
> 0.05140, 0.11201, 0.19756, 0.26992, 0.30311, 3.20459, |
> 0.05140d0, 0.11201d0, 0.19756d0, 0.26992d0, 0.30311d0, 3.20459d0, |
> 37.47848, -31.60155, 12.7851, 13.97016, -8.18266, -18.63125, |
> 37.47848d0, -31.60155d0, 12.7851d0, 13.97016d0, -8.18266d0, -18.63125d0, |
> 0.04703, 0.17203, 0.21807, 0.23841, 0.26405, 1.17976, |
> 0.04703d0, 0.17203d0, 0.21807d0, 0.23841d0, 0.26405d0, 1.17976d0, |
> 10.84169, -4.87255, 4.5850, 2.60271, -1.08836, -2.73421, |
> 10.84169d0, -4.87255d0, 4.5850d0, 2.60271d0, -1.08836d0, -2.73421d0, |
> 0.09651, -0.07097, -0.01140, 0.41291, 0.67944, 8.63177, |
> 0.09651d0, -0.07097d0, -0.01140d0, 0.41291d0, 0.67944d0, 8.63177d0, |
> 116.48230,-115.22170, 34.6775, 50.71781, -62.71606, -31.72991, |
> 116.48230d0,-115.22170d0, 34.6775d0, 50.71781d0, -62.71606d0, -31.72991d0, |
> 0.05407, -0.01853, 0.18484, 0.39172, 0.39883, 4.54709, |
> 0.05407d0, -0.01853d0, 0.18484d0, 0.39172d0, 0.39883d0, 4.54709d0, |
> 57.66357, -47.37784, 18.0009, 22.01861, -16.19442, -25.28908, |
> 57.66357d0, -47.37784d0, 18.0009d0, 22.01861d0, -16.19442d0, -25.28908d0, |
> 0.05646, -0.02226, 0.19051, 0.38210, 0.35937, 2.35732, |
> 0.05646d0, -0.02226d0, 0.19051d0, 0.38210d0, 0.35937d0, 2.35732d0, |
> 31.71869, -21.30330, 9.5296, 14.48209, -14.30969, -7.90519, |
> 31.71869d0, -21.30330d0, 9.5296d0, 14.48209d0, -14.30969d0, -7.90519d0, |
> 0.04477, 0.00480, 0.15076, 0.29189, 0.31529, 3.12081, |
> 0.04477d0, 0.00480d0, 0.15076d0, 0.29189d0, 0.31529d0, 3.12081d0, |
> 41.81401, -31.97399, 12.4581, 15.71344, -15.68024, -13.09181, |
> 41.81401d0, -31.97399d0, 12.4581d0, 15.71344d0, -15.68024d0, -13.09181d0, |
> 0.03143, 0.06733, 0.18032, 0.27882, 0.27721, 1.09240, |
> 0.03143d0, 0.06733d0, 0.18032d0, 0.27882d0, 0.27721d0, 1.09240d0, |
> 15.14939, -7.18577, 4.4895, 7.48708, -11.64642, 1.23596, |
> 15.14939d0, -7.18577d0, 4.4895d0, 7.48708d0, -11.64642d0, 1.23596d0, |
> 0.09641, -0.07037, -0.04597, 0.26390, 0.48362, 3.68185, |
> 0.09641d0, -0.07037d0, -0.04597d0, 0.26390d0, 0.48362d0, 3.68185d0, |
> 61.90438, -66.11752, 15.0943, 29.88299, -76.76054, 30.73642, |
> 61.90438d0, -66.11752d0, 15.0943d0, 29.88299d0, -76.76054d0, 30.73642d0, |
> 0.05348, -0.06567, 0.26431, 0.51672, -0.39508, -0.17103, |
> 0.05348d0, -0.06567d0, 0.26431d0, 0.51672d0, -0.39508d0, -0.17103d0, |
> 0.38781, 9.77783, -0.2076, 3.74785, -13.89592, 15.93324, |
> 0.38781d0, 9.77783d0, -0.2076d0, 3.74785d0, -13.89592d0, 15.93324d0, |
> 0.01978, 0.02381, 0.16282, 0.27720, 0.28227, 1.10671, |
> 0.01978d0, 0.02381d0, 0.16282d0, 0.27720d0, 0.28227d0, 1.10671d0, |
> 19.57988, -8.77431, 4.6923, 10.90089, -21.92500, 4.82789, |
> 19.57988d0, -8.77431d0, 4.6923d0, 10.90089d0, -21.92500d0, 4.82789d0, |
> 0.02816, 0.08035, 0.15177, 0.19596, 0.20729, 0.16046, |
> 0.02816d0, 0.08035d0, 0.15177d0, 0.19596d0, 0.20729d0, 0.16046d0, |
> 4.85393, 1.60025, 0.9924, 6.34038, -23.91168, 19.81668, |
> 4.85393d0, 1.60025d0, 0.9924d0, 6.34038d0, -23.91168d0, 19.81668d0, |
> 0.00810, 0.06904, 0.19159, 0.22518, 0.16841, 0.22909, |
> 0.00810d0, 0.06904d0, 0.19159d0, 0.22518d0, 0.16841d0, 0.22909d0, |
> 4.46808, -1.67027, 1.1692, 4.17188, -22.18001, 13.56817, |
> 4.46808d0, -1.67027d0, 1.1692d0, 4.17188d0, -22.18001d0, 13.56817d0, |
> -0.00459, 0.08266, 0.20067, 0.22163, 0.15112, 0.67785, |
> -0.00459d0, 0.08266d0, 0.20067d0, 0.22163d0, 0.15112d0, 0.67785d0, |
> 10.03290, -23.97605, 2.8133, 1.49803, -32.05834, 14.49983, |
> 10.03290d0, -23.97605d0, 2.8133d0, 1.49803d0, -32.05834d0, 14.49983d0, |
> -0.00869, 0.07988, 0.20642, 0.23439, 0.16011, 1.58356, |
> -0.00869d0, 0.07988d0, 0.20642d0, 0.23439d0, 0.16011d0, 1.58356d0, |
> 23.98258, -56.34460, 6.2639, -0.04471, -43.27986, 13.87953/ |
> 23.98258d0, -56.34460d0, 6.2639d0, -0.04471d0, -43.27986d0, 13.87953d0/ |
endif | endif |
| |
c if (ww .ge. 3.9) then | c if (ww .ge. 3.9) then |
|
|
| |
c if (ww.le.1.2 .or. ww.ge.3.85 | c if (ww.le.1.2 .or. ww.ge.3.85 |
c >.or. abs(ww-ww1).le.0.001) then | c >.or. abs(ww-ww1).le.0.001) then |
if (ww.le.1.2 .or. abs(ww-ww1).le.0.001) then |
if (ww.le.1.2d0 .or. abs(ww-ww1).le.0.001d0) then |
r = r1 | r = r1 |
f2 = f21 | f2 = f21 |
else | else |
|
|
rlog = fac/tmp2 | rlog = fac/tmp2 |
| |
r2 = rp(i)*rlog | r2 = rp(i)*rlog |
&+(rp(i+1)/qq+rp(i+2)/(q4+.3**2)) |
&+(rp(i+1)/qq+rp(i+2)/(q4+.3d0**2)) |
&*(1.+rp(i+3)*xx+rp(i+4)*(xx**2)) |
&*(1.d0+rp(i+3)*xx+rp(i+4)*(xx**2)) |
&*(xx**rp(i+5)) | &*(xx**rp(i+5)) |
| |
q = log(qq) | q = log(qq) |
qth = .2d0+3.2d0*xx | qth = .2d0+3.2d0*xx |
dq = q-qth | dq = q-qth |
y = 1.-xx |
y = 1.d0-xx |
f2th = f2p(i1)*y**1+f2p(i1+1)*y**2+f2p(i1+2)*y**3+ | f2th = f2p(i1)*y**1+f2p(i1+1)*y**2+f2p(i1+2)*y**3+ |
> f2p(i1+3)*y**4+f2p(i1+4)*y**5 | > f2p(i1+3)*y**4+f2p(i1+4)*y**5 |
quad = (f2p(i1+5)+f2p(i1+6)*xx+f2p(i1+7)*xx**2)*dq**2 | quad = (f2p(i1+5)+f2p(i1+6)*xx+f2p(i1+7)*xx**2)*dq**2 |