View Javadoc
1   /*
2    * Licensed to the Apache Software Foundation (ASF) under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * The ASF licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *      https://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  
18  /*
19   * This is not the original file distributed by the Apache Software Foundation
20   * It has been modified by the Hipparchus project
21   */
22  package org.hipparchus.special;
23  
24  import org.hipparchus.exception.MathIllegalArgumentException;
25  import org.junit.Assert;
26  import org.junit.Test;
27  
28  /**
29   * @version $Id$
30   */
31  public class BesselJTest {
32  
33      /**
34       * Reference data for the {@link BesselJ#value(double, double)} function. This data
35       * was generated with the following <a href="http://www.r-project.org/">R</a> script.
36       *
37       * <pre>
38       * smallxs = 10 ** (-(seq(0,8,.5)))
39       * medxs = seq(1,20)
40       * near.eight = 8 + seq(-.5,.5,.1)
41       * largexs = c(10,30,100,300,1000)
42       * xs = unique(sort(c(smallxs, medxs, near.eight,largexs)))
43       *
44       * for (n in c(0:15, 30, 100)) {
45       * for (x in xs) {
46       * val = format( besselJ(x,n), digits=20 )
47       * s = paste("{", n, ",", x, ",", val, "},")
48       * print(s)
49       * }
50       * }
51       * </pre>
52       */
53      private static final double[][] BESSEL_J_REF = {
54          { 0 , 1e-08 , 1 },
55          { 0 , 3.16227766016838e-08 , 0.99999999999999977796 },
56          { 0 , 1e-07 , 0.99999999999999744649 },
57          { 0 , 3.16227766016838e-07 , 0.99999999999997501998 },
58          { 0 , 1e-06 , 0.99999999999974997777 },
59          { 0 , 3.16227766016838e-06 , 0.99999999999749999979 },
60          { 0 , 1e-05 , 0.99999999997499999793 },
61          { 0 , 3.16227766016838e-05 , 0.99999999974999997931 },
62          { 0 , 1e-04 , 0.99999999750000001519 },
63          { 0 , 0.000316227766016838 , 0.99999997500000015194 },
64          { 0 , 0.001 , 0.9999997500000156192 },
65          { 0 , 0.00316227766016838 , 0.99999750000156251151 },
66          { 0 , 0.01 , 0.99997500015624951608 },
67          { 0 , 0.0316227766016838 , 0.9997500156245660019 },
68          { 0 , 0.1 , 0.99750156206604001508 },
69          { 0 , 0.316227766016838 , 0.97515581664971295872 },
70          { 0 , 1 , 0.76519768655796649437 },
71          { 0 , 2 , 0.22389077914123567403 },
72          { 0 , 3 , -0.26005195490193344643 },
73          { 0 , 4 , -0.39714980986384734729 },
74          { 0 , 5 , -0.177596771314338292 },
75          { 0 , 6 , 0.15064525725099692233 },
76          { 0 , 7 , 0.3000792705195556298 },
77          { 0 , 7.5 , 0.26633965788037838873 },
78          { 0 , 7.6 , 0.25160183384997636402 },
79          { 0 , 7.7 , 0.23455913958646437689 },
80          { 0 , 7.8 , 0.21540780774626291927 },
81          { 0 , 7.9 , 0.19436184484127824734 },
82          { 0 , 8 , 0.17165080713755390129 },
83          { 0 , 8.1 , 0.14751745404437766052 },
84          { 0 , 8.2 , 0.12221530178413773926 },
85          { 0 , 8.3 , 0.09600610089501022959 },
86          { 0 , 8.4 , 0.069157261656985186127 },
87          { 0 , 8.5 , 0.041939251842934503756 },
88          { 0 , 9 , -0.090333611182876139001 },
89          { 0 , 10 , -0.24593576445134832098 },
90          { 0 , 11 , -0.17119030040719609986 },
91          { 0 , 12 , 0.0476893107968335353 },
92          { 0 , 13 , 0.20692610237706782206 },
93          { 0 , 14 , 0.17107347611045864433 },
94          { 0 , 15 , -0.014224472826780772475 },
95          { 0 , 16 , -0.17489907398362919411 },
96          { 0 , 17 , -0.16985425215118354902 },
97          { 0 , 18 , -0.013355805721984111492 },
98          { 0 , 19 , 0.14662943965965119508 },
99          { 0 , 20 , 0.16702466434058313438 },
100         { 0 , 30 , -0.086367983581040239094 },
101         { 0 , 100 , 0.019985850304223118368 },
102         { 0 , 300 , -0.033298554876305661021 },
103         { 0 , 1000 , 0.024786686152420175921 },
104         { 1 , 1e-08 , 5.0000000000000001046e-09 },
105         { 1 , 3.16227766016838e-08 , 1.5811388300841892647e-08 },
106         { 1 , 1e-07 , 4.999999999999993818e-08 },
107         { 1 , 3.16227766016838e-07 , 1.5811388300841697432e-07 },
108         { 1 , 1e-06 , 4.9999999999993750869e-07 },
109         { 1 , 3.16227766016838e-06 , 1.5811388300822132559e-06 },
110         { 1 , 1e-05 , 4.9999999999375003889e-06 },
111         { 1 , 3.16227766016838e-05 , 1.5811388298865475016e-05 },
112         { 1 , 1e-04 , 4.9999999937500002646e-05 },
113         { 1 , 0.000316227766016838 , 0.00015811388103199544184 },
114         { 1 , 0.001 , 0.00049999993750000253697 },
115         { 1 , 0.00316227766016838 , 0.0015811368536614756226 },
116         { 1 , 0.01 , 0.0049999375002604158624 },
117         { 1 , 0.0316227766016838 , 0.015809411959653556917 },
118         { 1 , 0.1 , 0.049937526036241998428 },
119         { 1 , 0.316227766016838 , 0.15614567743386048582 },
120         { 1 , 1 , 0.44005058574493355339 },
121         { 1 , 2 , 0.57672480775687340326 },
122         { 1 , 3 , 0.33905895852593642692 },
123         { 1 , 4 , -0.066043328023549133232 },
124         { 1 , 5 , -0.32757913759146523036 },
125         { 1 , 6 , -0.27668385812756562947 },
126         { 1 , 7 , -0.0046828234823458325664 },
127         { 1 , 7.5 , 0.13524842757970551022 },
128         { 1 , 7.6 , 0.15921376839635667522 },
129         { 1 , 7.7 , 0.18131271532458800855 },
130         { 1 , 7.8 , 0.20135687275589611578 },
131         { 1 , 7.9 , 0.21917939992175122788 },
132         { 1 , 8 , 0.23463634685391462908 },
133         { 1 , 8.1 , 0.2476077669815928417 },
134         { 1 , 8.2 , 0.25799859764868082745 },
135         { 1 , 8.3 , 0.26573930204186430037 },
136         { 1 , 8.4 , 0.27078626827683538458 },
137         { 1 , 8.5 , 0.27312196367405372488 },
138         { 1 , 9 , 0.24531178657332528004 },
139         { 1 , 10 , 0.043472746168861438332 },
140         { 1 , 11 , -0.17678529895672151495 },
141         { 1 , 12 , -0.2234471044906276016 },
142         { 1 , 13 , -0.070318052121778371055 },
143         { 1 , 14 , 0.13337515469879324126 },
144         { 1 , 15 , 0.20510403861352277666 },
145         { 1 , 16 , 0.090397175661304188243 },
146         { 1 , 17 , -0.097668492757780639435 },
147         { 1 , 18 , -0.18799488548806958521 },
148         { 1 , 19 , -0.10570143114240927729 },
149         { 1 , 20 , 0.066833124175850036619 },
150         { 1 , 30 , -0.1187510626166229516 },
151         { 1 , 100 , -0.077145352014112156258 },
152         { 1 , 300 , -0.031887431377499955709 },
153         { 1 , 1000 , 0.0047283119070895248195 },
154         { 2 , 1e-08 , 1.2499999999999999739e-17 },
155         { 2 , 3.16227766016838e-08 , 1.2499999999999998506e-16 },
156         { 2 , 1e-07 , 1.2499999999999988152e-15 },
157         { 2 , 3.16227766016838e-07 , 1.2499999999999894672e-14 },
158         { 2 , 1e-06 , 1.249999999999895719e-13 },
159         { 2 , 3.16227766016838e-06 , 1.2499999999989582746e-12 },
160         { 2 , 1e-05 , 1.2499999999895835475e-11 },
161         { 2 , 3.16227766016838e-05 , 1.2499999998958334818e-10 },
162         { 2 , 1e-04 , 1.2499999989583335487e-09 },
163         { 2 , 0.000316227766016838 , 1.2499999895833333493e-08 },
164         { 2 , 0.001 , 1.2499998958333367811e-07 },
165         { 2 , 0.00316227766016838 , 1.2499989583336589057e-06 },
166         { 2 , 0.01 , 1.2499895833658854395e-05 },
167         { 2 , 0.0316227766016838 , 0.00012498958365884872863 },
168         { 2 , 0.1 , 0.0012489586587999190141 },
169         { 2 , 0.316227766016838 , 0.012396158312196680837 },
170         { 2 , 1 , 0.11490348493190047363 },
171         { 2 , 2 , 0.35283402861563772923 },
172         { 2 , 3 , 0.48609126058589108288 },
173         { 2 , 4 , 0.36412814585207281537 },
174         { 2 , 5 , 0.046565116277752213736 },
175         { 2 , 6 , -0.24287320996018546548 },
176         { 2 , 7 , -0.3014172200859401296 },
177         { 2 , 7.5 , -0.23027341052579025638 },
178         { 2 , 7.6 , -0.20970347374567196996 },
179         { 2 , 7.7 , -0.18746492781384410664 },
180         { 2 , 7.8 , -0.16377784037295622932 },
181         { 2 , 7.9 , -0.13887338916488553564 },
182         { 2 , 8 , -0.11299172042407523708 },
183         { 2 , 8.1 , -0.086379733802009056598 },
184         { 2 , 8.2 , -0.059288814552752158726 },
185         { 2 , 8.3 , -0.031972534137934507936 },
186         { 2 , 8.4 , -0.0046843406386910518141 },
187         { 2 , 8.5 , 0.022324739609784025052 },
188         { 2 , 9 , 0.14484734153250397592 },
189         { 2 , 10 , 0.25463031368512062391 },
190         { 2 , 11 , 0.13904751877870125121 },
191         { 2 , 12 , -0.084930494878604809172 },
192         { 2 , 13 , -0.21774426424195678087 },
193         { 2 , 14 , -0.15201988258205964555 },
194         { 2 , 15 , 0.041571677975250471981 },
195         { 2 , 16 , 0.18619872094129222284 },
196         { 2 , 17 , 0.15836384123850347216 },
197         { 2 , 18 , -0.0075325148878013998069 },
198         { 2 , 19 , -0.15775590609569428713 },
199         { 2 , 20 , -0.16034135192299814321 },
200         { 2 , 30 , 0.07845124607326538213 },
201         { 2 , 100 , -0.021528757344505360799 },
202         { 2 , 300 , 0.033085972000455661501 },
203         { 2 , 1000 , -0.024777229528605997089 },
204         { 3 , 1e-08 , 2.0833333333333334614e-26 },
205         { 3 , 3.16227766016838e-08 , 6.5880784586841223417e-25 },
206         { 3 , 1e-07 , 2.0833333333333317693e-23 },
207         { 3 , 3.16227766016838e-07 , 6.5880784586840819929e-22 },
208         { 3 , 1e-06 , 2.0833333333332027799e-20 },
209         { 3 , 3.16227766016838e-06 , 6.5880784586800051603e-19 },
210         { 3 , 1e-05 , 2.0833333333203129762e-17 },
211         { 3 , 3.16227766016838e-05 , 6.5880784582723696076e-16 },
212         { 3 , 1e-04 , 2.083333332031250315e-14 },
213         { 3 , 0.000316227766016838 , 6.5880784175086335665e-13 },
214         { 3 , 0.001 , 2.0833332031250032117e-11 },
215         { 3 , 0.00316227766016838 , 6.5880743411361163135e-10 },
216         { 3 , 0.01 , 2.083320312532551971e-08 },
217         { 3 , 0.0316227766016838 , 6.5876667140741846331e-07 },
218         { 3 , 0.1 , 2.0820315754756265453e-05 },
219         { 3 , 0.316227766016838 , 0.00065470057642003857534 },
220         { 3 , 1 , 0.019563353982668403586 },
221         { 3 , 2 , 0.1289432494744020552 },
222         { 3 , 3 , 0.30906272225525166508 },
223         { 3 , 4 , 0.43017147387562193472 },
224         { 3 , 5 , 0.36483123061366695694 },
225         { 3 , 6 , 0.11476838482077529602 },
226         { 3 , 7 , -0.16755558799533423753 },
227         { 3 , 7.5 , -0.25806091319346030621 },
228         { 3 , 7.6 , -0.26958401773618401176 },
229         { 3 , 7.7 , -0.27869709340970183487 },
230         { 3 , 7.8 , -0.28534550884459158882 },
231         { 3 , 7.9 , -0.2894950400052375139 },
232         { 3 , 8 , -0.29113220706595221987 },
233         { 3 , 8.1 , -0.29026442564925164502 },
234         { 3 , 8.2 , -0.28691997060124291297 },
235         { 3 , 8.3 , -0.28114775222882071315 },
236         { 3 , 8.4 , -0.27301690667621203445 },
237         { 3 , 8.5 , -0.26261620385768480457 },
238         { 3 , 9 , -0.1809351903366568648 },
239         { 3 , 10 , 0.058379379305186815396 },
240         { 3 , 11 , 0.22734803305806741691 },
241         { 3 , 12 , 0.19513693953109267909 },
242         { 3 , 13 , 0.0033198169704070513292 },
243         { 3 , 14 , -0.17680940686509599713 },
244         { 3 , 15 , -0.19401825782012263599 },
245         { 3 , 16 , -0.043847495425981139472 },
246         { 3 , 17 , 0.13493057304919323092 },
247         { 3 , 18 , 0.18632099329078039007 },
248         { 3 , 19 , 0.072489661438052577225 },
249         { 3 , 20 , -0.098901394560449676363 },
250         { 3 , 30 , 0.12921122875972501642 },
251         { 3 , 100 , 0.076284201720331942798 },
252         { 3 , 300 , 0.032328577670839367397 },
253         { 3 , 1000 , -0.0048274208252039483777 },
254         { 4 , 1e-08 , 2.6041666666666666342e-35 },
255         { 4 , 3.16227766016838e-08 , 2.6041666666666659714e-33 },
256         { 4 , 1e-07 , 2.6041666666666649861e-31 },
257         { 4 , 3.16227766016838e-07 , 2.6041666666666531276e-29 },
258         { 4 , 1e-06 , 2.6041666666665358894e-27 },
259         { 4 , 3.16227766016838e-06 , 2.6041666666653639536e-25 },
260         { 4 , 1e-05 , 2.6041666666536465525e-23 },
261         { 4 , 3.16227766016838e-05 , 2.604166666536458817e-21 },
262         { 4 , 1e-04 , 2.6041666653645840559e-19 },
263         { 4 , 0.000316227766016838 , 2.6041666536458338462e-17 },
264         { 4 , 0.001 , 2.6041665364583368871e-15 },
265         { 4 , 0.00316227766016838 , 2.604165364583604802e-13 },
266         { 4 , 0.01 , 2.6041536458604605458e-11 },
267         { 4 , 0.0316227766016838 , 2.6040364610459735901e-09 },
268         { 4 , 0.1 , 2.6028648545684040871e-07 },
269         { 4 , 0.316227766016838 , 2.5911729278009268374e-05 },
270         { 4 , 1 , 0.002476638964109955255 },
271         { 4 , 2 , 0.033995719807568429427 },
272         { 4 , 3 , 0.13203418392461221953 },
273         { 4 , 4 , 0.28112906496136008672 },
274         { 4 , 5 , 0.39123236045864817623 },
275         { 4 , 6 , 0.35764159478096080313 },
276         { 4 , 7 , 0.15779814466136793394 },
277         { 4 , 7.5 , 0.02382467997102201071 },
278         { 4 , 7.6 , -0.0031260139407891210546 },
279         { 4 , 7.7 , -0.029701638479430046702 },
280         { 4 , 7.8 , -0.055718704892114230554 },
281         { 4 , 7.9 , -0.080996261472003727722 },
282         { 4 , 8 , -0.10535743487538892782 },
283         { 4 , 8.1 , -0.12863095186410331006 },
284         { 4 , 8.2 , -0.15065262735059631316 },
285         { 4 , 8.3 , -0.17126680482265874139 },
286         { 4 , 8.4 , -0.19032773555860327264 },
287         { 4 , 8.5 , -0.2077008835093262229 },
288         { 4 , 9 , -0.26547080175694187654 },
289         { 4 , 10 , -0.21960268610200855965 },
290         { 4 , 11 , -0.015039500747028132846 },
291         { 4 , 12 , 0.18249896464415113484 },
292         { 4 , 13 , 0.21927648745906774819 },
293         { 4 , 14 , 0.076244422497018474183 },
294         { 4 , 15 , -0.11917898110329952499 },
295         { 4 , 16 , -0.20264153172603513453 },
296         { 4 , 17 , -0.11074128604467056713 },
297         { 4 , 18 , 0.069639512651394869236 },
298         { 4 , 19 , 0.18064737812876355272 },
299         { 4 , 20 , 0.13067093355486322781 },
300         { 4 , 30 , -0.05260900032132037607 },
301         { 4 , 100 , 0.026105809447725277644 },
302         { 4 , 300 , -0.032439400447038871378 },
303         { 4 , 1000 , 0.024748265003654772859 },
304         { 5 , 1e-08 , 2.6041666666666666817e-44 },
305         { 5 , 3.16227766016838e-08 , 8.2350980733551520153e-42 },
306         { 5 , 1e-07 , 2.6041666666666648818e-39 },
307         { 5 , 3.16227766016838e-07 , 8.2350980733551185479e-37 },
308         { 5 , 1e-06 , 2.6041666666665576923e-34 },
309         { 5 , 3.16227766016838e-06 , 8.2350980733517218878e-32 },
310         { 5 , 1e-05 , 2.6041666666558171668e-29 },
311         { 5 , 3.16227766016838e-05 , 8.2350980730120282499e-27 },
312         { 5 , 1e-04 , 2.6041666655815978351e-24 },
313         { 5 , 0.000316227766016838 , 8.2350980390422474771e-22 },
314         { 5 , 0.001 , 2.6041665581597245947e-19 },
315         { 5 , 0.00316227766016838 , 8.2350946420649040367e-17 },
316         { 5 , 0.01 , 2.6041558159915982186e-14 },
317         { 5 , 0.0316227766016838 , 8.2347549503960048977e-12 },
318         { 5 , 0.1 , 2.6030817909644421178e-09 },
319         { 5 , 0.316227766016838 , 8.2008463739855235578e-07 },
320         { 5 , 1 , 0.00024975773021123438876 },
321         { 5 , 2 , 0.0070396297558716850601 },
322         { 5 , 3 , 0.043028434877047584683 },
323         { 5 , 4 , 0.13208665604709826646 },
324         { 5 , 5 , 0.26114054612017006951 },
325         { 5 , 6 , 0.36208707488717239986 },
326         { 5 , 7 , 0.34789632475118331678 },
327         { 5 , 7.5 , 0.28347390516255044357 },
328         { 5 , 7.6 , 0.26629347674587972028 },
329         { 5 , 7.7 , 0.2478382482362680439 },
330         { 5 , 7.8 , 0.22819811921165392143 },
331         { 5 , 7.9 , 0.20747350940067682545 },
332         { 5 , 8 , 0.18577477219056331981 },
333         { 5 , 8.1 , 0.16322151022791506203 },
334         { 5 , 8.2 , 0.13994179757627084326 },
335         { 5 , 8.3 , 0.11607131384553516507 },
336         { 5 , 8.4 , 0.091752396620399440108 },
337         { 5 , 8.5 , 0.067133019378318919967 },
338         { 5 , 9 , -0.055038855669513706004 },
339         { 5 , 10 , -0.23406152818679362704 },
340         { 5 , 11 , -0.23828585178317879256 },
341         { 5 , 12 , -0.073470963101658584571 },
342         { 5 , 13 , 0.13161955992748081146 },
343         { 5 , 14 , 0.22037764829196368477 },
344         { 5 , 15 , 0.1304561345650295523 },
345         { 5 , 16 , -0.057473270437036434732 },
346         { 5 , 17 , -0.18704411942315585238 },
347         { 5 , 18 , -0.15537009877904933708 },
348         { 5 , 19 , 0.0035723925109004857348 },
349         { 5 , 20 , 0.15116976798239498136 },
350         { 5 , 30 , -0.14324029551207712041 },
351         { 5 , 100 , -0.074195736964513925304 },
352         { 5 , 300 , -0.03319362834942707341 },
353         { 5 , 1000 , 0.0050254069452331864842 },
354         { 6 , 1e-08 , 2.170138888888889163e-53 },
355         { 6 , 3.16227766016838e-08 , 2.1701388888888880947e-50 },
356         { 6 , 1e-07 , 2.1701388888888875174e-47 },
357         { 6 , 3.16227766016838e-07 , 2.170138888888880604e-44 },
358         { 6 , 1e-06 , 2.1701388888888106952e-41 },
359         { 6 , 3.16227766016838e-06 , 2.1701388888881133808e-38 },
360         { 6 , 1e-05 , 2.1701388888811393588e-35 },
361         { 6 , 3.16227766016838e-05 , 2.1701388888113848269e-32 },
362         { 6 , 1e-04 , 2.1701388881138396044e-29 },
363         { 6 , 0.000316227766016838 , 2.1701388811383932341e-26 },
364         { 6 , 0.001 , 2.1701388113839301844e-23 },
365         { 6 , 0.00316227766016838 , 2.1701381138394068717e-20 },
366         { 6 , 0.01 , 2.1701311384049674283e-17 },
367         { 6 , 0.0316227766016838 , 2.1700613851395740978e-14 },
368         { 6 , 0.1 , 2.1693639603760032489e-11 },
369         { 6 , 0.316227766016838 , 2.1624004918010960028e-08 },
370         { 6 , 1 , 2.0938338002389272967e-05 },
371         { 6 , 2 , 0.0012024289717899932714 },
372         { 6 , 3 , 0.011393932332213070266 },
373         { 6 , 4 , 0.049087575156385579445 },
374         { 6 , 5 , 0.13104873178169201831 },
375         { 6 , 6 , 0.24583686336432652997 },
376         { 6 , 7 , 0.33919660498317966146 },
377         { 6 , 7.5 , 0.35414052691237862813 },
378         { 6 , 7.6 , 0.35351216755378872536 },
379         { 6 , 7.7 , 0.35156949333172621275 },
380         { 6 , 7.8 , 0.3482803961891063893 },
381         { 6 , 7.9 , 0.34362095691589844559 },
382         { 6 , 8 , 0.33757590011359311921 },
383         { 6 , 8.1 , 0.33013898918251699532 },
384         { 6 , 8.2 , 0.32131335610214611931 },
385         { 6 , 8.3 , 0.3111117612630625584 },
386         { 6 , 8.4 , 0.29955677915431688785 },
387         { 6 , 8.5 , 0.28668090630734854862 },
388         { 6 , 9 , 0.20431651767970440692 },
389         { 6 , 10 , -0.014458842084785107282 },
390         { 6 , 11 , -0.20158400087404348966 },
391         { 6 , 12 , -0.24372476722886662892 },
392         { 6 , 13 , -0.1180306721302363665 },
393         { 6 , 14 , 0.081168183425812737153 },
394         { 6 , 15 , 0.20614973747998591169 },
395         { 6 , 16 , 0.16672073770288736716 },
396         { 6 , 17 , 0.00071533344281418307069 },
397         { 6 , 18 , -0.15595623419531115528 },
398         { 6 , 19 , -0.17876717154407903432 },
399         { 6 , 20 , -0.055086049563665764883 },
400         { 6 , 30 , 0.0048622351506280026001 },
401         { 6 , 100 , -0.033525383144176669481 },
402         { 6 , 300 , 0.031332946168724638836 },
403         { 6 , 1000 , -0.024698010934202440508 },
404         { 7 , 1e-08 , 1.5500992063492066701e-62 },
405         { 7 , 3.16227766016838e-08 , 4.9018440912828279875e-59 },
406         { 7 , 1e-07 , 1.5500992063492053031e-55 },
407         { 7 , 3.16227766016838e-07 , 4.9018440912828133382e-52 },
408         { 7 , 1e-06 , 1.55009920634915736e-48 },
409         { 7 , 3.16227766016838e-06 , 4.9018440912812964979e-45 },
410         { 7 , 1e-05 , 1.550099206344363137e-41 },
411         { 7 , 3.16227766016838e-05 , 4.9018440911296494971e-38 },
412         { 7 , 1e-04 , 1.5500992058648010339e-34 },
413         { 7 , 0.000316227766016838 , 4.9018440759645687969e-31 },
414         { 7 , 0.001 , 1.5500991579086071003e-27 },
415         { 7 , 0.00316227766016838 , 4.9018425594567649302e-24 },
416         { 7 , 0.01 , 1.550094362295914728e-20 },
417         { 7 , 0.0316227766016838 , 4.9016909107824929132e-17 },
418         { 7 , 0.1 , 1.5496148676202282287e-13 },
419         { 7 , 0.316227766016838 , 4.8865470861431505644e-10 },
420         { 7 , 1 , 1.5023258174368078499e-06 },
421         { 7 , 2 , 0.00017494407486827416175 },
422         { 7 , 3 , 0.0025472944518046929108 },
423         { 7 , 4 , 0.015176069422058449318 },
424         { 7 , 5 , 0.053376410155890716136 },
425         { 7 , 6 , 0.12958665184148068783 },
426         { 7 , 7 , 0.23358356950569605925 },
427         { 7 , 7.5 , 0.28315093789725531703 },
428         { 7 , 7.6 , 0.29188362991799726709 },
429         { 7 , 7.7 , 0.30006226085213638655 },
430         { 7 , 7.8 , 0.30761787492543296585 },
431         { 7 , 7.9 , 0.31448237452220684229 },
432         { 7 , 8 , 0.32058907797982633125 },
433         { 7 , 8.1 , 0.3258732885609990082 },
434         { 7 , 8.2 , 0.33027286989028453723 },
435         { 7 , 8.3 , 0.33372882292033839713 },
436         { 7 , 8.4 , 0.33618585931433897507 },
437         { 7 , 8.5 , 0.33759296599676130723 },
438         { 7 , 9 , 0.32746087924245292911 },
439         { 7 , 10 , 0.21671091768505151842 },
440         { 7 , 11 , 0.018376032647858614455 },
441         { 7 , 12 , -0.17025380412720803047 },
442         { 7 , 13 , -0.24057094958616048741 },
443         { 7 , 14 , -0.15080491964126707671 },
444         { 7 , 15 , 0.034463655418959161791 },
445         { 7 , 16 , 0.18251382371420196704 },
446         { 7 , 17 , 0.1875490606769070201 },
447         { 7 , 18 , 0.051399275982175231248 },
448         { 7 , 19 , -0.11647797453873988405 },
449         { 7 , 20 , -0.18422139772059445417 },
450         { 7 , 30 , 0.1451851895723283159 },
451         { 7 , 100 , 0.070172690987212724134 },
452         { 7 , 300 , 0.034446946196176060628 },
453         { 7 , 1000 , -0.0053217830764436153956 },
454         { 8 , 1e-08 , 9.6881200396825412359e-72 },
455         { 8 , 3.16227766016838e-08 , 9.6881200396825359082e-68 },
456         { 8 , 1e-07 , 9.6881200396825335915e-64 },
457         { 8 , 3.16227766016838e-07 , 9.6881200396825091166e-60 },
458         { 8 , 1e-06 , 9.6881200396822669073e-56 },
459         { 8 , 3.16227766016838e-06 , 9.6881200396798449492e-52 },
460         { 8 , 1e-05 , 9.6881200396556345156e-48 },
461         { 8 , 3.16227766016838e-05 , 9.6881200394134308774e-44 },
462         { 8 , 1e-04 , 9.6881200369913995322e-40 },
463         { 8 , 0.000316227766016838 , 9.688120012771098157e-36 },
464         { 8 , 0.001 , 9.6881197705681010442e-32 },
465         { 8 , 0.00316227766016838 , 9.688117348538421731e-28 },
466         { 8 , 0.01 , 9.6880931282716245736e-24 },
467         { 8 , 0.0316227766016838 , 9.6878509286008909493e-20 },
468         { 8 , 0.1 , 9.6854292315946525669e-16 },
469         { 8 , 0.316227766016838 , 9.6612422089625085973e-12 },
470         { 8 , 1 , 9.4223441726045005392e-08 },
471         { 8 , 2 , 2.2179552287925904881e-05 },
472         { 8 , 3 , 0.00049344177620883479096 },
473         { 8 , 4 , 0.0040286678208190035769 },
474         { 8 , 5 , 0.018405216654802002835 },
475         { 8 , 6 , 0.056531990932461785582 },
476         { 8 , 7 , 0.12797053402821254031 },
477         { 8 , 7.5 , 0.17440789049583127479 },
478         { 8 , 7.6 , 0.18416820334778524759 },
479         { 8 , 7.7 , 0.19399825367215817185 },
480         { 8 , 7.8 , 0.20385425111295268907 },
481         { 8 , 7.9 , 0.21368958021206302389 },
482         { 8 , 8 , 0.22345498635110294661 },
483         { 8 , 8.1 , 0.23309879351550599758 },
484         { 8 , 8.2 , 0.24256715346663235144 },
485         { 8 , 8.3 , 0.25180432559052018382 },
486         { 8 , 8.4 , 0.26075298636958132992 },
487         { 8 , 8.5 , 0.26935456709908189854 },
488         { 8 , 9 , 0.30506707225300011554 },
489         { 8 , 10 , 0.31785412684385727644 },
490         { 8 , 11 , 0.22497167878949989039 },
491         { 8 , 12 , 0.045095329080457241533 },
492         { 8 , 13 , -0.14104573511639803551 },
493         { 8 , 14 , -0.23197310306707982774 },
494         { 8 , 15 , -0.17398365908895732646 },
495         { 8 , 16 , -0.0070211419529606520704 },
496         { 8 , 17 , 0.1537368341734622057 },
497         { 8 , 18 , 0.19593344884811411677 },
498         { 8 , 19 , 0.092941295568165452345 },
499         { 8 , 20 , -0.073868928840750344711 },
500         { 8 , 30 , 0.062890853316458550371 },
501         { 8 , 100 , 0.043349559882386451415 },
502         { 8 , 300 , -0.029725422012903089664 },
503         { 8 , 1000 , 0.02462350597113223058 } ,
504         { 9 , 1e-08 , 5.382288910934745386e-81 },
505         { 9 , 3.16227766016838e-08 , 1.702029198362092975e-76 },
506         { 9 , 1e-07 , 5.3822889109347404393e-72 },
507         { 9 , 3.16227766016838e-07 , 1.7020291983620889989e-67 },
508         { 9 , 1e-06 , 5.3822889109346077433e-63 },
509         { 9 , 3.16227766016838e-06 , 1.7020291983616675345e-58 },
510         { 9 , 1e-05 , 5.3822889109212923968e-54 },
511         { 9 , 3.16227766016838e-05 , 1.7020291983195439884e-49 },
512         { 9 , 1e-04 , 5.3822889095891748719e-45 },
513         { 9 , 0.000316227766016838 , 1.7020291941070210695e-40 },
514         { 9 , 0.001 , 5.382288776377523226e-36 },
515         { 9 , 0.00316227766016838 , 1.7020287728548427703e-31 },
516         { 9 , 0.01 , 5.3822754552277587118e-27 },
517         { 9 , 0.0316227766016838 , 1.7019866481156611902e-22 },
518         { 9 , 0.1 , 5.3809434916023306372e-18 },
519         { 9 , 0.316227766016838 , 1.6977789573201714453e-13 },
520         { 9 , 1 , 5.2492501799118757129e-09 },
521         { 9 , 2 , 2.492343435133064173e-06 },
522         { 9 , 3 , 8.4395021309091773631e-05 },
523         { 9 , 4 , 0.00093860186121756401367 },
524         { 9 , 5 , 0.005520283139475687037 },
525         { 9 , 6 , 0.021165323978417364265 },
526         { 9 , 7 , 0.058920508273075426764 },
527         { 9 , 7.5 , 0.088919228493851462658 },
528         { 9 , 7.6 , 0.095838903445761125521 },
529         { 9 , 7.7 , 0.10305099353156887965 },
530         { 9 , 7.8 , 0.11054469146011103309 },
531         { 9 , 7.9 , 0.11830664869209804591 },
532         { 9 , 8 , 0.12632089472237958971 },
533         { 9 , 8.1 , 0.13456877270419806414 },
534         { 9 , 8.2 , 0.14302889297143717151 },
535         { 9 , 8.3 , 0.15167710592885710885 },
536         { 9 , 8.4 , 0.16048649567533976312 },
537         { 9 , 8.5 , 0.16942739560151048872 },
538         { 9 , 9 , 0.2148805825406584491 },
539         { 9 , 10 , 0.29185568526512006837 },
540         { 9 , 11 , 0.30885550013686852155 },
541         { 9 , 12 , 0.23038090956781773211 },
542         { 9 , 13 , 0.066976198673670619965 },
543         { 9 , 14 , -0.11430719814968128001 },
544         { 9 , 15 , -0.22004622511384699934 },
545         { 9 , 16 , -0.18953496566716260263 },
546         { 9 , 17 , -0.042855569690119083015 },
547         { 9 , 18 , 0.12276378966059287023 },
548         { 9 , 19 , 0.19474432870140553908 },
549         { 9 , 20 , 0.12512625464799415065 },
550         { 9 , 30 , -0.11164340113688375755 },
551         { 9 , 100 , -0.063236761406030891908 },
552         { 9 , 300 , -0.036032302036864222172 },
553         { 9 , 1000 , 0.0057157591719817308837 },
554         { 10 , 1e-08 , 2.6911444554673727331e-90 },
555         { 10 , 3.16227766016838e-08 , 2.6911444554673710334e-85 },
556         { 10 , 1e-07 , 2.6911444554673703522e-80 },
557         { 10 , 3.16227766016838e-07 , 2.6911444554673646193e-75 },
558         { 10 , 1e-06 , 2.6911444554673096152e-70 },
559         { 10 , 3.16227766016838e-06 , 2.6911444554667591352e-65 },
560         { 10 , 1e-05 , 2.6911444554612582946e-60 },
561         { 10 , 3.16227766016838e-05 , 2.6911444554062112799e-55 },
562         { 10 , 1e-04 , 2.6911444548557502587e-50 },
563         { 10 , 0.000316227766016838 , 2.691144449351135505e-45 },
564         { 10 , 0.001 , 2.6911443943049990395e-40 },
565         { 10 , 0.00316227766016838 , 2.6911438438436965201e-35 },
566         { 10 , 0.01 , 2.6911383392363445476e-30 },
567         { 10 , 0.0316227766016838 , 2.691083293730485964e-25 },
568         { 10 , 0.1 , 2.6905328954342172306e-20 },
569         { 10 , 0.316227766016838 , 2.6850345850670040022e-15 },
570         { 10 , 1 , 2.630615123687452921e-10 },
571         { 10 , 2 , 2.5153862827167368199e-07 },
572         { 10 , 3 , 1.292835164571588302e-05 },
573         { 10 , 4 , 0.00019504055466003448463 },
574         { 10 , 5 , 0.0014678026473104743583 },
575         { 10 , 6 , 0.0069639810027903158857 },
576         { 10 , 7 , 0.023539344388267140901 },
577         { 10 , 7.5 , 0.038998257889412211996 },
578         { 10 , 7.6 , 0.042818673234280582585 },
579         { 10 , 7.7 , 0.04690017276527555512 },
580         { 10 , 7.8 , 0.051248883025765065713 },
581         { 10 , 7.9 , 0.055869872504109699407 },
582         { 10 , 8 , 0.060767026774251164944 },
583         { 10 , 8.1 , 0.065942923604934144954 },
584         { 10 , 8.2 , 0.071398709153595626975 },
585         { 10 , 8.3 , 0.077133976423868738648 },
586         { 10 , 8.4 , 0.083146647220432454151 },
587         { 10 , 8.5 , 0.089432858880587384753 },
588         { 10 , 9 , 0.12469409282831672714 },
589         { 10 , 10 , 0.20748610663335886883 },
590         { 10 , 11 , 0.28042823052537591 },
591         { 10 , 12 , 0.300476035271269315 },
592         { 10 , 13 , 0.23378201020301889179 },
593         { 10 , 14 , 0.085006705446061009424 },
594         { 10 , 15 , -0.09007181104765905888 },
595         { 10 , 16 , -0.20620569442259728543 },
596         { 10 , 17 , -0.19911331972770593413 },
597         { 10 , 18 , -0.073169659187521246535 },
598         { 10 , 19 , 0.091553331622639774756 },
599         { 10 , 20 , 0.18648255802394508862 },
600         { 10 , 30 , -0.129876893998588816 },
601         { 10 , 100 , -0.054732176935472012791 },
602         { 10 , 300 , 0.027563483890691235778 },
603         { 10 , 1000 , -0.02452062230603655954 },
604         { 30 , 1e-08 , 3.511074584737334481e-282 },
605         { 30 , 3.16227766016838e-08 , 3.5110745847373276748e-267 },
606         { 30 , 1e-07 , 3.5110745847373271436e-252 },
607         { 30 , 3.16227766016838e-07 , 3.5110745847373244839e-237 },
608         { 30 , 1e-06 , 3.5110745847372989351e-222 },
609         { 30 , 3.16227766016838e-06 , 3.511074584737044636e-207 },
610         { 30 , 1e-05 , 3.5110745847345094386e-192 },
611         { 30 , 3.16227766016838e-05 , 3.5110745847090235522e-177 },
612         { 30 , 1e-04 , 3.5110745844541855471e-162 },
613         { 30 , 0.000316227766016838 , 3.5110745819058229075e-147 },
614         { 30 , 0.001 , 3.5110745564222159037e-132 },
615         { 30 , 0.00316227766016838 , 3.5110743015861690319e-117 },
616         { 30 , 0.01 , 3.5110717532266786188e-102 },
617         { 30 , 0.0316227766016838 , 3.5110462697303107185e-87 },
618         { 30 , 0.1 , 3.5107914446214635799e-72 },
619         { 30 , 0.316227766016838 , 3.5082441787554764315e-57 },
620         { 30 , 1 , 3.4828697942514824077e-42 },
621         { 30 , 2 , 3.6502562664740960186e-33 },
622         { 30 , 3 , 6.7223399381463293316e-28 },
623         { 30 , 4 , 3.5570357020361055268e-24 },
624         { 30 , 5 , 2.6711772782507989195e-21 },
625         { 30 , 6 , 5.7984683652785706951e-19 },
626         { 30 , 7 , 5.3172607940100176027e-17 },
627         { 30 , 7.5 , 3.9705139492720914996e-16 },
628         { 30 , 7.6 , 5.8351206236969734897e-16 },
629         { 30 , 7.7 , 8.5295046954365007979e-16 },
630         { 30 , 7.8 , 1.240300099862031423e-15 },
631         { 30 , 7.9 , 1.7943809060373146352e-15 },
632         { 30 , 8 , 2.5830997825663086363e-15 },
633         { 30 , 8.1 , 3.7004810818946501642e-15 },
634         { 30 , 8.2 , 5.2761304350589830578e-15 },
635         { 30 , 8.3 , 7.4879207291538333461e-15 },
636         { 30 , 8.4 , 1.057892772982890842e-14 },
637         { 30 , 8.5 , 1.4879948521285087748e-14 },
638         { 30 , 9 , 7.6921564693354977569e-14 },
639         { 30 , 10 , 1.5510960782574666161e-12 },
640         { 30 , 11 , 2.2735383676316185421e-11 },
641         { 30 , 12 , 2.5522590430344176732e-10 },
642         { 30 , 13 , 2.2828783239868354402e-09 },
643         { 30 , 14 , 1.6775399533577877891e-08 },
644         { 30 , 15 , 1.0374710201078721135e-07 },
645         { 30 , 16 , 5.5052386643076382366e-07 },
646         { 30 , 17 , 2.5460065118711982301e-06 },
647         { 30 , 18 , 1.0393652487465728599e-05 },
648         { 30 , 19 , 3.7849142225173515583e-05 },
649         { 30 , 20 , 0.00012401536360354329497 },
650         { 30 , 30 , 0.14393585001030734238 },
651         { 30 , 100 , 0.081460129581172213697 },
652         { 30 , 300 , -0.029514887800373371812 },
653         { 30 , 1000 , -0.020271896981075843147 },
654         { 100 , 1e-08 , 0 },
655         { 100 , 3.16227766016838e-08 , 0 },
656         { 100 , 1e-07 , 0 },
657         { 100 , 3.16227766016838e-07 , 0 },
658         { 100 , 1e-06 , 0 },
659         { 100 , 3.16227766016838e-06 , 0 },
660         { 100 , 1e-05 , 0 },
661         { 100 , 3.16227766016838e-05 , 0 },
662         { 100 , 1e-04 , 0 },
663         { 100 , 0.000316227766016838 , 0 },
664         { 100 , 0.001 , 0 },
665         { 100 , 0.00316227766016838 , 0 },
666         { 100 , 0.01 , 0 },
667         { 100 , 0.0316227766016838 , 0 },
668         { 100 , 0.1 , 8.4525165351217888791e-289 },
669         { 100 , 0.316227766016838 , 8.4506337559752745816e-239 },
670         { 100 , 1 , 8.4318287896267070128e-189 },
671         { 100 , 2 , 1.0609531124391718917e-158 },
672         { 100 , 3 , 4.260360181132621405e-141 },
673         { 100 , 4 , 1.305547836452271925e-128 },
674         { 100 , 5 , 6.2677893955418752099e-119 },
675         { 100 , 6 , 5.0513258541507019365e-111 },
676         { 100 , 7 , 2.4215591572118171706e-104 },
677         { 100 , 7.5 , 2.3583800455568589368e-101 },
678         { 100 , 7.6 , 8.8352979458474109476e-101 },
679         { 100 , 7.7 , 3.253025120751429903e-100 },
680         { 100 , 7.8 , 1.1776236102157393805e-99 },
681         { 100 , 7.9 , 4.1933885427120016432e-99 },
682         { 100 , 8 , 1.4694094093552327336e-98 },
683         { 100 , 8.1 , 5.0688862671208964077e-98 },
684         { 100 , 8.2 , 1.7220304874625643909e-97 },
685         { 100 , 8.3 , 5.7635248300942440709e-97 },
686         { 100 , 8.4 , 1.9011188242236321325e-96 },
687         { 100 , 8.5 , 6.182346491260611201e-96 },
688         { 100 , 9 , 1.8369106342703587456e-93 },
689         { 100 , 10 , 6.5973160641553802341e-89 },
690         { 100 , 11 , 8.6297901331738815878e-85 },
691         { 100 , 12 , 4.8983704457507876536e-81 },
692         { 100 , 13 , 1.3781127544328333402e-77 },
693         { 100 , 14 , 2.1310751903146119988e-74 },
694         { 100 , 15 , 1.9660095611249536378e-71 },
695         { 100 , 16 , 1.1559435724349575529e-68 },
696         { 100 , 17 , 4.5721265690179434188e-66 },
697         { 100 , 18 , 1.2722370655682102766e-63 },
698         { 100 , 19 , 2.5856336302772506687e-61 },
699         { 100 , 20 , 3.9617550943362506795e-59 },
700         { 100 , 30 , 4.5788015281752424119e-42 },
701         { 100 , 100 , 0.09636667329586150188 },
702         { 100 , 300 , -0.014491227064785699996 },
703         { 100 , 1000 , 0.011676135007802557891 },
704         { 300 , 1e-08 , 0 },
705         { 300 , 3.16227766016838e-08 , 0 },
706         { 300 , 1e-07 , 0 },
707         { 300 , 3.16227766016838e-07 , 0 },
708         { 300 , 1e-06 , 0 },
709         { 300 , 3.16227766016838e-06 , 0 },
710         { 300 , 1e-05 , 0 },
711         { 300 , 3.16227766016838e-05 , 0 },
712         { 300 , 1e-04 , 0 },
713         { 300 , 0.000316227766016838 , 0 },
714         { 300 , 0.001 , 0 },
715         { 300 , 0.00316227766016838 , 0 },
716         { 300 , 0.01 , 0 },
717         { 300 , 0.0316227766016838 , 0 },
718         { 300 , 0.1 , 0 },
719         { 300 , 0.316227766016838 , 0 },
720         { 300 , 1 , 0 },
721         { 300 , 2 , 0 },
722         { 300 , 3 , 0 },
723         { 300 , 4 , 0 },
724         { 300 , 5 , 0 },
725         { 300 , 6 , 0 },
726         { 300 , 7 , 0 },
727         { 300 , 7.5 , 0 },
728         { 300 , 7.6 , 0 },
729         { 300 , 7.7 , 0 },
730         { 300 , 7.8 , 0 },
731         { 300 , 7.9 , 0 },
732         { 300 , 8 , 0 },
733         { 300 , 8.1 , 0 },
734         { 300 , 8.2 , 0 },
735         { 300 , 8.3 , 0 },
736         { 300 , 8.4 , 0 },
737         { 300 , 8.5 , 0 },
738         { 300 , 9 , 0 },
739         { 300 , 10 , 0 },
740         { 300 , 11 , 0 },
741         { 300 , 12 , 0 },
742         { 300 , 13 , 0 },
743         { 300 , 14 , 0 },
744         { 300 , 15 , 0 },
745         { 300 , 16 , 0 },
746         { 300 , 17 , 0 },
747         { 300 , 18 , 0 },
748         { 300 , 19 , 0 },
749         { 300 , 20 , 0 },
750         { 300 , 30 , 1.0388021531643495593e-262 },
751         { 300 , 100 , 3.5203666218469330448e-109 },
752         { 300 , 300 , 0.066818398128979980544 },
753         { 300 , 1000 , 0.00046782803879124944908 }
754     };
755 
756     @Test
757     public void testBesselJ() {
758         final double tol = 1e-15;
759 
760         for (int i = 0; i < BESSEL_J_REF.length; i++) {
761             final double[] data = BESSEL_J_REF[i];
762             final double order = data[0];
763             final double x = data[1];
764             final double expected = data[2];
765             final double actual = BesselJ.value(order, x);
766 
767             String msg = "" + order + " @ " + x;
768             Assert.assertEquals(msg, expected, actual, tol);
769         }
770     }
771 
772     @Test(expected=MathIllegalArgumentException.class)
773     public void testIAEBadOrder() {
774         BesselJ.value(-1, 1);
775     }
776 
777     @Test(expected=MathIllegalArgumentException.class)
778     public void testIAEBadArgument() {
779         BesselJ.value(1, 100000);
780     }
781 }