hankelpt.m

这是在网上下的一个东东

%
%  Y=hankelpts(B)
%
%    B             The transform argument.
%
%  This routine returns a vector of points at which to evaluate the function.
%
%  Based on Walt Anderson's Fortran code, which was published in:
%
%  Anderson, W. L., 1979, Computer Program Numerical Integration of Related
%  Hankel Transforms of Orders 0 and 1 by Adaptive Digital Filtering.
%  Geophysic, 44(7):1287-1305.
%
function Y=hankelpts(B)
%
%  Points at which the function will be evaluated.  Note that these
%  are scaled by B before being used.
%
YBASE=[
   8.9170998013274418e-14
   9.8549193740052245e-14
   1.0891370292130841e-13
   1.2036825704856076e-13
   1.3302749714952345e-13
   1.4701812115404443e-13
   1.6248015192957209e-13
   1.7956833867707590e-13
   1.9845370571306282e-13
   2.1932526413842005e-13
   2.4239190352504162e-13
   2.6788448255287407e-13
   2.9605813952117967e-13
   3.2719484585839035e-13
   3.6160622818693727e-13
   3.9963668718722961e-13
   4.4166684447542096e-13
   4.8811735199247527e-13
   5.3945310203017795e-13
   5.9618788002944785e-13
   6.5888950671771899e-13
   7.2818552104963217e-13
   8.0476946082781583e-13
   8.8940780386232124e-13
   9.8294763913816720e-13
   1.0863251447666186e-12
   1.2005749575703847e-12
   1.3268405280766938e-12
   1.4663855645544968e-12
   1.6206066806315702e-12
   1.7910473730731204e-12
   1.9794134696161976e-12
   2.1875902014670364e-12
   2.4176610713286156e-12
   2.6719287057960000e-12
   2.9529379008172425e-12
   3.2635010908665675e-12
   3.6067264967338821e-12
   3.9860492336431492e-12
   4.4052656910401313e-12
   4.8685715281339743e-12
   5.3806036654647833e-12
   5.9464866927629101e-12
   6.5718841575654077e-12
   7.2630552479033658e-12
   8.0269174363595134e-12
   8.8711157124588670e-12
   9.8040990962934688e-12
   1.0835205199155280e-11
   1.1974753677488471e-11
   1.3234149515479672e-11
   1.4625997169973058e-11
   1.6164226720110948e-11
   1.7864233284247930e-11
   1.9743031099469828e-11
   2.1819423805797136e-11
   2.4114192639334462e-11
   2.6650304417866294e-11
   2.9453141400488783e-11
   3.2550755321790057e-11
   3.5974148143038491e-11
   3.9757582330231206e-11
   4.3938923764369768e-11
   4.8560020715924429e-11
   5.3667122676390673e-11
   5.9311343238745088e-11
   6.5549171659463761e-11
   7.2443038221987792e-11
   8.0061939059983487e-11
   8.8482126693838504e-11
   9.7787873191515278e-11
   1.0807231359173195e-10
   1.1943837803073370e-10
   1.3199982190169224e-10
   1.4588236435691521e-10
   1.6122494654737812e-10
   1.7818112219246310e-10
   1.9692059439719362e-10
   2.1763091409794871e-10
   2.4051935713527238e-10
   2.6581499874015355e-10
   2.9377100619593265e-10
   3.2466717262156566e-10
   3.5881271723520049e-10
   3.9654938012404428e-10
   4.3825484249401899e-10
   4.8434650663021337e-10
   5.3528567339924565e-10
   5.9158215910338558e-10
   6.5379939789346253e-10
   7.2256008080722362e-10
   7.9855238787053352e-10
   8.8253687563437822e-10
   9.7535408908045947e-10
   1.0779329740778886e-09
   1.1913001745856734e-09
   1.3165903076505279e-09
   1.4550573190356333e-09
   1.6080870331313015e-09
   1.7772110227512649e-09
   1.9641219376281762e-09
   2.1706904450210516e-09
   2.3989839519819520e-09
   2.6512872966606394e-09
   2.9301256157327411e-09
   3.2382896168163261e-09
   3.5788635088117363e-09
   3.9552558697009004e-09
   4.3712337607414382e-09
   4.8309604284818817e-09
   5.3390369719324459e-09
   5.9005483919104070e-09
   6.5211144834374116e-09
   7.2069460805369272e-09
   7.9649072163486862e-09
   8.8025838206794284e-09
   9.7283596425381269e-09
   1.0751500157513941e-08
   1.1882245299770154e-08
   1.3131911946747030e-08
   1.4513007182274982e-08
   1.6039353471673310e-08
   1.7726227001629019e-08
   1.9590510569407677e-08
   2.1650862551562963e-08
   2.3927903643240499e-08
   2.6444423237025737e-08
   2.9225607506844726e-08
   3.2299291479658129e-08
   3.5696237617766722e-08
   3.9450443699873716e-08
   4.3599483082281089e-08
   4.8184880745668262e-08
   5.3252528891055793e-08
   5.8853146244378082e-08
   6.5042785666539686e-08
   7.1883395149287240e-08
   7.9443437811532343e-08
   8.7798577101256825e-08
   9.7032434060731539e-08
   1.0723742423401342e-07
   1.1851568259277232e-07
   1.3098008573741623e-07
   1.4475538160404734e-07
   1.5997943798373573e-07
   1.7680462234971138e-07
   1.9539932680224871e-07
   2.1594965339340473e-07
   2.3866127669890702e-07
   2.6376150227843728e-07
   2.9150154162607256e-07
   3.2215902637935327e-07
   3.5604078695002664e-07
   3.9348592338593700e-07
   4.3486919919827997e-07
   4.8060479212078480e-07
   5.3115043933968360e-07
   5.8701201868132178e-07
   6.4874861160747570e-07
   7.1697809869053571e-07
   7.9238334356995170e-07
   8.7571902728105488e-07
   9.6781920135651654e-07
   1.0696056352944216e-06
   1.1820970419372223e-06
   1.3064192730922673e-06
   1.4438165874351013e-06
   1.5956641034684996e-06
   1.7634815621706371e-06
   1.9489485370736003e-06
   2.1539212439998211e-06
   2.3804511186939235e-06
   2.6308053482811660e-06
   2.9074895620382204e-06
   3.2132729085731429e-06
   3.5512157703953872e-06
   3.9247003932525885e-06
   4.3374647367828189e-06
   4.7936398852710157e-06
   5.2977913929290109e-06
   5.8549649774966190e-06
   6.4707370194807025e-06
   7.1512703724455683e-06
   7.9033760429228480e-06
   8.7345813572541237e-06
   9.6532052976029776e-06
   1.0668441761124589e-05
   1.1790451575578641e-05
   1.3030464192308714e-05
   1.4400890074365673e-05
   1.5915444904593194e-05
   1.7589286856791647e-05
   1.9439168303816350e-05
   2.1483603480955747e-05
   2.3743053782621043e-05
   2.6240132546858776e-05
   2.8999831377238596e-05
   3.2049770267221755e-05
   3.5420474030339066e-05
   3.9145677802784465e-05
   4.3262664675996812e-05
   4.7812638838370288e-05
   5.2841137960621061e-05
   5.8398488952101537e-05
   6.4540311649424626e-05
   7.1328075478483034e-05
   7.8829714661124184e-05
   8.7120308123675967e-05
   9.6282830912076270e-05
   1.0640898463402168e-04
   1.1760011523947924e-04
   1.2996822732501724e-04
   1.4363710511345378e-04
   1.5874355132796403e-04
   1.7543875635971472e-04
   1.9388981143211579e-04
   2.1428138090594562e-04
   2.3681755046234149e-04
   2.6172386966089197e-04
   2.8924960931543914e-04
   3.1967025628016629e-04
   3.5329027061462893e-04
   3.9044613272236350e-04
   4.3150971095986065e-04
   4.7689198342006658e-04
   5.2704715113927156e-04
   5.8247718389374341e-04
   6.4373684408196633e-04
   7.1143923897318674e-04
   7.8626195689103691e-04
   8.6895384874522254e-04
   9.6034252278312515e-04
   1.0613426275713101e-03
   1.1729650061058051e-03
   1.2963268126685603e-03
   1.4326626936829910e-03
   1.5833371444703617e-03
   1.7498581655775839e-03
   1.9338923553535471e-03
   2.1372815898255564e-03
   2.3620614568136901e-03
   2.6104816287778878e-03
   2.8850283782960702e-03
   3.1884494615157647e-03
   3.5237816186211822e-03
   3.8943809665496639e-03
   4.3039565881380203e-03
   4.7566076538702283e-03
   5.2568644477534125e-03
   5.8097337079228714e-03
   6.4207487357601564e-03
   7.0960247750331065e-03
   7.8423202153108801e-03
   8.6671042321983371e-03
   9.5786315413559676e-03
   1.0586025014468731e-02
   1.1699366984012178e-02
   1.2929800150624660e-02
   1.4289639103000504e-02
   1.5792493566432742e-02
   1.7453404613518231e-02
   1.9288995200267688e-02
   2.1317636534236600e-02
   2.3559631939745231e-02
   2.6037420060372587e-02
   2.8775799432443259e-02
   3.1802176677114019e-02
   3.5146840795050052e-02
   3.8843266308924096e-02
   4.2928448287690518e-02
   4.7443272605669898e-02
   5.2432925142121424e-02
   5.7947344016710048e-02
   6.4041719386992838e-02
   7.0777045810065886e-02
   7.8220732696592687e-02
   8.6447278966843177e-02
   9.5539018660927705e-02
   1.0558694496554391e-01
   1.1669162090437306e-01
   1.2896418580662142e-01
   1.4252746762678220e-01
   1.5751721224808804e-01
   1.7408344207293611e-01
   1.9239195749751564e-01
   2.1262599629790035e-01
   2.3498806753529980e-01
   2.5970197833480957e-01
   2.8701507382234348e-01
   3.1720071263778915e-01
   3.5056100280015512e-01
   3.8742982530616715e-01
   4.2817617572350458e-01
   4.7320785722246539e-01
   5.2297556200716211e-01
   5.7797738199458315e-01
   6.3876379388591276e-01
   7.0594316852237804e-01
   7.8018785966510817e-01
   8.6224093313756223e-01
   9.5292360367804285e-01
   1.0531434539328173e+00
   1.1639035178482904e+00
   1.2863123193718711e+00
   1.4215949669322265e+00
   1.5711054147362089e+00
   1.7363400135976372e+00
   1.9189524869191834e+00
   2.1207704817120212e+00
   2.3438138603014083e+00
   2.5903149157877352e+00
   2.8627407135861755e+00
   3.1638177826465683e+00
   3.4965594034715681e+00
   3.8642957660407120e+00
   4.2707072994710522e+00
   4.7198615069887930e+00
   5.2162536748687147e+00
   5.7648518627701284e+00
   6.3711466257477705e+00
   7.0412059655722290e+00
   7.7817360613311877e+00
   8.6001483871237632e+00
   9.5046338885843706e+00
   1.0504244960619703e+01
   1.1608986046819572e+01
   1.2829913767290970e+01
   1.4179247577028352e+01
   1.5670492062326328e+01
   1.7318572099218340e+01
   1.9139982226652432e+01
   2.1152951729381048e+01
   2.3377627082769912e+01
   2.5836273585494951e+01
   2.8553498198135060e+01
   3.1556495817904278e+01
   3.4875321454323611e+01
   3.8543191029858157e+01
   4.2596813816033411e+01
   4.7076759832163077e+01
   5.2027865883738443e+01
   5.7499684304247886e+01
   6.3546978891585546e+01
   7.0230273002547406e+01
   7.7616455290928698e+01
   8.5779449151653139e+01
   9.4800952570955843e+01
   1.0477125578728921e+02
   1.1579014494637693e+02
   1.2796790079449971e+02
   1.4142640240527066e+02
   1.5630034698636896e+02
   1.7273859797446769e+02
   1.9090567491054267e+02
   2.1098340000673559e+02
   2.3317271788416559e+02
   2.5769570669423729e+02
   2.8479780075142304e+02
   3.1475024692237560e+02
   3.4785281935573863e+02
   3.8443681972258412e+02
   4.2486839299489054e+02
   4.6955219194748827e+02
   5.1893542705903837e+02
   5.7351234234481581e+02
   6.3382916191693528e+02
   7.0048955677885772e+02
   7.7416068656769369e+02
   8.5557987671209173e+02
   9.4556199783295187e+02
   1.0450076212424869e+03
   1.1549120321646080e+03
   1.2763751908839718e+03
   1.4106127415182191e+03
   1.5589681785928965e+03
   1.7229262931862318e+03
   1.9041280332173003e+03
   2.1043869266043412e+03
   2.3257072316617105e+03
   2.5703039963907459e+03
   2.8406252274246667e+03
   3.1393763905017645e+03
   3.4695474876758481e+03
   3.8344429822617740e+03
   4.2377148710149695e+03
   4.6833992345424385e+03
   5.1759566317540530e+03
   5.7203167426353639e+03
   6.3219277061418234e+03
   6.9868106470046323e+03
   7.7216199371708199e+03
   8.5337097949943000e+03
   9.4312078887249972e+03
   1.0423096680944496e+04
   1.1519303328070666e+04