source: MML/trunk/machine/SOLEIL/common/toolbox/chebfun_v2_0501/chebfun/chebtests/besseljextrema.m @ 4

Last change on this file since 4 was 4, checked in by zhangj, 10 years ago

Initial import--MML version from SOLEIL@2013
File size: 837 bytes
Line 
1function pass = besseljextrema
2
3% TAD
4
5J2 = chebfun( @(x) besselj(2,x), [0 100] );
6extrema = roots(diff(J2));
7
8mma = [ 0
9   3.054236928227140
10   6.706133194158459
11   9.969467823087596
12  13.170370856016124
13  16.347522318321783
14  19.512912782488204
15  22.671581772477424
16  25.826037141785264
17  28.977672772993678
18  32.127327020443474
19  35.275535050674691
20  38.422654817555909
21  41.568934936074314
22  44.714553532819735
23  47.859641607992096
24  51.004297672458868
25  54.148597242671237
26  57.292599186428227
27  60.436350075253564
28  63.579887238154626
29  66.723240947717301
30  69.866436013337719
31  73.009492961171489
32  76.152428920759021
33  79.295258300056489
34  82.437993305560212
35  85.580644347487763
36  88.723220358610419
37  91.865729047477643
38  95.008177101267663
39  98.150570349583958];
40
41pass = norm( extrema-mma, Inf) < 1e-12*chebfunpref('eps')/eps;
Note: See TracBrowser for help on using the repository browser.