|
|
CALL GINIT(400,800)
SET WINDOW 0,5,-1,49
DRAW GRID(1,10)
FOR N=0 TO 5
SET LINE COLOR N+1
FOR X=0 TO 5 STEP 1/16
LET Y=BESSELI(N,X)
PLOT LINES:X,Y;
NEXT X
PLOT LINES
NEXT N
END
EXTERNAL SUB GINIT(XSIZE,YSIZE)
SET BITMAP SIZE XSIZE,YSIZE
SET WINDOW 0,XSIZE-1,YSIZE-1,0
SET POINT STYLE 1
SET COLOR MODE "REGULAR"
SET COLOR MIX(0) 0,0,0
SET COLOR MIX(1) 0,0,1
SET COLOR MIX(2) 1,0,0
SET COLOR MIX(3) 1,0,1
SET COLOR MIX(4) 0,1,0
SET COLOR MIX(5) 0,1,1
SET COLOR MIX(6) 1,1,0
SET COLOR MIX(7) 1,1,1
CLEAR
END SUB
EXTERNAL FUNCTION INTEGRAL(R,Z) !'ガウス・ルジャンドル
LET A=0
LET B=PI
LET U=(B+A)/2
LET V=(B-A)/2
FOR I=1 TO 400
READ X,C
LET S=S+C*F(R,Z,U+V*X)*V
NEXT I
LET INTEGRAL=S/PI
DATA -.9999819727039624507107997,4.6263724177190118157440220E-05
DATA -.9999050164693312795671775,1.0769038102183646230364914E-04
DATA -.9997665712667108498622001,1.6920137336541413264932958E-04
DATA -.9995666148626408000179844,2.3070985544248655269465132E-04
DATA -.9993051557937030335325083,2.9220572803523217847223343E-04
DATA -.9989822092480797992371162,3.5368409613102448239639817E-04
DATA -.9985977947971757433952425,4.1514088233419025171823795E-04
DATA -.9981519359718441275385541,4.7657220501003854013133061E-04
DATA -.9976446601486403940368916,5.3797424394626279849881322E-04
DATA -.9970759985107923450247946,5.9934320270831884264722444E-04
DATA -.9964459860317571225964113,6.6067529607738161575318920E-04
DATA -.9957546614667239851797758,7.2196674533702894531306739E-04
DATA -.9950020673471962444031388,7.8321377638390310244817342E-04
DATA -.9941882499768158116647105,8.4441261897085745246513294E-04
DATA -.9933132594276647813762767,9.0555950644549534547199634E-04
DATA -.9923771495366982433589461,9.6665067572141452408517998E-04
DATA -.9913799779021414417532216,1.0276823673654497630797823E-03
DATA -.9903218058797661642331360,1.0886508257457038100719739E-03
DATA -.9892026985790008762180358,1.1495522992128106211745720E-03
DATA -.9880227248588493179127030,1.2103830403000240185954168E-03
DATA -.9867819573236030434401229,1.2711393059342868032230354E-03
DATA -.9854804723183393394746898,1.3318173576538514742077422E-03
DATA -.9841183499241993752558070,1.3924134618298698080684263E-03
DATA -.9826956739534434545500623,1.4529238898903999088514186E-03
DATA -.9812125319442814682424151,1.5133449185458729415321843E-03
DATA -.9796690151554774136346128,1.5736728300154127016475733E-03
DATA -.9780652185607273384489774,1.6339039122536137611976595E-03
DATA -.9764012408428103908462605,1.6940344591775155933476319E-03
DATA -.9746771843875128746551103,1.7540607708935932918534161E-03
DATA -.9728931552773253596853138,1.8139791539246390831213050E-03
DATA -.9710492632849130040397740,1.8737859214364439296000032E-03
DATA -.9691456218663593235845560,1.9334773934642118869862805E-03
DATA -.9671823481541837026768633,1.9930498971386556619114334E-03
DATA -.9651595629501329860638605,2.0524997669117326296387785E-03
DATA -.9630773907177475286294707,2.1118233447819880761749695E-03
DATA -.9609359595747021101129220,2.1710169805194777030769764E-03
DATA -.9587354012849221479049303,2.2300770318902451834521819E-03
DATA -.9564758512504756638175957,2.2889998648803332549695712E-03
DATA -.9541574485032414812063787,2.3477818539193087973829979E-03
DATA -.9517803356963541476176046,2.4064193821032837844401297E-03
DATA -.9493446590954260956830096,2.4649088414174150722575315E-03
DATA -.9468505685695475715990344,2.5232466329578667921869215E-03
DATA -.9442982175820648764404846,2.5814291671532197289820420E-03
DATA -.9416877631811374809315652,2.6394528639853125366129890E-03
DATA -.9390193659900745892537874,2.6973141532095000115583202E-03
DATA -.9362931901974517420997422,2.7550094745743139337060721E-03
DATA -.9335094035470080635511228,2.8125352780405122176822307E-03
DATA -.9306681773273247705187315,2.8698880239995023068591947E-03
DATA -.9277696863612855774691976,2.9270641834911248990910487E-03
DATA -.9248141089953196430062722,2.9840602384207842252760188E-03
DATA -.9218016270884277185954871,3.0408726817759112150766292E-03
DATA -.9187324260009921733361549,3.0974980178417459829955126E-03
DATA -.9156066945833715822067555,3.1539327624164261558923710E-03
DATA -.9124246251642805786483078,3.2101734430253676425153532E-03
DATA -.9091864135389556857127581,3.2662165991349245186867942E-03
DATA -.9058922589571078532953371,3.3220587823653147699542163E-03
DATA -.9025423641106624421955255,3.3776965567027986979745151E-03
DATA -.8991369351212874089139762,3.4331264987110968585623238E-03
DATA -.8956761815277104581948603,3.4883451977420344589220515E-03
DATA -.8921603162728259433664469,3.5433492561453991996620077E-03
DATA -.8885895556905923075185429,3.5981352894779996042038747E-03
DATA -.8849641194927208714846091,3.6526999267119109345042264E-03
DATA -.8812842307551567874694922,3.7070398104418958478769910E-03
DATA -.8775501159043529899810783,3.7611515970919870053700346E-03
DATA -.8737620047033379884859100,3.8150319571212188977839493E-03
DATA -.8699201302375783589148594,3.8686775752284962111698272E-03
DATA -.8660247289006368037951485,3.9220851505565861096203416E-03
DATA -.8620760403796266633790730,3.9752513968952218694674985E-03
DATA -.8580743076404637726773631,4.0281730428833053556957251E-03
DATA -.8540197769129165717857668,4.0808468322101958885308878E-03
DATA -.8499126976754553893166989,4.1332695238160731058231723E-03
DATA -.8457533226399018311131247,4.1854378920913614850439286E-03
DATA -.8415419077358792187286900,4.2373487270752042474954035E-03
DATA -.8372787120950650344058650,4.2889988346529744267133843E-03
DATA -.8329639980352463414719476,4.3403850367528109430467473E-03
DATA -.8285980310441791612005141,4.3915041715411675870411756E-03
DATA -.8241810797632527992526874,4.4423530936173628755499398E-03
DATA -.8197134159709601268177470,4.4929286742071188064528328E-03
DATA -.8151953145661748335154631,4.5432278013550766004930588E-03
DATA -.8106270535512366810024356,4.5932473801162775820470507E-03
DATA -.8060089140148457980409629,4.6429843327465974146282383E-03
DATA -.8013411801147670695408899,4.6924355988921219715957326E-03
DATA -.7966241390603456837717921,4.7415981357774531878946569E-03
DATA -.7918580810948349135640607,4.7904689183929333046974074E-03
DATA -.7870432994775372188722786,4.8390449396807759855446397E-03
DATA -.7821800904657597695620265,4.8873232107200928500007496E-03
DATA -.7772687532965854987012501,4.9353007609108040389399652E-03
DATA -.7723095901684608079888676,4.9829746381564214943643207E-03
DATA -.7673029062226010582357195,5.0303419090456937061217566E-03
DATA -.7622490095242149890255755,5.0773996590331007480390167E-03
DATA -.7571482110435492228260470,5.1241449926181884968072417E-03
DATA -.7520008246367540198902183,5.1705750335237309984552282E-03
DATA -.7468071670265714612889527,5.2166869248727100194130774E-03
DATA -.7415675577828472483404577,5.2624778293641008920039892E-03
DATA -.7362823193028673175571692,5.3079449294474538377007114E-03
DATA -.7309517767915204810096465,5.3530854274962600256418876E-03
DATA -.7255762582412883127123272,5.3978965459800916987183857E-03
DATA -.7201560944120635122659979,5.4423755276355057750066148E-03
DATA -.7146916188107979875460589,5.4865196356357004084407400E-03
DATA -.7091831676709819087034412,5.5303261537589140693743281E-03
DATA -.7036310799319549961457379,5.5737923865555567830799909E-03
DATA -.6980356972180513154891026,5.6169156595140632422686060E-03
DATA -.6923973638175788627161007,5.6596933192254575883731863E-03
DATA -.6867164266616352329403680,5.7021227335466197356318722E-03
DATA -.6809932353027606762649893,5.7442012917622431919151687E-03
DATA -.6752281418934298542273508,5.7859264047454744107697479E-03
DATA -.6694215011643836202482341,5.8272955051172237902900934E-03
DATA -.6635736704028021573464804,5.8683060474041385161751495E-03
DATA -.6576850094303208161420752,5.9089555081952275286750641E-03
DATA -.6517558805808900058493764,5.9492413862971289760781243E-03
DATA -.6457866486784805005578351,5.9891612028880106009251001E-03
DATA -.6397776810146355326093649,6.0287125016700935892623685E-03
DATA -.6337293473258710543088905,6.0678928490207904979513170E-03
DATA -.6276420197709255585470036,6.1066998341424479603344647E-03
DATA -.6215160729078608581704682,6.1451310692106849564133301E-03
DATA -.6153518836710152331070009,6.1831841895213175201140730E-03
DATA -.6091498313478103633347320,6.2208568536358608431990998E-03
DATA -.6029102975554134747834702,6.2581467435255998229208227E-03
DATA -.5966336662172561341638003,6.2950515647142191886022727E-03
DATA -.5903203235394111375405864,6.3315690464189844309628949E-03
DATA -.5839706579868289461990941,6.3676969416904648471811825E-03
DATA -.5775850602594351319941408,6.4034330275507901043933791E-03
DATA -.5711639232680903029249140,6.4387751051304318145637988E-03
DATA -.5647076421104139881398226,6.4737209998035017044218583E-03
DATA -.5582166140464739699464555,6.5082685613215580554381159E-03
DATA -.5516912384743425586808994,6.5424156639459121806008761E-03
DATA -.5451319169055213144777974,6.5761602065784267970506106E-03
DATA -.5385390529402357280771185,6.6095001128907982464259091E-03
DATA -.5319130522426013808051484,6.6424333314523146080662055E-03
DATA -.5252543225156631117752170,6.6749578358560818439973984E-03
DATA -.5185632734763087281676594,6.7070716248437102088909494E-03
DATA -.5118403168300588021679862,6.7387727224284532529292952E-03
DATA -.5050858662457341057667365,6.7700591780167918407246379E-03
DATA -.4983003373300022421535401,6.8009290665284557051185368E-03
DATA -.4914841476018050398710538,6.8313804885148751508303160E-03
DATA -.4846377164666682832312103,6.8614115702760556195172418E-03
DATA -.4777614651908953597361706,6.8910204639758679248527525E-03
DATA -.4708558168756464123890576,6.9202053477557470637147884E-03
DATA -.4639211964309045918245310,6.9489644258467926074984745E-03
DATA -.4569580305493310101361121,6.9772959286802637759200314E-03
DATA -.4499667476800100051254408,7.0051981129964623944557826E-03
DATA -.4429477780020863304479467,7.0326692619519970357554193E-03
DATA -.4359015533982958937793024,7.0597076852254217449761731E-03
DATA -.4288285074283916716771105,7.0863117191212428489981040E-03
DATA -.4217290753024664362621394,7.1124797266722874498941965E-03
DATA -.4146036938541739351926808,7.1382100977404273038361762E-03
DATA -.4074528015138501726538528,7.1635012491156518878117343E-03
DATA -.4002768382815364452305435,7.1883516246134845581049307E-03
DATA -.3930762456999057925777893,7.2127596951707358064427126E-03
DATA -.3858514668270945287453553,7.2367239589395877220304415E-03
DATA -.3786029462094405258537453,7.2602429413800038703817938E-03
DATA -.3713311298541299275564774,7.2833151953504589028870525E-03
DATA -.3640364652017539753578509,7.3059393011969823144523193E-03
DATA -.3567194010987776363862590,7.3281138668405108702741736E-03
DATA -.3493803877699217266500321,7.3498375278625443268834112E-03
DATA -.3420198767904602291254888,7.3711089475890991769913129E-03
DATA -.3346383210584345112450044,7.3919268171729552523959906E-03
DATA -.3272361747667861514661529,7.4122898556741901242482968E-03
DATA -.3198138933754100896110293,7.4321968101389963453301042E-03
DATA -.3123719335831298205674013,7.4516464556767766846559753E-03
DATA -.3049107532995963557400760,7.4706375955355126106658549E-03
DATA -.2974308116171126813314970,7.4891690611754013855249070E-03
DATA -.2899325687823854471148313,7.5072397123407572395804471E-03
DATA -.2824164861682056238403682,7.5248484371301722018385363E-03
DATA -.2748830262450598717866731,7.5419941520649322694076350E-03
DATA -.2673326525526743672313363,7.5586758021556847062071669E-03
DATA -.2597658296714928377720767,7.5748923609673523688483221E-03
DATA -.2521830231940905614771423,7.5906428306822910654563041E-03
DATA -.2445846996965260887841458,7.6059262421616860613108628E-03
DATA -.2369713267096324498984454,7.6207416550051839535287106E-03
DATA -.2293433726902496141656847,7.6350881576087562455906106E-03
DATA -.2217013069923999715079185,7.6489648672207910613208914E-03
DATA -.2140455998384086095186518,7.6623709299964095469511727E-03
DATA -.2063767222899701632088853,7.6753055210500036191364789E-03
DATA -.1986951462191640176836869,7.6877678445059918262339545E-03
DATA -.1910013442794196472066982,7.6997571335477901997953341E-03
DATA -.1832957898764338771781292,7.7112726504649950830574263E-03
DATA -.1755789571390418585100171,7.7223136866987750332333719E-03
DATA -.1678513208900435467306408,7.7328795628854690046045805E-03
DATA -.1601133566169874808878088,7.7429696288983881297822398E-03
DATA -.1523655404429136599481171,7.7525832638878185270413483E-03
DATA -.1446083490970573169060391,7.7617198763192226723225395E-03
DATA -.1368422598855153932227026,7.7703789040096369853407298E-03
DATA -.1290677506618775185092873,7.7785598141622633902280246E-03
DATA -.1212852997978233025539997,7.7862621033992527222645108E-03
DATA -.1134953861536877488644070,7.7934852977926779635077296E-03
DATA -.1056984890489966008584237,7.8002289528936954015129025E-03
DATA -.0978950882329734336873092,7.8064926537598919168345279E-03
DATA -.0900856638550203064125526,7.8122760149808167166089168E-03
DATA -.0822706964351737908853691,7.8175786807016959432297203E-03
DATA -.0744506668345381951926272,7.8224003246453286989376515E-03
DATA -.0666260562256978009362603,7.8267406501321631390445391E-03
DATA -.0587973460631099349045127,7.8305993900985513984936792E-03
DATA -.0509650180534806968726444,7.8339763071131822285162928E-03
DATA -.0431295541261251663379026,7.8368711933916903322708496E-03
DATA -.0352914364033139119485902,7.8392838708094415005411850E-03
DATA -.0274511471706076282298676,7.8412141909124927608138163E-03
DATA -.0196091688471817249394630,7.8426620349267268653477464E-03
DATA -.0117659839561426950046856,7.8436273137651605561844209E-03
DATA -.0039220750948380874980108,7.8441099680334261574144661E-03
DATA .0039220750948380874980108,7.8441099680334261574144661E-03
DATA .0117659839561426950046856,7.8436273137651605561844209E-03
DATA .0196091688471817249394630,7.8426620349267268653477464E-03
DATA .0274511471706076282298676,7.8412141909124927608138163E-03
DATA .0352914364033139119485902,7.8392838708094415005411850E-03
DATA .0431295541261251663379026,7.8368711933916903322708496E-03
DATA .0509650180534806968726444,7.8339763071131822285162928E-03
DATA .0587973460631099349045127,7.8305993900985513984936792E-03
DATA .0666260562256978009362603,7.8267406501321631390445391E-03
DATA .0744506668345381951926272,7.8224003246453286989376515E-03
DATA .0822706964351737908853691,7.8175786807016959432297203E-03
DATA .0900856638550203064125526,7.8122760149808167166089168E-03
DATA .0978950882329734336873092,7.8064926537598919168345279E-03
DATA .1056984890489966008584237,7.8002289528936954015129025E-03
DATA .1134953861536877488644070,7.7934852977926779635077296E-03
DATA .1212852997978233025539997,7.7862621033992527222645108E-03
DATA .1290677506618775185092873,7.7785598141622633902280246E-03
DATA .1368422598855153932227026,7.7703789040096369853407298E-03
DATA .1446083490970573169060391,7.7617198763192226723225395E-03
DATA .1523655404429136599481171,7.7525832638878185270413483E-03
DATA .1601133566169874808878088,7.7429696288983881297822398E-03
DATA .1678513208900435467306408,7.7328795628854690046045805E-03
DATA .1755789571390418585100171,7.7223136866987750332333719E-03
DATA .1832957898764338771781292,7.7112726504649950830574263E-03
DATA .1910013442794196472066982,7.6997571335477901997953341E-03
DATA .1986951462191640176836869,7.6877678445059918262339545E-03
DATA .2063767222899701632088853,7.6753055210500036191364789E-03
DATA .2140455998384086095186518,7.6623709299964095469511727E-03
DATA .2217013069923999715079185,7.6489648672207910613208914E-03
DATA .2293433726902496141656847,7.6350881576087562455906106E-03
DATA .2369713267096324498984454,7.6207416550051839535287106E-03
DATA .2445846996965260887841458,7.6059262421616860613108628E-03
DATA .2521830231940905614771423,7.5906428306822910654563041E-03
DATA .2597658296714928377720767,7.5748923609673523688483221E-03
DATA .2673326525526743672313363,7.5586758021556847062071669E-03
DATA .2748830262450598717866731,7.5419941520649322694076350E-03
DATA .2824164861682056238403682,7.5248484371301722018385363E-03
DATA .2899325687823854471148313,7.5072397123407572395804471E-03
DATA .2974308116171126813314970,7.4891690611754013855249070E-03
DATA .3049107532995963557400760,7.4706375955355126106658549E-03
DATA .3123719335831298205674013,7.4516464556767766846559753E-03
DATA .3198138933754100896110293,7.4321968101389963453301042E-03
DATA .3272361747667861514661529,7.4122898556741901242482968E-03
DATA .3346383210584345112450044,7.3919268171729552523959906E-03
DATA .3420198767904602291254888,7.3711089475890991769913129E-03
DATA .3493803877699217266500321,7.3498375278625443268834112E-03
DATA .3567194010987776363862590,7.3281138668405108702741736E-03
DATA .3640364652017539753578509,7.3059393011969823144523193E-03
DATA .3713311298541299275564774,7.2833151953504589028870525E-03
DATA .3786029462094405258537453,7.2602429413800038703817938E-03
DATA .3858514668270945287453553,7.2367239589395877220304415E-03
DATA .3930762456999057925777893,7.2127596951707358064427126E-03
DATA .4002768382815364452305435,7.1883516246134845581049307E-03
DATA .4074528015138501726538528,7.1635012491156518878117343E-03
DATA .4146036938541739351926808,7.1382100977404273038361762E-03
DATA .4217290753024664362621394,7.1124797266722874498941965E-03
DATA .4288285074283916716771105,7.0863117191212428489981040E-03
DATA .4359015533982958937793024,7.0597076852254217449761731E-03
DATA .4429477780020863304479467,7.0326692619519970357554193E-03
DATA .4499667476800100051254408,7.0051981129964623944557826E-03
DATA .4569580305493310101361121,6.9772959286802637759200314E-03
DATA .4639211964309045918245310,6.9489644258467926074984745E-03
DATA .4708558168756464123890576,6.9202053477557470637147884E-03
DATA .4777614651908953597361706,6.8910204639758679248527525E-03
DATA .4846377164666682832312103,6.8614115702760556195172418E-03
DATA .4914841476018050398710538,6.8313804885148751508303160E-03
DATA .4983003373300022421535401,6.8009290665284557051185368E-03
DATA .5050858662457341057667365,6.7700591780167918407246379E-03
DATA .5118403168300588021679862,6.7387727224284532529292952E-03
DATA .5185632734763087281676594,6.7070716248437102088909494E-03
DATA .5252543225156631117752170,6.6749578358560818439973984E-03
DATA .5319130522426013808051484,6.6424333314523146080662055E-03
DATA .5385390529402357280771185,6.6095001128907982464259091E-03
DATA .5451319169055213144777974,6.5761602065784267970506106E-03
DATA .5516912384743425586808994,6.5424156639459121806008761E-03
DATA .5582166140464739699464555,6.5082685613215580554381159E-03
DATA .5647076421104139881398226,6.4737209998035017044218583E-03
DATA .5711639232680903029249140,6.4387751051304318145637988E-03
DATA .5775850602594351319941408,6.4034330275507901043933791E-03
DATA .5839706579868289461990941,6.3676969416904648471811825E-03
DATA .5903203235394111375405864,6.3315690464189844309628949E-03
DATA .5966336662172561341638003,6.2950515647142191886022727E-03
DATA .6029102975554134747834702,6.2581467435255998229208227E-03
DATA .6091498313478103633347320,6.2208568536358608431990998E-03
DATA .6153518836710152331070009,6.1831841895213175201140730E-03
DATA .6215160729078608581704682,6.1451310692106849564133301E-03
DATA .6276420197709255585470036,6.1066998341424479603344647E-03
DATA .6337293473258710543088905,6.0678928490207904979513170E-03
DATA .6397776810146355326093649,6.0287125016700935892623685E-03
DATA .6457866486784805005578351,5.9891612028880106009251001E-03
DATA .6517558805808900058493764,5.9492413862971289760781243E-03
DATA .6576850094303208161420752,5.9089555081952275286750641E-03
DATA .6635736704028021573464804,5.8683060474041385161751495E-03
DATA .6694215011643836202482341,5.8272955051172237902900934E-03
DATA .6752281418934298542273508,5.7859264047454744107697479E-03
DATA .6809932353027606762649893,5.7442012917622431919151687E-03
DATA .6867164266616352329403680,5.7021227335466197356318722E-03
DATA .6923973638175788627161007,5.6596933192254575883731863E-03
DATA .6980356972180513154891026,5.6169156595140632422686060E-03
DATA .7036310799319549961457379,5.5737923865555567830799909E-03
DATA .7091831676709819087034412,5.5303261537589140693743281E-03
DATA .7146916188107979875460589,5.4865196356357004084407400E-03
DATA .7201560944120635122659979,5.4423755276355057750066148E-03
DATA .7255762582412883127123272,5.3978965459800916987183857E-03
DATA .7309517767915204810096465,5.3530854274962600256418876E-03
DATA .7362823193028673175571692,5.3079449294474538377007114E-03
DATA .7415675577828472483404577,5.2624778293641008920039892E-03
DATA .7468071670265714612889527,5.2166869248727100194130774E-03
DATA .7520008246367540198902183,5.1705750335237309984552282E-03
DATA .7571482110435492228260470,5.1241449926181884968072417E-03
DATA .7622490095242149890255755,5.0773996590331007480390167E-03
DATA .7673029062226010582357195,5.0303419090456937061217566E-03
DATA .7723095901684608079888676,4.9829746381564214943643207E-03
DATA .7772687532965854987012501,4.9353007609108040389399652E-03
DATA .7821800904657597695620265,4.8873232107200928500007496E-03
DATA .7870432994775372188722786,4.8390449396807759855446397E-03
DATA .7918580810948349135640607,4.7904689183929333046974074E-03
DATA .7966241390603456837717921,4.7415981357774531878946569E-03
DATA .8013411801147670695408899,4.6924355988921219715957326E-03
DATA .8060089140148457980409629,4.6429843327465974146282383E-03
DATA .8106270535512366810024356,4.5932473801162775820470507E-03
DATA .8151953145661748335154631,4.5432278013550766004930588E-03
DATA .8197134159709601268177470,4.4929286742071188064528328E-03
DATA .8241810797632527992526874,4.4423530936173628755499398E-03
DATA .8285980310441791612005141,4.3915041715411675870411756E-03
DATA .8329639980352463414719476,4.3403850367528109430467473E-03
DATA .8372787120950650344058650,4.2889988346529744267133843E-03
DATA .8415419077358792187286900,4.2373487270752042474954035E-03
DATA .8457533226399018311131247,4.1854378920913614850439286E-03
DATA .8499126976754553893166989,4.1332695238160731058231723E-03
DATA .8540197769129165717857668,4.0808468322101958885308878E-03
DATA .8580743076404637726773631,4.0281730428833053556957251E-03
DATA .8620760403796266633790730,3.9752513968952218694674985E-03
DATA .8660247289006368037951485,3.9220851505565861096203416E-03
DATA .8699201302375783589148594,3.8686775752284962111698272E-03
DATA .8737620047033379884859100,3.8150319571212188977839493E-03
DATA .8775501159043529899810783,3.7611515970919870053700346E-03
DATA .8812842307551567874694922,3.7070398104418958478769910E-03
DATA .8849641194927208714846091,3.6526999267119109345042264E-03
DATA .8885895556905923075185429,3.5981352894779996042038747E-03
DATA .8921603162728259433664469,3.5433492561453991996620077E-03
DATA .8956761815277104581948603,3.4883451977420344589220515E-03
DATA .8991369351212874089139762,3.4331264987110968585623238E-03
DATA .9025423641106624421955255,3.3776965567027986979745151E-03
DATA .9058922589571078532953371,3.3220587823653147699542163E-03
DATA .9091864135389556857127581,3.2662165991349245186867942E-03
DATA .9124246251642805786483078,3.2101734430253676425153532E-03
DATA .9156066945833715822067555,3.1539327624164261558923710E-03
DATA .9187324260009921733361549,3.0974980178417459829955126E-03
DATA .9218016270884277185954871,3.0408726817759112150766292E-03
DATA .9248141089953196430062722,2.9840602384207842252760188E-03
DATA .9277696863612855774691976,2.9270641834911248990910487E-03
DATA .9306681773273247705187315,2.8698880239995023068591947E-03
DATA .9335094035470080635511228,2.8125352780405122176822307E-03
DATA .9362931901974517420997422,2.7550094745743139337060721E-03
DATA .9390193659900745892537874,2.6973141532095000115583202E-03
DATA .9416877631811374809315652,2.6394528639853125366129890E-03
DATA .9442982175820648764404846,2.5814291671532197289820420E-03
DATA .9468505685695475715990344,2.5232466329578667921869215E-03
DATA .9493446590954260956830096,2.4649088414174150722575315E-03
DATA .9517803356963541476176046,2.4064193821032837844401297E-03
DATA .9541574485032414812063787,2.3477818539193087973829979E-03
DATA .9564758512504756638175957,2.2889998648803332549695712E-03
DATA .9587354012849221479049303,2.2300770318902451834521819E-03
DATA .9609359595747021101129220,2.1710169805194777030769764E-03
DATA .9630773907177475286294707,2.1118233447819880761749695E-03
DATA .9651595629501329860638605,2.0524997669117326296387785E-03
DATA .9671823481541837026768633,1.9930498971386556619114334E-03
DATA .9691456218663593235845560,1.9334773934642118869862805E-03
DATA .9710492632849130040397740,1.8737859214364439296000032E-03
DATA .9728931552773253596853138,1.8139791539246390831213050E-03
DATA .9746771843875128746551103,1.7540607708935932918534161E-03
DATA .9764012408428103908462605,1.6940344591775155933476319E-03
DATA .9780652185607273384489774,1.6339039122536137611976595E-03
DATA .9796690151554774136346128,1.5736728300154127016475733E-03
DATA .9812125319442814682424151,1.5133449185458729415321843E-03
DATA .9826956739534434545500623,1.4529238898903999088514186E-03
DATA .9841183499241993752558070,1.3924134618298698080684263E-03
DATA .9854804723183393394746898,1.3318173576538514742077422E-03
DATA .9867819573236030434401229,1.2711393059342868032230354E-03
DATA .9880227248588493179127030,1.2103830403000240185954168E-03
DATA .9892026985790008762180358,1.1495522992128106211745720E-03
DATA .9903218058797661642331360,1.0886508257457038100719739E-03
DATA .9913799779021414417532216,1.0276823673654497630797823E-03
DATA .9923771495366982433589461,9.6665067572141452408517998E-04
DATA .9933132594276647813762767,9.0555950644549534547199634E-04
DATA .9941882499768158116647105,8.4441261897085745246513294E-04
DATA .9950020673471962444031388,7.8321377638390310244817342E-04
DATA .9957546614667239851797758,7.2196674533702894531306739E-04
DATA .9964459860317571225964113,6.6067529607738161575318920E-04
DATA .9970759985107923450247946,5.9934320270831884264722444E-04
DATA .9976446601486403940368916,5.3797424394626279849881322E-04
DATA .9981519359718441275385541,4.7657220501003854013133061E-04
DATA .9985977947971757433952425,4.1514088233419025171823795E-04
DATA .9989822092480797992371162,3.5368409613102448239639817E-04
DATA .9993051557937030335325083,2.9220572803523217847223343E-04
DATA .9995666148626408000179844,2.3070985544248655269465132E-04
DATA .9997665712667108498622001,1.6920137336541413264932958E-04
DATA .9999050164693312795671775,1.0769038102183646230364914E-04
DATA .9999819727039624507107997,4.6263724177190118157440220E-05
END FUNCTION
EXTERNAL FUNCTION F(R,X,T)
LET F=EXP(X*COS(T))*COS(R*T)
END FUNCTION
EXTERNAL FUNCTION INTEGRAL2(R,Z) !'ガウス・ラゲール
FOR I=1 TO 6
READ X,W
LET S=S+FF(R,Z,X)*EXP(X)*W
NEXT I
LET INTEGRAL2=S
DATA .2228466041792606894643548,4.5896467394996359356828488E-01
DATA 1.1889321016726230307431509,4.1700083077212099411337757E-01
DATA 2.9927363260593140776913253,1.1337338207404497573870619E-01
DATA 5.7751435691045105018398304,1.0399197453149074898913303E-02
DATA 9.8374674183825899177155470,2.6101720281493205947924286E-04
DATA 15.9828739806017017825457916,8.9854790642962123882529205E-07
END FUNCTION
EXTERNAL FUNCTION FF(R,X,T)
LET FF=EXP(-X*COSH(T)-R*T)
END FUNCTION
EXTERNAL FUNCTION BESSELI(R,X) !'Rは実数
LET S=INTEGRAL(ABS(R),X)-SIN(ABS(R)*PI)/PI*INTEGRAL2(ABS(R),X)
IF R<0 THEN LET S=S*(-1)^INT(ABS(R))
LET BESSELI=S
END FUNCTION
|
|