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

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

Initial import--MML version from SOLEIL@2013
File size: 1011 bytes
Line 
1function pass = besseljroots
2
3% TAD
4
5g = chebfun(@(x) besselj(0,x),[0 100]);
6r = roots(g);
7exact = [ 2.40482555769577276862163; 5.52007811028631064959660
88.65372791291101221695437; 11.7915344390142816137431
914.9309177084877859477626; 18.0710639679109225431479
1021.2116366298792589590784; 24.3524715307493027370579
1127.4934791320402547958773; 30.6346064684319751175496
1233.7758202135735686842385; 36.9170983536640439797695
1340.0584257646282392947993; 43.1997917131767303575241
1446.3411883716618140186858; 49.4826098973978171736028
1552.6240518411149960292513; 55.7655107550199793116835
1658.9069839260809421328344; 62.0484691902271698828525
1765.1899648002068604406360; 68.3314693298567982709923
1871.4729816035937328250631; 74.6145006437018378838205
1977.7560256303880550377394; 80.8975558711376278637723
2084.0390907769381901578795; 87.1806298436411536512617
2190.3221726372104800557177; 93.4637187819447741711905
2296.6052679509962687781216; 99.7468198586805964702799
23];
24
25pass = norm(r-exact,Inf) < 4e-13*chebfunpref('eps')/eps;
26
Note: See TracBrowser for help on using the repository browser.