Add references
Pierre Marechal [Tue, 15 Apr 2008 06:44:55 +0000 (06:44 +0000)]
13 files changed:
scilab/modules/polynomials/tests/nonreg_tests/bug_2834.dia.ref
scilab/modules/polynomials/tests/nonreg_tests/bug_2834.tst
scilab/modules/sound/tests/nonreg_tests/bug_141.dia.ref [new file with mode: 0644]
scilab/modules/sparse/tests/nonreg_tests/bug_119.dia.ref [new file with mode: 0644]
scilab/modules/sparse/tests/nonreg_tests/bug_119.tst
scilab/modules/sparse/tests/nonreg_tests/bug_2277.dia.ref [new file with mode: 0644]
scilab/modules/special_functions/tests/nonreg_tests/bug_132.dia.ref [new file with mode: 0644]
scilab/modules/special_functions/tests/nonreg_tests/bug_132.tst
scilab/modules/special_functions/tests/nonreg_tests/bug_421.dia.ref [new file with mode: 0644]
scilab/modules/special_functions/tests/nonreg_tests/bug_421.tst
scilab/modules/special_functions/tests/unit_tests/besselj.dia.ref [new file with mode: 0644]
scilab/modules/special_functions/tests/unit_tests/erf.dia.ref [new file with mode: 0644]
scilab/modules/special_functions/tests/unit_tests/erfinv.dia.ref [new file with mode: 0644]

index 526237e..f0062e8 100644 (file)
 // the code below produced n=2 instead of 1
 n1=poly([4/3 1-%eps/2],'x','c');
 d1=poly([4/3 1],'x','c');
-[n,d]=simp(n1,d1)
- d  =
-    1   
- n  =
-    2   
+[n,d]=simp(n1,d1);
 if or(abs(coeff(n)-1)>10*%eps) then bugmes();quit;end
index e4edba5..9cd951f 100644 (file)
@@ -16,5 +16,5 @@
 
 n1=poly([4/3 1-%eps/2],'x','c');
 d1=poly([4/3 1],'x','c');
-[n,d]=simp(n1,d1) 
+[n,d]=simp(n1,d1);
 if or(abs(coeff(n)-1)>10*%eps) then pause,end
diff --git a/scilab/modules/sound/tests/nonreg_tests/bug_141.dia.ref b/scilab/modules/sound/tests/nonreg_tests/bug_141.dia.ref
new file mode 100644 (file)
index 0000000..bb76184
--- /dev/null
@@ -0,0 +1,17 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA - Vincent COUVERT
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 141 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=141
+//
+// <-- Short Description -->
+//    savewave crashes my machine
+if execstr("savewave(TMPDIR+""/toto.WAV"", rand(1:10000), 11025);", "errcatch")<>0 then pause; end
+Writing Wave file: Microsoft PCM format, 1 channel, 11025 samp/sec
+        22050 byte/sec, 2 block align, 16 bits/samp
+Finished writing Wave file, 20000 data bytes
diff --git a/scilab/modules/sparse/tests/nonreg_tests/bug_119.dia.ref b/scilab/modules/sparse/tests/nonreg_tests/bug_119.dia.ref
new file mode 100644 (file)
index 0000000..4e90716
--- /dev/null
@@ -0,0 +1,42 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA - Vincent COUVERT
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 119 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=119
+//
+// <-- Short Description -->
+//    when running sample from readmps help it crashes complete SCILAB when readmps is invoked.
+//Generate MPS file
+txt=['NAME          TESTPROB'
+     'ROWS'
+     ' N  COST'
+     ' L  LIM1'
+     ' G  LIM2'
+     ' E  MYEQN'
+     'COLUMNS'
+     '    XONE      COST                 1   LIM1                 1'
+     '    XONE      LIM2                 1'
+     '    YTWO      COST                 4   LIM1                 1'
+     '    YTWO      MYEQN               -1'
+     '    ZTHREE    COST                 9   LIM2                 1'
+     '    ZTHREE    MYEQN                1'
+     'RHS'
+     '    RHS1      LIM1                 5   LIM2                10'
+     '    RHS1      MYEQN                7'
+     'BOUNDS'
+     ' UP BND1      XONE                 4'
+     ' LO BND1      YTWO                -1'
+     ' UP BND1      YTWO                 1'
+     'ENDATA'];
+mputl(txt,TMPDIR+'/test.mps');
+//Read the MPS file
+P=readmps(TMPDIR+'/test.mps',[0 10^30]);
+//Convert it to linpro format
+LP=mps2linpro(P);
+//Solve it with linpro
+[x,lagr,f]=linpro(LP(2:$));
index 346ee39..82b7a9c 100644 (file)
@@ -13,7 +13,7 @@
 // <-- Short Description -->
 //    when running sample from readmps help it crashes complete SCILAB when readmps is invoked.
 
-//Generate MPS file 
+//Generate MPS file
 txt=['NAME          TESTPROB'
      'ROWS'
      ' N  COST'
@@ -35,10 +35,10 @@ txt=['NAME          TESTPROB'
      ' LO BND1      YTWO                -1'
      ' UP BND1      YTWO                 1'
      'ENDATA'];
-mputl(txt,TMPDIR+'/test.mps')
+mputl(txt,TMPDIR+'/test.mps');
 //Read the MPS file
-P=readmps(TMPDIR+'/test.mps',[0 10^30])
+P=readmps(TMPDIR+'/test.mps',[0 10^30]);
 //Convert it to linpro format
-LP=mps2linpro(P)
+LP=mps2linpro(P);
 //Solve it with linpro
-[x,lagr,f]=linpro(LP(2:$))
+[x,lagr,f]=linpro(LP(2:$));
diff --git a/scilab/modules/sparse/tests/nonreg_tests/bug_2277.dia.ref b/scilab/modules/sparse/tests/nonreg_tests/bug_2277.dia.ref
new file mode 100644 (file)
index 0000000..a3bc130
--- /dev/null
@@ -0,0 +1,26 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 2277 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=2277
+//
+// <-- Short Description -->
+//    Scilab has a bug with the SPARSE command. When the "values vector" in the
+//    row/column representation is COMPLEX the routine fails to build the sparse
+//    matrix correctly. This seems to occur always when the first AND third entry of
+//    this vector are zero.
+//    The other bug is that when other elements are zero (not the third, the fifth for
+//    example) they are still stored in the sparse representation.
+//    This is very problematic for people doing finite element calculations with scilab!
+ij          = [1,3;1,7;3,3;4,4;5,3;6,10;7,1;7,9;8,10;9,6;9,7;9,9;10,10] ;
+v           = [1 0 0 4:13]'*%i;
+A           = sparse(ij,v,[10 10]);
+[ij2,v2,mn] = spget(A);
+if or(v2<>[1 4:13]'*%i) then bugmes();quit;end
+ij(2:3,:)   = [];
+if or(ij<>ij2) then bugmes();quit;end
diff --git a/scilab/modules/special_functions/tests/nonreg_tests/bug_132.dia.ref b/scilab/modules/special_functions/tests/nonreg_tests/bug_132.dia.ref
new file mode 100644 (file)
index 0000000..4ae9575
--- /dev/null
@@ -0,0 +1,14 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA - Vincent COUVERT
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 132 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=132
+//
+// <-- Short Description -->
+//    bessely(0,1) produces convergence error
+if execstr("bessely(0,1);", "errcatch")<>0 then bugmes();quit;end
index f15161a..e0108b5 100644 (file)
@@ -14,5 +14,5 @@
 //    bessely(0,1) produces convergence error
 
 
-if execstr("bessely(0,1);", "errcatch")<>0 then pause; end
+if execstr("bessely(0,1);", "errcatch")<>0 then pause,end
 
diff --git a/scilab/modules/special_functions/tests/nonreg_tests/bug_421.dia.ref b/scilab/modules/special_functions/tests/nonreg_tests/bug_421.dia.ref
new file mode 100644 (file)
index 0000000..d4088fc
--- /dev/null
@@ -0,0 +1,15 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) 2008 - INRIA - Vincent COUVERT
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+// <-- Non-regression test for bug 421 -->
+//
+// <-- Bugzilla URL -->
+// http://bugzilla.scilab.org/show_bug.cgi?id=421
+//
+// <-- Short Description -->
+//    gamma gives wrong results
+if gamma(200)<>%inf then bugmes();quit;end
+if gamma(1000)<>%inf then bugmes();quit;end
index bdeb47d..23549f4 100644 (file)
@@ -14,7 +14,7 @@
 //    gamma gives wrong results
 
 
-if gamma(200)<>%inf then pause; end
+if gamma(200)<>%inf then pause,end
 
-if gamma(1000)<>%inf then pause; end
+if gamma(1000)<>%inf then pause,end
 
diff --git a/scilab/modules/special_functions/tests/unit_tests/besselj.dia.ref b/scilab/modules/special_functions/tests/unit_tests/besselj.dia.ref
new file mode 100644 (file)
index 0000000..2084184
--- /dev/null
@@ -0,0 +1,1211 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) ????-2008 - INRIA
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+function r=testresult(eps)
+  if argn(2)<=0 then eps=1d-10;end
+  if or(isnan(real(y))<>isnan(real(yref)))  then r=%t,return,end;
+  if or(isnan(imag(y))<>isnan(imag(yref)))  then r=%t,return,end;
+  if or(size(y)<>size(yref)) then r=%t,return,end;
+  if yref==[] then r=%f,return,end;
+  k=find(abs(yref)==0);
+  r=or(abs(y(k)-yref(k))>eps);
+  k=find(abs(yref)<>0);
+  r=r|or(max(abs(y(k)-yref(k))./abs(yref(k)))>eps);
+endfunction
+ieee(1)
+// ------------------------------ besselj(alpha,x)------------------------------
+//test of numerical precision reference computed by maple
+besselJ=besselj(alpha',x);
+if norm((besselJ-besselJ_ref)./abs(besselJ_ref),1)>1d6*%eps then bugmes();quit;end
+alpha=[];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[0];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[9.99999999999999955e-08];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[1];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[3 4 5];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[3
+4
+5];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[0];
+yref=[1];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[1];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[9.99999999999999955e-08];
+yref=[0.999999999999997446];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.999999999999997446];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[1];
+yref=[0.765197686557966605];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.765197686557966605];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[3 4 5];
+yref=[-0.260051954901933557 -0.397149809863847236 -0.177596771314338348];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.260051954901933557 -0.397149809863847236 -0.177596771314338348];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[0.288216947635014376-1.63247121237385119e-16*%i 0.0691572616569850473+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.288216947635014376-1.63247121237385119e-16*%i 0.0691572616569850473+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[3
+4
+5];
+yref=[-0.260051954901933557
+-0.397149809863847236
+-0.177596771314338348];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.260051954901933557
+-0.397149809863847236
+-0.177596771314338348];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[ 0.288216947635014376-1.63247121237385119e-16*%i
+0.0691572616569850473+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.288216947635014376-1.63247121237385119e-16*%i
+0.0691572616569850473+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[0.288272330491939344+0.00165156721805876858*%i 0.0691646458471141812-0.00541608716438180676*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.282564155899245439+0.0016188639960181446*%i 0.0677950940990729028-0.00530884145304078969*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[0];
+yref=[0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[9.99999999999999955e-08];
+yref=[4.99999999999999647e-08];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[4.99999999999999647e-08];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[1];
+yref=[0.440050585744933553];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.440050585744933553];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[3 4 5];
+yref=[0.339058958525936538 -0.0660433280235490638 -0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.339058958525936538 -0.0660433280235490638 -0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i 0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i 0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[3
+4
+5];
+yref=[ 0.339058958525936538
+-0.0660433280235490638
+ -0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.339058958525936538
+-0.0660433280235490638
+ -0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i
+ 0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i
+ 0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[-0.0825942221022646084+0.00553845279173677783*%i 0.270840538797723684+0.000738422002625521083*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0809587469274432592+0.0054287840786324832*%i 0.265477536807215531+0.000723800267314064468*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[0];
+yref=[0 0 0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0 0 0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[9.99999999999999955e-08];
+yref=[4.99999999999999647e-08 1.25e-15 2.08333333333333265e-23];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[4.99999999999999647e-08 1.25e-15 2.08333333333333265e-23];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[1];
+yref=[0.440050585744933553 0.114903484931900501 0.0195633539826684105];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.440050585744933553 0.114903484931900501 0.0195633539826684105];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[3 4 5];
+yref=[0.339058958525936538 0.36412814585207276 0.364831230613666901];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.339058958525936538 0.36412814585207276 0.364831230613666901];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[3
+4
+5];
+yref=[ 0.339058958525936649 0.486091260585891194 0.30906272225525161
+-0.0660433280235490083 0.364128145852072926 0.430171473875621935
+ -0.327579137591465064 0.046565116277752297 0.364831230613666901];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.339058958525936649 0.486091260585891194 0.30906272225525161
+-0.0660433280235490083 0.364128145852072926 0.430171473875621935
+ -0.327579137591465064 0.046565116277752297 0.364831230613666901];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1 2 3];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i -0.265594911883436868+1.66696379225679121e-16*%i 0.228101889059524748-1.43164485837168552e-16*%i
+ 0.270786268276835385+0*%i -0.00468434063869093212+0*%i -0.273016906676211979+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0825704304932580041+5.18240040696731648e-17*%i -0.265594911883436868+1.66696379225679121e-16*%i 0.228101889059524748-1.43164485837168552e-16*%i
+ 0.270786268276835385+0*%i -0.00468434063869093212+0*%i -0.273016906676211979+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[0];
+yref=[-0 0 0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0 0 0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[9.99999999999999955e-08];
+yref=[1.57275203335608714e-09 3.57443643944564709e-17 5.585056936633812e-25];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[1.57275203335608714e-09 3.57443643944564709e-17 5.585056936633812e-25];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[1];
+yref=[0.351883806529722176 0.0829767227596755474 0.0132137736128501541];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.351883806529722176 0.0829767227596755474 0.0132137736128501541];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[3 4 5];
+yref=[0.412044162415463822 0.406986020132752813 0.391959608420348327];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.412044162415463822 0.406986020132752813 0.391959608420348327];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[3
+4
+5];
+yref=[ 0.412044162415463822 0.463422297657139715 0.267641874148340286
+0.0388438196375475928 0.406986020132752813 0.40884080250847965
+ -0.2772936869328117 0.130302183508563874 0.391959608420348327];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.412044162415463822 0.463422297657139715 0.267641874148340286
+0.0388438196375475928 0.406986020132752813 0.40884080250847965
+ -0.2772936869328117 0.130302183508563874 0.391959608420348327];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[1.19999999999999996 2.20000000000000018 3.20000000000000018];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[0.000439741362268266983+0.000319490801010906766*%i -0.237705022442928571-0.172702807925255664*%i 0.142834518740318772+0.103775352332019871*%i
+ 0.244808700524772438+0*%i -0.0775662890440063368+0*%i -0.285438661452585507+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.000439741362268266983+0.000319490801010906766*%i -0.237705022442928571-0.172702807925255664*%i 0.142834518740318772+0.103775352332019871*%i
+ 0.244808700524772438+0*%i -0.0775662890440063368+0*%i -0.285438661452585507+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[0];
+yref=[1
+0
+0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[1
+0
+0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[9.99999999999999955e-08];
+yref=[ 0.999999999999997446
+4.99999999999999647e-08
+ 1.25e-15];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.999999999999997446
+4.99999999999999647e-08
+ 1.25e-15];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[1];
+yref=[0.765197686557966605
+0.440050585744933553
+0.114903484931900501];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.765197686557966605
+0.440050585744933553
+0.114903484931900501];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[3 4 5];
+yref=[-0.260051954901933557 0.339058958525936482 0.486091260585891194
+-0.397149809863847236 -0.0660433280235490638 0.36412814585207276
+-0.177596771314338348 -0.327579137591465064 0.046565116277752297];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.260051954901933557 0.339058958525936482 0.486091260585891194
+-0.397149809863847236 -0.0660433280235490638 0.36412814585207276
+-0.177596771314338348 -0.327579137591465064 0.046565116277752297];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[ 0.288216947635014376-1.63247121237385119e-16*%i -0.0825704304932580041+4.67681903786785266e-17*%i -0.265594911883436868+1.5043391839391156e-16*%i
+0.0691572616569850473+0*%i 0.270786268276835385+0*%i -0.00468434063869093212+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.288216947635014376-1.63247121237385119e-16*%i -0.0825704304932580041+4.67681903786785266e-17*%i -0.265594911883436868+1.5043391839391156e-16*%i
+0.0691572616569850473+0*%i 0.270786268276835385+0*%i -0.00468434063869093212+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[3
+4
+5];
+yref=[ -0.260051954901933557
+-0.0660433280235490638
+ 0.046565116277752297];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -0.260051954901933557
+-0.0660433280235490638
+ 0.046565116277752297];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0
+1
+2];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[ 0.288272330491939344+0.00165156721805876858*%i -0.0825942221022646084+0.00553845279173677263*%i -0.265639789156075279-0.0031069445819979849*%i
+0.0691646458471141812-0.00541608716438180676*%i 0.270840538797723684+0.000738422002625537671*%i -0.00467875023595177783+0.00543836408021845449*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.282564155899245439+0.0016188639960181446*%i -0.0809587469274432592+0.005428784078632478*%i -0.260379768908271136-0.00304542295731204116*%i
+0.0677950940990729028-0.00530884145304078969*%i 0.265477536807215531+0.000723800267314080731*%i -0.00458610477401360205+0.00533067725638924125*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[0];
+yref=[ 0
+-0
+ 0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0
+-0
+ 0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[9.99999999999999955e-08];
+yref=[ 0.722522753107868621
+3.54177820150916157e-08
+8.76677772650784535e-16];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.722522753107868621
+3.54177820150916157e-08
+8.76677772650784535e-16];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[1];
+yref=[0.767624213200681971
+0.430847058235354519
+0.111303785599441221];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.767624213200681971
+0.430847058235354519
+0.111303785599441221];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[3 4 5];
+yref=[-0.248113342531625397 0.347410321977723313 0.4843523614764772
+-0.397488017750937228 -0.0555197528009771721 0.369172943822438926
+-0.187188738237360264 -0.323527150348924653 0.0551896608949990139];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.248113342531625397 0.347410321977723313 0.4843523614764772
+-0.397488017750937228 -0.0555197528009771721 0.369172943822438926
+-0.187188738237360264 -0.323527150348924653 0.0551896608949990139];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[ 0.289473237714468368+0.0182121124296419773*%i -0.0742460772355384713-0.00467116724416256467*%i -0.268725018870564547-0.0169067451449718964*%i
+0.0774792372615417391+0*%i 0.269158266698860282+0*%i -0.0121122296346756614+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.289473237714468368+0.0182121124296419773*%i -0.0742460772355384713-0.00467116724416256467*%i -0.268725018870564547-0.0169067451449718964*%i
+0.0774792372615417391+0*%i 0.269158266698860282+0*%i -0.0121122296346756614+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[3
+4
+5];
+yref=[ -0.24811334253162734
+-0.0555197528009771929
+ 0.0551896608949990139];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -0.24811334253162734
+-0.0555197528009771929
+ 0.0551896608949990139];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[0.0200000000000000004
+ 1.02000000000000002
+ 2.02000000000000002];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[ 0.289436686533186605+0.0196848354538640968*%i -0.0746194938467206886+0.000909756330523188092*%i -0.2685835756459789-0.0198819369191412312*%i
+0.0774883289690142946-0.00537983661590852799*%i 0.269212372909699293+0.0008959309889286695*%i -0.0121080338436329836+0.00544175248673077396*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.283705456146132684+0.0192950495961393617*%i -0.0731419288713772014+0.000891741948211250118*%i -0.263265264520173126-0.0194882481909108424*%i
+0.0759539572521852335-0.00527330871352064186*%i 0.263881610763850749+0.000878190366722298409*%i -0.0118682787098823468+0.00533399856795724152*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[0];
+yref=[-0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[9.99999999999999955e-08];
+yref=[-4.99999999999999647e-08];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-4.99999999999999647e-08];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[1];
+yref=[-0.440050585744933553];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.440050585744933553];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[3 4 5];
+yref=[-0.339058958525936538 0.0660433280235490638 0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.339058958525936538 0.0660433280235490638 0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i -0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i -0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[3
+4
+5];
+yref=[-0.339058958525936538
+0.0660433280235490638
+ 0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.339058958525936538
+0.0660433280235490638
+ 0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i
+-0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i
+-0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-1];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[0.0825942221022646084-0.00553845279173677783*%i -0.270840538797723684-0.000738422002625521083*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0809587469274432592-0.0054287840786324832*%i -0.265477536807215531-0.000723800267314064468*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[0];
+yref=[ -0
+-%inf];
+y=besselj(alpha,x);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+yref=[ -0
+-%inf];
+y=besselj(alpha,x,1);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[9.99999999999999955e-08];
+yref=[-4.99999999999999647e-08
+ -25231325220.2017288];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-4.99999999999999647e-08
+ -25231325220.2017288];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[1];
+yref=[-0.440050585744933553
+ -1.10249557516017949];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.440050585744933553
+ -1.10249557516017949];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[3 4 5];
+yref=[-0.339058958525936538 0.0870080907208351517
+0.0660433280235490638 0.367112032460934101
+ 0.327579137591465064 0.3219244429611402];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.339058958525936538 0.0870080907208351378
+0.0660433280235490638 0.367112032460934101
+ 0.327579137591465064 0.321924442961140256];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i 9.50628464835290288e-16-0.272424375376841277*%i
+-0.270786268276835385+0*%i -0.218249100772301557+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i 9.43689570931383059e-16-0.272424375376841332*%i
+-0.270786268276835385+0*%i -0.218249100772301585+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[-7.29999999999999982
+ 8.40000000000000036];
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i
+-0.218249100772301557+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0825704304932580041-5.18240040696731648e-17*%i
+-0.218249100772301585+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1
+-1.5];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[0.0825942221022646084-0.00553845279173677783*%i -0.00198764870679882157-0.272479284229019014*%i
+-0.270840538797723684-0.000738422002625521083*%i -0.218287028059932248+0.0036388855639192105*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.0809587469274432592-0.0054287840786324832*%i -0.00194829062540415629-0.267083832904858764*%i
+-0.265477536807215531-0.000723800267314064468*%i -0.213964655304420243+0.00356683080206891736*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[0];
+yref=[-0
+ 0
+-0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0
+ 0
+-0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[9.99999999999999955e-08];
+yref=[-2.08333333333333265e-23
+ 1.25e-15
+-4.99999999999999647e-08];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-2.08333333333333265e-23
+ 1.25e-15
+-4.99999999999999647e-08];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[1];
+yref=[-0.0195633539826684105
+ 0.114903484931900501
+ -0.440050585744933553];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0195633539826684105
+ 0.114903484931900501
+ -0.440050585744933553];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[3 4 5];
+yref=[ -0.30906272225525161 0.486091260585891194 -0.339058958525936538
+-0.430171473875621935 0.36412814585207276 0.0660433280235490638
+-0.364831230613666901 0.046565116277752297 0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -0.30906272225525161 0.486091260585891194 -0.339058958525936538
+-0.430171473875621935 0.36412814585207276 0.0660433280235490638
+-0.364831230613666901 0.046565116277752297 0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[-0.228101889059524748+1.43164485837168552e-16*%i -0.265594911883436868+1.5043391839391156e-16*%i 0.0825704304932580041-5.18240040696731648e-17*%i
+ 0.273016906676211979+0*%i -0.00468434063869093212+0*%i -0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.228101889059524748+1.43164485837168552e-16*%i -0.265594911883436868+1.5043391839391156e-16*%i 0.0825704304932580041-5.18240040696731648e-17*%i
+ 0.273016906676211979+0*%i -0.00468434063869093212+0*%i -0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[3
+4
+5];
+yref=[-0.30906272225525161
+ 0.36412814585207276
+0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.30906272225525161
+ 0.36412814585207276
+0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-3
+-2
+-1];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[-0.228144514233508455+0.00343724948069338773*%i -0.265639789156075279-0.0031069445819979849*%i 0.0825942221022646084-0.00553845279173677783*%i
+ 0.273062336559663221-0.00185656517310216836*%i -0.00467875023595177783+0.00543836408021845449*%i -0.270840538797723684-0.000738422002625521083*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.223626950173899136+0.00336918738079998795*%i -0.260379768908271136-0.00304542295731204116*%i 0.0809587469274432592-0.0054287840786324832*%i
+ 0.267655340025824695-0.00181980271958226327*%i -0.00458610477401360205+0.00533067725638924125*%i -0.265477536807215531-0.000723800267314064468*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[0];
+yref=[ 0
+-0
+ 1
+ 0];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0
+-0
+ 1
+ 0];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[9.99999999999999955e-08];
+yref=[ 1.25e-15
+-4.99999999999999647e-08
+ 0.999999999999997446
+ 4.99999999999999647e-08];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 1.25e-15
+-4.99999999999999647e-08
+ 0.999999999999997446
+ 4.99999999999999647e-08];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[1];
+yref=[ 0.114903484931900501
+-0.440050585744933553
+ 0.765197686557966605
+ 0.440050585744933553];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.114903484931900501
+-0.440050585744933553
+ 0.765197686557966605
+ 0.440050585744933553];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[3 4 5];
+yref=[0.486091260585891194 -0.339058958525936538 -0.260051954901933557 0.339058958525936538
+ 0.36412814585207276 0.0660433280235490638 -0.397149809863847236 -0.0660433280235490638
+0.046565116277752297 0.327579137591465064 -0.177596771314338348 -0.327579137591465064];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[0.486091260585891194 -0.339058958525936538 -0.260051954901933557 0.339058958525936538
+ 0.36412814585207276 0.0660433280235490638 -0.397149809863847236 -0.0660433280235490638
+0.046565116277752297 0.327579137591465064 -0.177596771314338348 -0.327579137591465064];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[ -0.265594911883436868+1.5043391839391156e-16*%i 0.0825704304932580041-5.18240040696731648e-17*%i 0.288216947635014376-1.63247121237385119e-16*%i -0.0825704304932580041+5.18240040696731648e-17*%i
+-0.00468434063869093212+0*%i -0.270786268276835385+0*%i 0.0691572616569850473+0*%i 0.270786268276835385+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -0.265594911883436868+1.5043391839391156e-16*%i 0.0825704304932580041-5.18240040696731648e-17*%i 0.288216947635014376-1.63247121237385119e-16*%i -0.0825704304932580041+5.18240040696731648e-17*%i
+-0.00468434063869093212+0*%i -0.270786268276835385+0*%i 0.0691572616569850473+0*%i 0.270786268276835385+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[-2
+-1
+ 0
+ 1];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[ -0.265639789156075279-0.0031069445819979849*%i 0.0825942221022646084-0.00553845279173677783*%i 0.288272330491939344+0.00165156721805876858*%i -0.0825942221022646084+0.00553845279173677783*%i
+-0.00467875023595177783+0.00543836408021845449*%i -0.270840538797723684-0.000738422002625521083*%i 0.0691646458471141812-0.00541608716438180676*%i 0.270840538797723684+0.000738422002625521083*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -0.260379768908271136-0.00304542295731204116*%i 0.0809587469274432592-0.0054287840786324832*%i 0.282564155899245439+0.0016188639960181446*%i -0.0809587469274432592+0.0054287840786324832*%i
+-0.00458610477401360205+0.00533067725638924125*%i -0.265477536807215531-0.000723800267314064468*%i 0.0677950940990729028-0.00530884145304078969*%i 0.265477536807215531+0.000723800267314064468*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[0];
+yref=[-%inf
+-%inf
+ 0
+ -0];
+y=besselj(alpha,x);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+yref=[-%inf
+-%inf
+ 0
+ -0];
+y=besselj(alpha,x,1);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[9.99999999999999955e-08];
+yref=[ -716418914.617646098
+ 119.403152436272933
+8.52879660259094751e-06
+2.50846958899734603e-13];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ -716418914.617646098
+ 119.403152436272933
+8.52879660259094751e-06
+2.50846958899734603e-13];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[1];
+yref=[-0.962671286840223628
+ 0.633870726369385018
+ 0.582348851018592018
+ 0.181417665056644584];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.962671286840223739
+ 0.633870726369385018
+ 0.582348851018592018
+ 0.181417665056644584];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[3 4 5];
+yref=[-0.106536302179580183 -0.40675205644906709 0.187886713469395089 0.494432522734784052
+ 0.268869311811566036 -0.344459534063007056 -0.217200381702116863 0.268439400467269551
+ 0.359445845778418371 -0.0150494093195678302 -0.357639916660074997 -0.0850897673452510606];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.106536302179580156 -0.406752056449067145 0.187886713469395089 0.494432522734784052
+ 0.268869311811565981 -0.344459534063007056 -0.217200381702116863 0.268439400467269551
+ 0.359445845778418316 -0.0150494093195678025 -0.357639916660074997 -0.0850897673452510606];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[ 0.12542315847504959-0.172630167733447554*%i 0.133575862270551238-0.183851401840542522*%i -0.114444320480207362+0.157519093609566369*%i -0.111627636425035953+0.153642260600349023*%i
+-0.26289292206878756+0*%i -0.0606892109991253315+0*%i 0.267227865711584134+0*%i 0.105227188617724732+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.125423158475049479-0.172630167733447665*%i 0.13357586227055121-0.183851401840542578*%i -0.114444320480207362+0.157519093609566369*%i -0.111627636425035953+0.153642260600349023*%i
+-0.262892922068787505+0*%i -0.0606892109991253037+0*%i 0.267227865711584134+0*%i 0.105227188617724732+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[ 0.122388153838117283-0.174892716810286619*%i 0.136900984596519226-0.181484735590922708*%i -0.111095235613447238+0.160007102702057979*%i -0.115511137516278492+0.150857044242800997*%i
+-0.262941829098775337+0.00202766723967471967*%i -0.0607076451088290336-0.00530154123317297127*%i 0.267278966499405468-0.00165931224063856075*%i 0.10525322854217678+0.00491892828012802791*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.119964706020585635-0.171429608988457022*%i 0.134190163475896684-0.177891097051649821*%i -0.108895402559002416+0.156838749788214965*%i -0.113223863745610387+0.147869874625772002*%i
+-0.257735232039471018+0.00198751673823663755*%i -0.0595055531952515537-0.00519656368323724409*%i 0.261986488365517967-0.0016264556568755667*%i 0.103169074978294403+0.00482152697427259275*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[];
+yref=[];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[0];
+yref=[-%inf
+-%inf
+-%inf
+ 0
+ -0];
+y=besselj(alpha,x);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+yref=[-%inf
+-%inf
+-%inf
+ 0
+ -0];
+y=besselj(alpha,x,1);
+Warning : division by zero...
+
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[9.99999999999999955e-08];
+yref=[ 18626891780058596
+ -716418914.617646098
+ 119.403152436272933
+8.52879660259094751e-06
+2.50846958899734603e-13];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 18626891780058600
+ -716418914.617646098
+ 119.403152436272933
+8.52879660259094751e-06
+2.50846958899734603e-13];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[1];
+yref=[ 1.86907461941519504
+-0.962671286840223628
+ 0.633870726369385018
+ 0.582348851018592018
+ 0.181417665056644584];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 1.86907461941519526
+-0.962671286840223739
+ 0.633870726369385018
+ 0.582348851018592018
+ 0.181417665056644584];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[3 4 5];
+yref=[ 0.499083518338036658 -0.106536302179580183 -0.40675205644906709 0.187886713469395089 0.494432522734784052
+ 0.16969448138548901 0.268869311811566036 -0.344459534063007056 -0.217200381702116863 0.268439400467269551
+-0.171862430485209827 0.359445845778418371 -0.0150494093195678302 -0.357639916660074997 -0.0850897673452510606];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[ 0.499083518338036658 -0.106536302179580156 -0.406752056449067145 0.187886713469395089 0.494432522734784052
+ 0.169694481385489065 0.268869311811565981 -0.344459534063007056 -0.217200381702116863 0.268439400467269551
+-0.171862430485209772 0.359445845778418316 -0.0150494093195678025 -0.357639916660074997 -0.0850897673452510606];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[-7.29999999999999982 8.40000000000000036];
+yref=[-0.0889046003479301572+0.122366684565615697*%i 0.12542315847504959-0.172630167733447554*%i 0.133575862270551238-0.183851401840542522*%i -0.114444320480207362+0.157519093609566369*%i -0.111627636425035953+0.153642260600349023*%i
+ 0.142060829734702354+0*%i -0.26289292206878756+0*%i -0.0606892109991253315+0*%i 0.267227865711584134+0*%i 0.105227188617724732+0*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0889046003479302266+0.122366684565615863*%i 0.125423158475049479-0.172630167733447665*%i 0.13357586227055121-0.183851401840542578*%i -0.114444320480207362+0.157519093609566369*%i -0.111627636425035953+0.153642260600349023*%i
+ 0.142060829734702354+0*%i -0.262892922068787505+0*%i -0.0606892109991253037+0*%i 0.267227865711584134+0*%i 0.105227188617724732+0*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
+alpha=[ -2.29999999999999982
+ -1.30000000000000004
+-0.300000000000000044
+ 0.699999999999999956
+ 1.69999999999999996];
+x=[-7.29999999999999982+0.0200000000000000004*%i 8.40000000000000036+0.0200000000000000004*%i ];
+yref=[-0.0931403526417651506+0.119314071061107668*%i 0.122388153838117283-0.174892716810286619*%i 0.136900984596519226-0.181484735590922708*%i -0.111095235613447238+0.160007102702057979*%i -0.115511137516278492+0.150857044242800997*%i
+ 0.142092446056964028+0.00448015660911126913*%i -0.262941829098775337+0.00202766723967471967*%i -0.0607076451088290336-0.00530154123317297127*%i 0.267278966499405468-0.00165931224063856075*%i 0.10525322854217678+0.00491892828012802791*%i ];
+y=besselj(alpha,x);
+if testresult() then bugmes();quit;end
+yref=[-0.0912960500907815109+0.116951494160925568*%i 0.119964706020585635-0.171429608988457022*%i 0.134190163475896684-0.177891097051649821*%i -0.108895402559002416+0.156838749788214965*%i -0.113223863745610387+0.147869874625772002*%i
+ 0.139278827111947917+0.00439144356445757625*%i -0.257735232039471018+0.00198751673823663755*%i -0.0595055531952515537-0.00519656368323724409*%i 0.261986488365517967-0.0016264556568755667*%i 0.103169074978294403+0.00482152697427259275*%i ];
+y=besselj(alpha,x,1);
+if testresult() then bugmes();quit;end
diff --git a/scilab/modules/special_functions/tests/unit_tests/erf.dia.ref b/scilab/modules/special_functions/tests/unit_tests/erf.dia.ref
new file mode 100644 (file)
index 0000000..1cf935a
--- /dev/null
@@ -0,0 +1,1018 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) ????-2008 - INRIA Serge Steer
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+function [x,y]=erf_ref()
+//reference Yr=erf(X) generee par Maple
+//pm:=array(1..1001):x:=-5:
+//for k from 1  to 1001 do
+//pm[k]:=evalf[24](erf(x)):
+//x:=x+0.01:end do;
+  x=-5:0.01:5;
+  y=[-0.999999999998462540205572
+     -0.999999999998297733894900
+     -0.999999999998115631480261
+     -0.999999999997914458007001
+     -0.999999999997692260661541
+     -0.999999999997446891397149
+     -0.999999999997175987908301
+     -0.999999999996876952801260
+     -0.999999999996546930794861
+     -0.999999999996182783770717
+     -0.999999999995781063475994
+     -0.999999999995337981664472
+     -0.999999999994849377442675
+     -0.999999999994310681567322
+     -0.999999999993716877418075
+     -0.999999999993062458345374
+     -0.999999999992341381066975
+     -0.999999999991547014758366
+     -0.999999999990672085451468
+     -0.999999999989708615322670
+     -0.999999999988647856415078
+     -0.999999999987480218300757
+     -0.999999999986195189146334
+     -0.999999999984781249599477
+     -0.999999999983225778864110
+     -0.999999999981514952278515
+     -0.999999999979633629652374
+     -0.999999999977565233555991
+     -0.999999999975291616686970
+     -0.999999999972792917366275
+     -0.999999999970047402136203
+     -0.999999999967031294347156
+     -0.999999999963718587527495
+     -0.999999999960080842230859
+     -0.999999999956086964947412
+     -0.999999999951702967549089
+     -0.999999999946891705613252
+     -0.999999999941612593833703
+     -0.999999999935821296581860
+     -0.999999999929469391523360
+     -0.999999999922504004025582
+     -0.999999999914867409908580
+     -0.999999999906496603894840
+     -0.999999999897322830900966
+     -0.999999999887271077085832
+     -0.999999999876259517323697
+     -0.999999999864198915505960
+     -0.999999999850991973790366
+     -0.999999999836532626610029
+     -0.999999999820705274925103
+     -0.999999999803383955845711
+     -0.999999999784431442373954
+     -0.999999999763698267603729
+     -0.999999999741021667277623
+     -0.999999999716224434128131
+     -0.999999999689113676923764
+     -0.999999999659479476596644
+     -0.999999999627093431244556
+     -0.999999999591707081174197
+     -0.999999999553050204480830
+     -0.999999999510828972939411
+     -0.999999999464723957210342
+     -0.999999999414387969535759
+     -0.999999999359443731215952
+     -0.999999999299481351206210
+     -0.999999999234055601157792
+     -0.999999999162682971138542
+     -0.999999999084838489103913
+     -0.999999998999952285943071
+     -0.999999998907405886591731
+     -0.999999998806528206277959
+     -0.999999998696591229443244
+     -0.999999998576805347252390
+     -0.999999998446314327865434
+     -0.999999998304189891785759
+     -0.999999998149425862613257
+     -0.999999997980931861411821
+     -0.999999997797526510638154
+     -0.999999997597930111164948
+     -0.999999997380756753356393
+     -0.999999997144505820407811
+     -0.999999996887552839233372
+     -0.999999996608139631065173
+     -0.999999996304363710601733
+     -0.999999995974166879001692
+     -0.999999995615322952246180
+     -0.999999995225424562376994
+     -0.999999994801868964843021
+     -0.999999994341842780638679
+     -0.999999993842305597079599
+     -0.999999993299972345915102
+     -0.999999992711294372006456
+     -0.999999992072439099985541
+     -0.999999991379268200130508
+     -0.999999990627314148132540
+     -0.999999989811755066458970
+     -0.999999988927387727619661
+     -0.999999987968598591791738
+     -0.999999986929332742926881
+     -0.999999985803060578628974
+     -0.999999984582742099719981
+     -0.999999983260788635479157
+     -0.999999981829021830014495
+     -0.999999980278629704073472
+     -0.999999978600119594788978
+     -0.999999976783267763350693
+     -0.999999974817065447355071
+     -0.999999972689661120580137
+     -0.999999970388298708114035
+     -0.999999967899251489096719
+     -0.999999965207751402768258
+     -0.999999962297913456009131
+     -0.999999959152654912059634
+     -0.999999955753608920567185
+     -0.999999952081032228479841
+     -0.999999948113706589527419
+     -0.999999943828833467051722
+     -0.999999939201921600705774
+     -0.999999934206666981977303
+     -0.999999928814824756540487
+     -0.999999922996072543035871
+     -0.999999916717864627952817
+     -0.999999909945276464770342
+     -0.999999902640838872326778
+     -0.999999894764361292459493
+     -0.999999886272743430203347
+     -0.999999877119774561178120
+     -0.999999867255919750145400
+     -0.999999856628092181986012
+     -0.999999845179410761448597
+     -0.999999832848942090853798
+     -0.999999819571425885409310
+     -0.999999805276982833797809
+     -0.999999789890803857138683
+     -0.999999773332819662188394
+     -0.999999755517349424622676
+     -0.999999736352727375323111
+     -0.999999715740905996653838
+     -0.999999693577034466641037
+     -0.999999669749010916634964
+     -0.999999644137006992314701
+     -0.999999616612963128659530
+     -0.999999587040052866624547
+     -0.999999555272114452585247
+     -0.999999521153047871016725
+     -0.999999484516175366204968
+     -0.999999445183563409904602
+     -0.999999402965303968610622
+     -0.999999357658752816349255
+     -0.999999309047722526460751
+     -0.999999256901627658587254
+     -0.999999200974579534832405
+     -0.999999141004427871663589
+     -0.999999076711746401418634
+     -0.999999007798759479089598
+     -0.999998933948206526218906
+     -0.999998854822141014087435
+     -0.999998770060660532728832
+     -0.999998679280564330496742
+     -0.999998582073934540768369
+     -0.999998478006637137714638
+     -0.999998366616738481729930
+     -0.999998247412833126918517
+     -0.999998119872278367806604
+     -0.999997983439330800015059
+     -0.999997837523179959816750
+     -0.999997681495873890143584
+     -0.999997514690131255533768
+     -0.999997336397034395554135
+     -0.999997145863597465233101
+     -0.999996942290203561838538
+     -0.999996724827904479777092
+     -0.999996492575576469327953
+     -0.999996244576925100210857
+     -0.999995979817332047487635
+     -0.999995697220536324878170
+     -0.999995395645142189112255
+     -0.999995073880945628324028
+     -0.999994730645071027620088
+     -0.999994364577909275721057
+     -0.999993974238848237905028
+     -0.999993558101786172297790
+     -0.999993114550418308799710
+     -0.999992641873286442566995
+     -0.999992138258581016945169
+     -0.999991601788684784070191
+     -0.999991030434446735009890
+     -0.999990422049174585335565
+     -0.999989774362333686430315
+     -0.999989084972939807716761
+     -0.999988351342632800403966
+     -0.999987570788417709415870
+     -0.999986740475059447000171
+     -0.999985857407116679281056
+     -0.999984918420600105892450
+     -0.999983920174239833019164
+     -0.999982859140346051919802
+     -0.999981731595246738576625
+     -0.999980533609285585815045
+     -0.999979261036362867394040
+     -0.999977909503001414558627
+     -0.999976474396919359773184
+     -0.999974950855090770264459
+     -0.999973333751274756079605
+     -0.999971617682993094132439
+     -0.999969796957935861737626
+     -0.999967865579774021028779
+     -0.999965817233357340078774
+     -0.999963645269275478192029
+     -0.999961342687759502470105
+     -0.999958902121900541164316
+     -0.999956315820161717373451
+     -0.999953575628158945224767
+     -0.999950672969685610748981
+     -0.999947598826955602236751
+     -0.999944343720038601012541
+     -0.999940897685460994406321
+     -0.999937250253945229426789
+     -0.999933390427259889483441
+     -0.999929306654152248770470
+     -0.999924986805334540975777
+     -0.999920418147494672237463
+     -0.999915587316301614224693
+     -0.999910480288375233419033
+     -0.999905082352189848729210
+     -0.999899378077880363163096
+     -0.999893351285919388145981
+     -0.999886985014633373018311
+     -0.999880261486525369137321
+     -0.999873162073371699777589
+     -0.999865667260059475670860
+     -0.999857756607131593604098
+     -0.999849408712005584127087
+     -0.999840601168832436289035
+     -0.999831310526961324669729
+     -0.999821512247975999096039
+     -0.999811180661268472697860
+     -0.999800288918115561773568
+     -0.999788808944223793772524
+     -0.999776711390708210081766
+     -0.999763965583470650796009
+     -0.999750539470943221871744
+     -0.999736399570162813678909
+     -0.999721510911142766670001
+     -0.999705836979508067426713
+     -0.999689339657360809492885
+     -0.999671979162343070964282
+     -0.999653713984864847611602
+     -0.999634500823465239215587
+     -0.999614294518275720661921
+     -0.999593047982555041060436
+     -0.999570712132266086606323
+     -0.999547235813665918980236
+     -0.999522565728881163674247
+     -0.999496646359441974605568
+     -0.999469419887748945596554
+     -0.999440826116448578590456
+     -0.999410802385694255639232
+     -0.999379283488271099505731
+     -0.999346201582564648884856
+     -0.999311486103354921430255
+     -0.999275063670419193575791
+     -0.999236857994928693095374
+     -0.999196789783626380902110
+     -0.999154776640775095111632
+     -0.999110732967867545150889
+     -0.999064569861091978842036
+     -0.999016195006549802974103
+     -0.998965512573224019809696
+     -0.998912423103700050040183
+     -0.998856823402643348534653
+     -0.998798606423041184302075
+     -0.998737661150219051712021
+     -0.998673872483645407328131
+     -0.998607121116541786642636
+     -0.998537283413318848302089
+     -0.998464231284862520615845
+     -0.998387832061698186549948
+     -0.998307948365064739107359
+     -0.998224437975934368814823
+     -0.998137153702018108556548
+     -0.998045943242801457529124
+     -0.997950649052658834678793
+     -0.997851108202100171372239
+     -0.997747152237207641699493
+     -0.997638607037325344858751
+     -0.997525292671069695356701
+     -0.997407023250733340776017
+     -0.997283606785160610778704
+     -0.997154845031177801648948
+     -0.997020533343667014496115
+     -0.996880460524377788338765
+     -0.996734408669576397404104
+     -0.996582153016638410456033
+     -0.996423461789695933775340
+     -0.996258096044456873132539
+     -0.996085809512319547895015
+     -0.995906348443912066070672
+     -0.995719451452192015978668
+     -0.995524849355248241312063
+     -0.995322265018952734162069
+     -0.995111413199616997238347
+     -0.994892000386813583757223
+     -0.994663724646529912575485
+     -0.994426275464827868000969
+     -0.994179333592189118776836
+     -0.993922570888732519066544
+     -0.993655650170496375504198
+     -0.993378225056984767759469
+     -0.993089939820183484438217
+     -0.992790429235257469948358
+     -0.992479318433147959270642
+     -0.992156222755293694076721
+     -0.991820747610706752658708
+     -0.991472488335639574649461
+     -0.991111030056085706155210
+     -0.990735947553362618000491
+     -0.990346805133030645297762
+     -0.989943156497407646205095
+     -0.989524544621944366953582
+     -0.989090501635730714183733
+     -0.988640548706408159263284
+     -0.988174195929768317289394
+     -0.987690942224322340437041
+     -0.987190275231130125566791
+     -0.986671671219182443772211
+     -0.986134594996632938939331
+     -0.985578499828180497431188
+     -0.985002827358905745666562
+     -0.984407007544868370574912
+     -0.983790458590774563626243
+     -0.983152586895026146341339
+     -0.982492787002464827943096
+     -0.981810441565126558280198
+     -0.981104921311322055650339
+     -0.980375585023360294162452
+     -0.979621779524232013515387
+     -0.978842839673570150500483
+     -0.978038088373203471420329
+     -0.977206836582618593554779
+     -0.976348383344644007774283
+     -0.975462015821667639794134
+     -0.974547009342696901235481
+     -0.973602627461567070965916
+     -0.972628122026600200033799
+     -0.971622733262012538372566
+     -0.970585689861363727448927
+     -0.969516209093335679948252
+     -0.968413496920123165931304
+     -0.967276748128711635913826
+     -0.966105146475310727066976
+     -0.964897864843204212102907
+     -0.963654065414268855166142
+     -0.962372899854405729886811
+     -0.961053509513118027313721
+     -0.959695025637459232377443
+     -0.958296569600564774513511
+     -0.956857253144968859850497
+     -0.955376178640896168788291
+     -0.953852439359705454847226
+     -0.952285119762648810516483
+     -0.950673295805096476827758
+     -0.949016035256362570111835
+     -0.947312398035251987699776
+     -0.945561436561433041323207
+     -0.943762196122724061065830
+     -0.941913715258365323657307
+     -0.940015026158330197754595
+     -0.938065155078711378050909
+     -0.936063122773199513379509
+     -0.934007944940652436603893
+     -0.931898632688733592554243
+     -0.929734193013578152514172
+     -0.927513629295424719100128
+     -0.925235941810129484510736
+     -0.922900128256458230136523
+     -0.920505184299029669863949
+     -0.918050104126761367892733
+     -0.915533881026646830667193
+     -0.912955507972669409081356
+     -0.910313978229635380238406
+     -0.907608285971685037790485
+     -0.904837426915216837567545
+     -0.902000398965935653353894
+     -0.899096202879712030000041
+     -0.896123842936915012419506
+     -0.893082327629856715369185
+     -0.889970670362962317254718
+     -0.886787890165254649549837
+     -0.883533012414718050732089
+     -0.880205069574081699771868
+     -0.876803101937538279856830
+     -0.873326158387889589649984
+     -0.869773297163586659132690
+     -0.866143586635108082774537
+     -0.862436106090096697765226
+     -0.858649946526651453191925
+     -0.854784211454148381612389
+     -0.850838017700942042052444
+     -0.846810496228276697886304
+     -0.842700792949714869341221
+     -0.838508069555369803579790
+     -0.834231504340207880514292
+     -0.829870293035667064651115
+     -0.825423649643818267103596
+     -0.820890807273277941907934
+     -0.816271018976062450999077
+     -0.811563558584557736234438
+     -0.806767721547761677691787
+     -0.801882825765941237429626
+     -0.796908212422832128518725
+     -0.791843246814495351087002
+     -0.786687319173932540042657
+     -0.781439845490550711651092
+     -0.776100268323556715663465
+     -0.770668057608352532380083
+     -0.765142711454994534663544
+     -0.759523756937772996814901
+     -0.753810750874962507155164
+     -0.748003280597789559230994
+     -0.742100964707660486167111
+     -0.736103453820691090374864
+     -0.730010431298578831244359
+     -0.723821613964859289254273
+     -0.717536752805590846718625
+     -0.711155633653515131598938
+     -0.704678077854745777673723
+     -0.698103942917044475821896
+     -0.691433123138751138842000
+     -0.684665550217444285287115
+     -0.677801193837418472975629
+     -0.670840062235077784356887
+     -0.663782202741357985424868
+     -0.656627702300305046439084
+     -0.649376687962954222750275
+     -0.642029327355671841292946
+     -0.634585829122141313859706
+     -0.627046443338195690464736
+     -0.619411461898721260562292
+     -0.611681218875880289504687
+     -0.603856090847925922562622
+     -0.595936497197908576518189
+     -0.587922900381600743987037
+     -0.579815806163996030236782
+     -0.571615763823768395668684
+     -0.563323366325108955991688
+     -0.554939250456390260497539
+     -0.546464096935141688222533
+     -0.537898630478854431106614
+     -0.529243619841170429010070
+     -0.520499877813046537682747
+     -0.511668261188523099211697
+     -0.502749670694764894533266
+     -0.493745050886082136517092
+     -0.484655390001679655026703
+     -0.475481719786923675319682
+     -0.466225115277957539319705
+     -0.456886694549540304987442
+     -0.447467618426025319819411
+     -0.437969090155439536373130
+     -0.428392355046668455103604
+     -0.418738700069796075759558
+     -0.409009453419694044918239
+     -0.399205984042999234582144
+     -0.389329701128664204753628
+     -0.379382053562310320297300
+     -0.369364529344658637204124
+     -0.359278654974358972988254
+     -0.349125994795582757315695
+     -0.338908150310790248699499
+     -0.328626759459127427638914
+     -0.318283495860952262338992
+     -0.307880068029034016649000
+     -0.297418218547012757893882
+     -0.286899723215749152048332
+     -0.276326390168236932985068
+     -0.265700058953792029795712
+     -0.255022599592273160770179
+     -0.244295911599128684869262
+     -0.233521922982103572925263
+     -0.222702589210478454140139
+     -0.211839892157749742936947
+     -0.200935839018695792648413
+     -0.189992461201808793395082
+     -0.179011813198105671155184
+     -0.167995971427363494770210
+     -0.156947033062855800871902
+     -0.145867114835695749082958
+     -0.134758351819920069295014
+     -0.123622896199474309769839
+     -0.112462916018284892203275
+     -0.101280593914626883352498
+     -0.0900781258410181607233922
+     -0.0788577197708907433569970
+     -0.0676215943933084420794315
+     -0.0563719777970166238312711
+     -0.0451111061451247520897491
+     -0.0338412223417354333022167
+     -0.0225645746918449442243659
+     -0.0112834155558496169159095
+     0.
+     0.0112834155558496169159095
+     0.0225645746918449442243659
+     0.0338412223417354333022167
+     0.0451111061451247520897491
+     0.0563719777970166238312711
+     0.0676215943933084420794315
+     0.0788577197708907433569970
+     0.0900781258410181607233922
+     0.101280593914626883352498
+     0.112462916018284892203275
+     0.123622896199474309769839
+     0.134758351819920069295014
+     0.145867114835695749082958
+     0.156947033062855800871902
+     0.167995971427363494770210
+     0.179011813198105671155184
+     0.189992461201808793395082
+     0.200935839018695792648413
+     0.211839892157749742936947
+     0.222702589210478454140139
+     0.233521922982103572925263
+     0.244295911599128684869262
+     0.255022599592273160770179
+     0.265700058953792029795712
+     0.276326390168236932985068
+     0.286899723215749152048332
+     0.297418218547012757893882
+     0.307880068029034016649000
+     0.318283495860952262338992
+     0.328626759459127427638914
+     0.338908150310790248699499
+     0.349125994795582757315695
+     0.359278654974358972988254
+     0.369364529344658637204124
+     0.379382053562310320297300
+     0.389329701128664204753628
+     0.399205984042999234582144
+     0.409009453419694044918239
+     0.418738700069796075759558
+     0.428392355046668455103604
+     0.437969090155439536373130
+     0.447467618426025319819411
+     0.456886694549540304987442
+     0.466225115277957539319705
+     0.475481719786923675319682
+     0.484655390001679655026703
+     0.493745050886082136517092
+     0.502749670694764894533266
+     0.511668261188523099211697
+     0.520499877813046537682747
+     0.529243619841170429010070
+     0.537898630478854431106614
+     0.546464096935141688222533
+     0.554939250456390260497539
+     0.563323366325108955991688
+     0.571615763823768395668684
+     0.579815806163996030236782
+     0.587922900381600743987037
+     0.595936497197908576518189
+     0.603856090847925922562622
+     0.611681218875880289504687
+     0.619411461898721260562292
+     0.627046443338195690464736
+     0.634585829122141313859706
+     0.642029327355671841292946
+     0.649376687962954222750275
+     0.656627702300305046439084
+     0.663782202741357985424868
+     0.670840062235077784356887
+     0.677801193837418472975629
+     0.684665550217444285287115
+     0.691433123138751138842000
+     0.698103942917044475821896
+     0.704678077854745777673723
+     0.711155633653515131598938
+     0.717536752805590846718625
+     0.723821613964859289254273
+     0.730010431298578831244359
+     0.736103453820691090374864
+     0.742100964707660486167111
+     0.748003280597789559230994
+     0.753810750874962507155164
+     0.759523756937772996814901
+     0.765142711454994534663544
+     0.770668057608352532380083
+     0.776100268323556715663465
+     0.781439845490550711651092
+     0.786687319173932540042657
+     0.791843246814495351087002
+     0.796908212422832128518725
+     0.801882825765941237429626
+     0.806767721547761677691787
+     0.811563558584557736234438
+     0.816271018976062450999077
+     0.820890807273277941907934
+     0.825423649643818267103596
+     0.829870293035667064651115
+     0.834231504340207880514292
+     0.838508069555369803579790
+     0.842700792949714869341221
+     0.846810496228276697886304
+     0.850838017700942042052444
+     0.854784211454148381612389
+     0.858649946526651453191925
+     0.862436106090096697765226
+     0.866143586635108082774537
+     0.869773297163586659132690
+     0.873326158387889589649984
+     0.876803101937538279856830
+     0.880205069574081699771868
+     0.883533012414718050732089
+     0.886787890165254649549837
+     0.889970670362962317254718
+     0.893082327629856715369185
+     0.896123842936915012419506
+     0.899096202879712030000041
+     0.902000398965935653353894
+     0.904837426915216837567545
+     0.907608285971685037790485
+     0.910313978229635380238406
+     0.912955507972669409081356
+     0.915533881026646830667193
+     0.918050104126761367892733
+     0.920505184299029669863949
+     0.922900128256458230136523
+     0.925235941810129484510736
+     0.927513629295424719100128
+     0.929734193013578152514172
+     0.931898632688733592554243
+     0.934007944940652436603893
+     0.936063122773199513379509
+     0.938065155078711378050909
+     0.940015026158330197754595
+     0.941913715258365323657307
+     0.943762196122724061065830
+     0.945561436561433041323207
+     0.947312398035251987699776
+     0.949016035256362570111835
+     0.950673295805096476827758
+     0.952285119762648810516483
+     0.953852439359705454847226
+     0.955376178640896168788291
+     0.956857253144968859850497
+     0.958296569600564774513511
+     0.959695025637459232377443
+     0.961053509513118027313721
+     0.962372899854405729886811
+     0.963654065414268855166142
+     0.964897864843204212102907
+     0.966105146475310727066976
+     0.967276748128711635913826
+     0.968413496920123165931304
+     0.969516209093335679948252
+     0.970585689861363727448927
+     0.971622733262012538372566
+     0.972628122026600200033799
+     0.973602627461567070965916
+     0.974547009342696901235481
+     0.975462015821667639794134
+     0.976348383344644007774283
+     0.977206836582618593554779
+     0.978038088373203471420329
+     0.978842839673570150500483
+     0.979621779524232013515387
+     0.980375585023360294162452
+     0.981104921311322055650339
+     0.981810441565126558280198
+     0.982492787002464827943096
+     0.983152586895026146341339
+     0.983790458590774563626243
+     0.984407007544868370574912
+     0.985002827358905745666562
+     0.985578499828180497431188
+     0.986134594996632938939331
+     0.986671671219182443772211
+     0.987190275231130125566791
+     0.987690942224322340437041
+     0.988174195929768317289394
+     0.988640548706408159263284
+     0.989090501635730714183733
+     0.989524544621944366953582
+     0.989943156497407646205095
+     0.990346805133030645297762
+     0.990735947553362618000491
+     0.991111030056085706155210
+     0.991472488335639574649461
+     0.991820747610706752658708
+     0.992156222755293694076721
+     0.992479318433147959270642
+     0.992790429235257469948358
+     0.993089939820183484438217
+     0.993378225056984767759469
+     0.993655650170496375504198
+     0.993922570888732519066544
+     0.994179333592189118776836
+     0.994426275464827868000969
+     0.994663724646529912575485
+     0.994892000386813583757223
+     0.995111413199616997238347
+     0.995322265018952734162069
+     0.995524849355248241312063
+     0.995719451452192015978668
+     0.995906348443912066070672
+     0.996085809512319547895015
+     0.996258096044456873132539
+     0.996423461789695933775340
+     0.996582153016638410456033
+     0.996734408669576397404104
+     0.996880460524377788338765
+     0.997020533343667014496115
+     0.997154845031177801648948
+     0.997283606785160610778704
+     0.997407023250733340776017
+     0.997525292671069695356701
+     0.997638607037325344858751
+     0.997747152237207641699493
+     0.997851108202100171372239
+     0.997950649052658834678793
+     0.998045943242801457529124
+     0.998137153702018108556548
+     0.998224437975934368814823
+     0.998307948365064739107359
+     0.998387832061698186549948
+     0.998464231284862520615845
+     0.998537283413318848302089
+     0.998607121116541786642636
+     0.998673872483645407328131
+     0.998737661150219051712021
+     0.998798606423041184302075
+     0.998856823402643348534653
+     0.998912423103700050040183
+     0.998965512573224019809696
+     0.999016195006549802974103
+     0.999064569861091978842036
+     0.999110732967867545150889
+     0.999154776640775095111632
+     0.999196789783626380902110
+     0.999236857994928693095374
+     0.999275063670419193575791
+     0.999311486103354921430255
+     0.999346201582564648884856
+     0.999379283488271099505731
+     0.999410802385694255639232
+     0.999440826116448578590456
+     0.999469419887748945596554
+     0.999496646359441974605568
+     0.999522565728881163674247
+     0.999547235813665918980236
+     0.999570712132266086606323
+     0.999593047982555041060436
+     0.999614294518275720661921
+     0.999634500823465239215587
+     0.999653713984864847611602
+     0.999671979162343070964282
+     0.999689339657360809492885
+     0.999705836979508067426713
+     0.999721510911142766670001
+     0.999736399570162813678909
+     0.999750539470943221871744
+     0.999763965583470650796009
+     0.999776711390708210081766
+     0.999788808944223793772524
+     0.999800288918115561773568
+     0.999811180661268472697860
+     0.999821512247975999096039
+     0.999831310526961324669729
+     0.999840601168832436289035
+     0.999849408712005584127087
+     0.999857756607131593604098
+     0.999865667260059475670860
+     0.999873162073371699777589
+     0.999880261486525369137321
+     0.999886985014633373018311
+     0.999893351285919388145981
+     0.999899378077880363163096
+     0.999905082352189848729210
+     0.999910480288375233419033
+     0.999915587316301614224693
+     0.999920418147494672237463
+     0.999924986805334540975777
+     0.999929306654152248770470
+     0.999933390427259889483441
+     0.999937250253945229426789
+     0.999940897685460994406321
+     0.999944343720038601012541
+     0.999947598826955602236751
+     0.999950672969685610748981
+     0.999953575628158945224767
+     0.999956315820161717373451
+     0.999958902121900541164316
+     0.999961342687759502470105
+     0.999963645269275478192029
+     0.999965817233357340078774
+     0.999967865579774021028779
+     0.999969796957935861737626
+     0.999971617682993094132439
+     0.999973333751274756079605
+     0.999974950855090770264459
+     0.999976474396919359773184
+     0.999977909503001414558627
+     0.999979261036362867394040
+     0.999980533609285585815045
+     0.999981731595246738576625
+     0.999982859140346051919802
+     0.999983920174239833019164
+     0.999984918420600105892450
+     0.999985857407116679281056
+     0.999986740475059447000171
+     0.999987570788417709415870
+     0.999988351342632800403966
+     0.999989084972939807716761
+     0.999989774362333686430315
+     0.999990422049174585335565
+     0.999991030434446735009890
+     0.999991601788684784070191
+     0.999992138258581016945169
+     0.999992641873286442566995
+     0.999993114550418308799710
+     0.999993558101786172297790
+     0.999993974238848237905028
+     0.999994364577909275721057
+     0.999994730645071027620088
+     0.999995073880945628324028
+     0.999995395645142189112255
+     0.999995697220536324878170
+     0.999995979817332047487635
+     0.999996244576925100210857
+     0.999996492575576469327953
+     0.999996724827904479777092
+     0.999996942290203561838538
+     0.999997145863597465233101
+     0.999997336397034395554135
+     0.999997514690131255533768
+     0.999997681495873890143584
+     0.999997837523179959816750
+     0.999997983439330800015059
+     0.999998119872278367806604
+     0.999998247412833126918517
+     0.999998366616738481729930
+     0.999998478006637137714638
+     0.999998582073934540768369
+     0.999998679280564330496742
+     0.999998770060660532728832
+     0.999998854822141014087435
+     0.999998933948206526218906
+     0.999999007798759479089598
+     0.999999076711746401418634
+     0.999999141004427871663589
+     0.999999200974579534832405
+     0.999999256901627658587254
+     0.999999309047722526460751
+     0.999999357658752816349255
+     0.999999402965303968610622
+     0.999999445183563409904602
+     0.999999484516175366204968
+     0.999999521153047871016725
+     0.999999555272114452585247
+     0.999999587040052866624547
+     0.999999616612963128659530
+     0.999999644137006992314701
+     0.999999669749010916634964
+     0.999999693577034466641037
+     0.999999715740905996653838
+     0.999999736352727375323111
+     0.999999755517349424622676
+     0.999999773332819662188394
+     0.999999789890803857138683
+     0.999999805276982833797809
+     0.999999819571425885409310
+     0.999999832848942090853798
+     0.999999845179410761448597
+     0.999999856628092181986012
+     0.999999867255919750145400
+     0.999999877119774561178120
+     0.999999886272743430203347
+     0.999999894764361292459493
+     0.999999902640838872326778
+     0.999999909945276464770342
+     0.999999916717864627952817
+     0.999999922996072543035871
+     0.999999928814824756540487
+     0.999999934206666981977303
+     0.999999939201921600705774
+     0.999999943828833467051722
+     0.999999948113706589527419
+     0.999999952081032228479841
+     0.999999955753608920567185
+     0.999999959152654912059634
+     0.999999962297913456009131
+     0.999999965207751402768258
+     0.999999967899251489096719
+     0.999999970388298708114035
+     0.999999972689661120580137
+     0.999999974817065447355071
+     0.999999976783267763350693
+     0.999999978600119594788978
+     0.999999980278629704073472
+     0.999999981829021830014495
+     0.999999983260788635479157
+     0.999999984582742099719981
+     0.999999985803060578628974
+     0.999999986929332742926881
+     0.999999987968598591791738
+     0.999999988927387727619661
+     0.999999989811755066458970
+     0.999999990627314148132540
+     0.999999991379268200130508
+     0.999999992072439099985541
+     0.999999992711294372006456
+     0.999999993299972345915102
+     0.999999993842305597079599
+     0.999999994341842780638679
+     0.999999994801868964843021
+     0.999999995225424562376994
+     0.999999995615322952246180
+     0.999999995974166879001692
+     0.999999996304363710601733
+     0.999999996608139631065173
+     0.999999996887552839233372
+     0.999999997144505820407811
+     0.999999997380756753356393
+     0.999999997597930111164948
+     0.999999997797526510638154
+     0.999999997980931861411821
+     0.999999998149425862613257
+     0.999999998304189891785759
+     0.999999998446314327865434
+     0.999999998576805347252390
+     0.999999998696591229443244
+     0.999999998806528206277959
+     0.999999998907405886591731
+     0.999999998999952285943071
+     0.999999999084838489103913
+     0.999999999162682971138542
+     0.999999999234055601157792
+     0.999999999299481351206210
+     0.999999999359443731215952
+     0.999999999414387969535759
+     0.999999999464723957210342
+     0.999999999510828972939411
+     0.999999999553050204480830
+     0.999999999591707081174197
+     0.999999999627093431244556
+     0.999999999659479476596644
+     0.999999999689113676923764
+     0.999999999716224434128131
+     0.999999999741021667277623
+     0.999999999763698267603729
+     0.999999999784431442373954
+     0.999999999803383955845711
+     0.999999999820705274925103
+     0.999999999836532626610029
+     0.999999999850991973790366
+     0.999999999864198915505960
+     0.999999999876259517323697
+     0.999999999887271077085832
+     0.999999999897322830900966
+     0.999999999906496603894840
+     0.999999999914867409908580
+     0.999999999922504004025582
+     0.999999999929469391523360
+     0.999999999935821296581860
+     0.999999999941612593833703
+     0.999999999946891705613252
+     0.999999999951702967549089
+     0.999999999956086964947412
+     0.999999999960080842230859
+     0.999999999963718587527495
+     0.999999999967031294347156
+     0.999999999970047402136203
+     0.999999999972792917366275
+     0.999999999975291616686970
+     0.999999999977565233555991
+     0.999999999979633629652374
+     0.999999999981514952278515
+     0.999999999983225778864110
+     0.999999999984781249599477
+     0.999999999986195189146334
+     0.999999999987480218300757
+     0.999999999988647856415078
+     0.999999999989708615322670
+     0.999999999990672085451468
+     0.999999999991547014758366
+     0.999999999992341381066975
+     0.999999999993062458345374
+     0.999999999993716877418075
+     0.999999999994310681567322
+     0.999999999994849377442675
+     0.999999999995337981664472
+     0.999999999995781063475994
+     0.999999999996182783770717
+     0.999999999996546930794861
+     0.999999999996876952801260
+     0.999999999997175987908301
+     0.999999999997446891397149
+     0.999999999997692260661541
+     0.999999999997914458007001
+     0.999999999998115631480261
+     0.999999999998297733894900
+     0.999999999998462540205572]';
+endfunction
+[Xr,Yr]=erf_ref();
+y=erf(Xr);
+if norm((Yr-y),1) >1d-12 then bugmes();quit;end
diff --git a/scilab/modules/special_functions/tests/unit_tests/erfinv.dia.ref b/scilab/modules/special_functions/tests/unit_tests/erfinv.dia.ref
new file mode 100644 (file)
index 0000000..01a3969
--- /dev/null
@@ -0,0 +1,518 @@
+// =============================================================================
+// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
+// Copyright (C) ????-2008 - INRIA Serge Steer
+//
+//  This file is distributed under the same license as the Scilab package.
+// =============================================================================
+function [x,y]=erfinv_ref()
+//reference Yr=erfinv(X) generee par Maple
+//pm:=array(1..1999):x:=-0.999:
+//for k from 1  to 1999 do
+//pm[k]:=evalf[24](RootOf(-erf(_Z)+x)):
+//x:=x+0.001:end do;
+  x=-0.999:0.001:0.999;
+  y=[-2.32675376551352467056022;-2.18512421913300426570596;-2.09850761179367089535179;-2.03516768306608309650651;
+     -1.98487261261553223101737;-1.94297485086632095708177;-1.90695686467094561133550;-1.87529654534774725251058;
+     -1.84700119608962619338729;-1.82138636771844967304021;-1.79795957771209517163847;-1.77635428959895142455381;
+     -1.75629011016036985182043;-1.73754760637571012247642;-1.71995172674121274921475;-1.70336051775180604164031;
+     -1.68765722376985222676591;-1.67274461920355980864179;-1.65854085494046514649335;-1.64497635713318705017720;
+     -1.63199147310755983780432;-1.61953465784330182137717;-1.60756105826610791904552;-1.59603139479067878351924;
+     -1.58491106805948089387746;-1.57416943844147646813481;-1.56377923959071921153068;-1.55371609713170636373557;
+     -1.54395813058290834226583;-1.53448562177771813914331;-1.52528073684951304945905;-1.51632729169521713246473;
+     -1.50761055298380584861693;-1.49911706841829737690302;-1.49083452122409072994342;-1.48275160481818015318184;
+     -1.47485791438202173394033;-1.46714385266652228379940;-1.45960054783849765418572;-1.45221978156224685014342;
+     -1.44499392581891844119832;-1.43791588721631257139615;-1.43097905774508171768046;-1.42417727110354626495924;
+     -1.41750476384996789516378;-1.41095614075394083115303;-1.40452634381214160594658;-1.39821062447163920701438;
+     -1.39200451866918531861457;-1.38590382434967794527797;-1.37990458117316776469996;-1.37400305215884578004754;
+     -1.36819570704762188366770;-1.36247920719316111112837;-1.35685039181539261030877;-1.35130626547120660708750;
+     -1.34584398661485169413881;-1.34046085713588930726969;-1.33515431277582849530455;-1.32992191433606380401593;
+     -1.32476133959973289089009;-1.31967037589881690640162;-1.31464691326540954816980;-1.30968893811273510877283;
+     -1.30479452739733301983674;-1.29996184321895768212866;-1.29518912781926248355837;-1.29047469894432663626360;
+     -1.28581694553961189638512;-1.28121432374906258624426;-1.27666535319283750241650;-1.27216861350063008607206;
+     -1.26772274107973147711542;-1.26332642609895248768275;-1.25897840967127347229233;-1.25467748121966017307579;
+     -1.25042247601189031583104;-1.24621227285149871837430;-1.24204579191308427645764;-1.23792199271124470601816;
+     -1.23383987219332687396791;-1.22979846294701294800680;-1.22579683151451512054446;-1.22183407680583284397770;
+     -1.21790932860414388392794;-1.21402174615696070959820;-1.21017051684719268883733;-1.20635485493871746442485;
+     -1.20257400039148639653512;-1.19882721774157318917829;-1.19511379504192544842046;-1.19143304285989922198847;
+     -1.18778429332794946696112;-1.18416689924411749800915;-1.18058023321920212215387;-1.17702368686772646342933;
+     -1.17349667004001930287996;-1.16999861009291979542268;-1.16652895119678919456914;-1.16308715367667408672625;
+     -1.15967269338561384167854;-1.15628506110822163735251;-1.15292376199279451989056;-1.14958831501032442296701;
+     -1.14627825243888971549360;-1.14299311937200641942412;-1.13973247324961041494949;-1.13649588341042734231966;
+     -1.13328293066456607513313;-1.13009320688524508602561;-1.12692631461862921361960;-1.12378186671081768969028;
+     -1.12065948595108318104233;-1.11755880473051639137566;-1.11447946471528177395617;-1.11142111653373741828454;
+     -1.10838341947671646039726;-1.10536604121030867162096;-1.10236865750051942874064;-1.09939095194921926524473;
+     -1.09643261574083083726659;-1.09349334739923158235067;-1.09057285255437976350432;-1.08767084371819912164160;
+     -1.08478704006928314130699;-1.08192116724600409166043;-1.07907295714763465155357;-1.07624214774311116967833;
+     -1.07342848288708754562668;-1.07063171214294743521462;-1.06785159061246006671034;-1.06508787877178148051829;
+     -1.06234034231351854451039;-1.05960875199458771638643;-1.05689288348961428416044;-1.05419251724963077262264;
+     -1.05150743836584540982523;-1.04883743643826305191789;-1.04618230544895181216668;-1.04354184363975887267130;
+     -1.04091585339428861413581;-1.03830414112396531630198;-1.03570651715801129305399;-1.03312279563717946315342;
+     -1.03055279441108704931261;-1.02799633493900437212031;-1.02545324219395958661019;-1.02292334457002672171526;
+     -1.02040647379267054860440;-1.01790246483202764360873;-1.01541115581900854441199;-1.01293238796411114243671;
+     -1.01046600547884042676895;-1.00801185549963441130450;-1.00556978801420055183020;-1.00313965579017120729917;
+     -1.00072131430599073256377;-0.998314621683950620428471;-0.995919438625292750450167;-0.993535628347304261113059;
+     -0.991163056522330850848242;-0.988801591218638441249043;-0.986451102843056111574170;-0.984111464085336045518998;
+     -0.981782549864168927059146;-0.979464237274795789251756;-0.977156405538159765104812;-0.974858935951543519472376;
+     -0.972571711840640361501646;-0.970294618513009154187884;-0.968027543212865156493316;-0.965770375077160859351606;
+     -0.963523005092912714513298;-0.961285326055731409118287;-0.959057232529515013380354;-0.956838620807265927866568;
+     -0.954629388872994084354930;-0.952429436364670313748928;-0.950238664538195189410099;-0.948056976232349987744631;
+     -0.945884275834697682976616;-0.943720469248403112623036;-0.941565463859942616962252;-0.939419168507674572321134;
+     -0.937281493451243306720676;-0.935152350341789909611329;-0.933031652192944427270628;-0.930919313352574873985593;
+     -0.928815249475269388351157;-0.926719377495528725731954;-0.924631615601647103910482;-0.922551883210260210841589;
+     -0.920480100941539942827582;-0.918416190595016169815512;-0.916360075126006523317668;-0.914311678622635873015377;
+     -0.912270926283427801704433;-0.910237744395451006093112;-0.908212060313004144225512;-0.906193802436823220071163;
+     -0.904182900193796143139771;-0.902179284017169626839285;-0.900182885327234094639705;-0.898193636512472748831602;
+     -0.896211470911161423631094;-0.894236322793406293395002;-0.892268127343606938544748;-0.890306820643332687192334;
+     -0.888352339654600550118256;-0.886404622203543451338165;-0.884463606964457826649897;-0.882529233444220018882342;
+     -0.880601441967061241651857;-0.878680173659691213822623;-0.876765370436760885091340;-0.874856974986654979676502;
+     -0.872954930757605380467648;-0.871059181944116660638341;-0.869169673473695344084982;-0.867286350993874740539321;
+     -0.865409160859527456214132;-0.863538050120457926759686;-0.861672966509267556497199;-0.859813858429485276701938;
+     -0.857960674943956556465451;-0.856113365763484112691231;-0.854271881235713771374789;-0.852436172334259130778193;
+     -0.850606190648058868708733;-0.848781888370960721117400;-0.846963218291526337899961;-0.845150133783051394355273;
+     -0.843342588793795503465446;-0.841540537837416635233851;-0.839743935983604904963956;-0.837952738848910742789354;
+     -0.836166902587762602169693;-0.834386383883669505636634;-0.832611139940603861988961;-0.830841128474560120569394;
+     -0.829076307705284955373421;-0.827316636348174794701583;-0.825562073606336631023609;-0.823812579162808160822003;
+     -0.822068113172933415564498;-0.820328636256890152754014;-0.818594109492365380350606;-0.816864494407375488876719;
+     -0.815139752973227563323866;-0.813419847597618541690289;-0.811704741117868978705118;-0.809994396794288262140596;
+     -0.808288778303668215181754;-0.806587849732902101709751;-0.804891575572726132154983;-0.803199920711580645878829;
+     -0.801512850429588221935367;-0.799830330392646043629732;-0.798152326646629913608153;-0.796478805611707385363143;
+     -0.794809734076757544089917;-0.793145079193895034858203;-0.791484808473095999135830;-0.789828889776923641882848;
+     -0.788177291315351210791011;-0.786529981640680226834519;-0.784886929642551861182764;-0.783248104543049407761361;
+     -0.781613475891889853388483;-0.779983013561702598512212;-0.778356687743393431181906;-0.776734468941591905051428;
+     -0.775116327970180318981467;-0.773502235947902541227630;-0.771892164294050965314105;-0.770286084724229927541830;
+     -0.768683969246193957705777;-0.767085790155759275037507;-0.765491520032786980684311;-0.763901131737236436221627;
+     -0.762314598405287354806035;-0.760731893445529167646167;-0.759152990535216263530909;-0.757577863616587733240031;
+     -0.756006486893250283802416;-0.754438834826623019790749;-0.752874882132442820177428;-0.751314603777329069751936;
+     -0.749757974975406533741602;-0.748204971184985193111121;-0.746655568105295886066185;-0.745109741673280628576998;
+     -0.743567468060436513291432;-0.742028723669712112047386;-0.740493485132455332341197;-0.738961729305411702584492;
+     -0.737433433267772084805852;-0.735908574318268836645561;-0.734387129972319467070368;-0.732869077959216852218817;
+     -0.731354396219365099193895;-0.729843062901560166465567;-0.728335056360314369847696;-0.726830355153223922787785;
+     -0.725328938038378678969375;-0.723830783971813263990758;-0.722335872104998801164239;-0.720844181782374454291622;
+     -0.719355692538918027627315;-0.717870384097754880153633;-0.716388236367804427776112;-0.714909229441463523112322;
+     -0.713433343592326018207544;-0.711960559272937830776354;-0.710490857112586849451786;-0.709024217915127028034019;
+     -0.707560622656836032879092;-0.706100052484305821364947;-0.704642488712365542827071;-0.703187912822036166478639;
+     -0.701736306458516253629577;-0.700287651429198304004255;-0.698841929701715118137356;-0.697399123402015629710152;
+     -0.695959214812469673283114;-0.694522186370001164193449;-0.693088020664249178425344;-0.691656700435756431033936;
+     -0.690228208574184662218446;-0.688802528116556450402505;-0.687379642245522981697320;-0.685959534287657314902413;
+     -0.684542187711772690745695;-0.683127586127265443385831;-0.681715713282482081301025;-0.680306553063110113575629;
+     -0.678900089490592206274711;-0.677496306720563262072639;-0.676095189041310024580015;-0.674696720872252816899203;
+     -0.673300886762449031837214;-0.671907671389117998920711;-0.670517059556186860896197;-0.669129036192857099763483;
+     -0.667743586352191359587010;-0.666360695209720220361604;-0.664980348062068584081204;-0.663602530325601340875000;
+     -0.662227227535087989639255;-0.660854425342385894008748;-0.659484109515141860783048;-0.658116265935511734053255;
+     -0.656750880598897704268154;-0.655387939612703037338166;-0.654027429195103934604611;-0.652669335673838240103679;
+     -0.651313645485010717032493;-0.649960345171914620681705;-0.648609421383869300338197;-0.647260860875073567785696;
+     -0.645914650503474575043100;-0.644570777229651948883019;-0.643229228115716934468944;-0.641889990324226305141360;
+     -0.640553051117110799973398;-0.639218397854617855207882;-0.637886017994268400082143;-0.636555899089827491847212;
+     -0.635228028790288568996146;-0.633902394838871105834534;-0.632578985072031455556846;-0.631257787418486672937314;
+     -0.629938789898251111605494;-0.628621980621685594656557;-0.627307347788558961046672;-0.625994879687121793846348;
+     -0.624684564693192139971279;-0.623376391269253034482742;-0.622070347963561645949760;-0.620766423409269862694708;
+     -0.619464606323556143004491;-0.618164885506768455582459;-0.616867249841578139643369;-0.615571688292144517116603;
+     -0.614278189903290092422814;-0.612986743799686178227838;-0.611697339185048788456327;-0.610409965341344642667636;
+     -0.609124611628007128659254;-0.607841267481162072869936;-0.606559922412863170806831;-0.605280566010336932319637;
+     -0.604003187935236999091285;-0.602727777922907694210153;-0.601454325781656666134284;-0.600182821392036491754959;
+     -0.598913254706135105615973;-0.597645615746874924647529;-0.596379894607320540030513;-0.595116081449994850019300;
+     -0.593854166506203509719984;-0.592594140075367575947128;-0.591335992524364227366604;-0.590079714286875442175869;
+     -0.588825295862744517576917;-0.587572727817340317262067;-0.586322000780929135059570;-0.585073105448054064775488;
+     -0.583826032576921768121352;-0.582580772988796534434417;-0.581337317567401527679750;-0.580095657258327117971639;
+     -0.578855783068446196566639;-0.577617686065336374962660;-0.576381357376708970388645;-0.575146788189844681588133;
+     -0.573913969751035860388153;-0.572682893365035286103023;-0.571453550394511351351370;-0.570225932259509569364741;
+     -0.569000030436920314338033;-0.567775836459952707816332;-0.566553341917614565530137;-0.565332538454198320481909;
+     -0.564113417768772839452082;-0.562895971614681051432447;-0.561680191799043307809948;-0.560466070182266395414699;
+     -0.559253598677558124813090;-0.558042769250447417470635;-0.556833573918309816630205;-0.555626004749898347949980;
+     -0.554420053864879657122263;-0.553215713433375352849737;-0.552012975675508484690151;-0.550811832860955086394373;
+     -0.549612277308500716456468;-0.548414301385601928668550;-0.547217897507952606527904;-0.546023058139055096379661;
+     -0.544829775789796075195604;-0.543638043018027089888778;-0.542447852428149706044881;-0.541259196670705204915226;
+     -0.540072068441968768462873;-0.538886460483548093183440;-0.537702365581986374335740;-0.536519776568369603114798;
+     -0.535338686317938120181512;-0.534159087749702369829439;-0.532980973826062799920226;-0.531804337552433853555395;
+     -0.530629171976871999273835;-0.529455470189707747371634;-0.528283225323181600734249;-0.527112430551083889350577;
+     -0.525943079088398438444614;-0.524775164190950020913310;-0.523608679155055545499181;-0.522443617317178932853529;
+     -0.521279972053589632360901;-0.520117736780024733298055;-0.518956904951354624591287;-0.517797470061252158114871;
+     -0.516639425641865271140692;-0.515482765263493024205197;-0.514327482534265011304754;-0.513173571099824099964556;
+     -0.512021024643012459349639;-0.510869836883560835199485;-0.509720001577781030970368;-0.508571512518261555162144;
+     -0.507424363533566395388924;-0.506278548487936880326002;-0.505134061280996591228925;-0.503990895847459285274674;
+     -0.502849046156839793519900;-0.501708506213167856807092;-0.500569270054704863476697;-0.499431331753663453261643;
+     -0.498294685415929952250655;-0.497159325180789604308388;-0.496025245220654564833729;-0.494892439740794623223041;
+     -0.493760902979070620882506;-0.492630629205670532103474;-0.491501612722848175576772;-0.490373847864664524776547;
+     -0.489247328996731585891494;-0.488122050515958812421373;-0.486998006850302025989727;-0.485875192458514813349731;
+     -0.484753601829902369979366;-0.483633229484077761074596;-0.482514069970720571155200;-0.481396117869337913897394;
+     -0.480279367789027774200510;-0.479163814368244654881929;-0.478049452274567500775245;-0.476936276204469873381418;
+     -0.475824280883092349591560;-0.474713461064017118363061;-0.473603811529044749588163;-0.472495327087973109745849;
+     -0.471388002578378399274204;-0.470281832865398286941320;-0.469176812841517116828335;-0.468072937426353163868614;
+     -0.466970201566447914212286;-0.465868600235057347005579;-0.464768128431945194489628;-0.463668781183178157633870;
+     -0.462570553540923054824727;-0.461473440583245881431235;-0.460377437413912758365590;-0.459282539162192748048363;
+     -0.458188740982662516475480;-0.457096038055012820367009;-0.456004425583856798656448;-0.454913898798540047853642;
+     -0.453824452952952461084706;-0.452736083325341810878516;-0.451648785218129056031497;-0.450562553957725353140641;
+     -0.449477384894350753649022;-0.448393273401854567498595;-0.447310214877537374731816;-0.446228204741974666626704;
+     -0.445147238438842098189415;-0.444067311434742334064264;-0.442988419219033470153526;-0.441910557303659013468255;
+     -0.440833721222979402956917;-0.439757906533605054280818;-0.438683108814230911724245;-0.437609323665472490642935;
+     -0.436536546709703394067020;-0.435464773590894287284015;-0.434393999974453314433743;-0.433324221547067941350446;
+     -0.432255434016548209087678;-0.431187633111671382759028;-0.430120814582027980522297;-0.429054974197869167726505;
+     -0.427990107749955501430070;-0.426926211049407010684775;-0.425863279927554598163639;-0.424801310235792748891791;
+     -0.423740297845433532017696;-0.422680238647561881737857;-0.421621128552892143661369;-0.420562963491625873071421;
+     -0.419505739413310871709193;-0.418449452286701449871496;-0.417394098099619900777074;-0.416339672858819174317723;
+     -0.415286172589846737469336;-0.414233593336909608794685;-0.413181931162740554624251;-0.412131182148465434653718;
+     -0.411081342393471684846923;-0.410032408015277925681117;-0.408984375149404683917372;-0.407937239949246216222914;
+     -0.406890998585943423114100;-0.405845647248257841828688;-0.404801182142446706874082;-0.403757599492139067134282;
+     -0.402714895538212948552500;-0.401673066538673551538719;-0.400632108768532472381979;-0.399592018519687938075904;
+     -0.398552792100806044092875;-0.397514425837202984767494;-0.396476916070728266073398;-0.395440259159648890699302;
+     -0.394404451478534505450229;-0.393369489418143501118366;-0.392335369385310055084850;-0.391302087802832107029012;
+     -0.390269641109360258235345;-0.389238025759287585100599;-0.388207238222640357554016;-0.387177274984969653212906;
+     -0.386148132547243858203376;-0.385119807425742045682276;-0.384092296151948223201199;-0.383065595272446440156734;
+     -0.382039701348816746673194;-0.381014610957531995364634;-0.379990320689855477522284;-0.378966827151739385371447;
+     -0.377944126963724092138613;-0.376922216760838241764837;-0.375901093192499640195609;-0.374880752922416940270228;
+     -0.373861192628492112325372;-0.372842409002723692717928;-0.371824398751110802561401;-0.370807158593557929058249;
+     -0.369790685263780461897400;-0.368774975509210977271926;-0.367760026090906262156518;-0.366745833783455071567895;
+     -0.365732395374886611613717;-0.364719707666579741216936;-0.363707767473172885482797;-0.362696571622474653754979;
+     -0.361686116955375155485562;-0.360676400325758007120738;-0.359667418600413023280381;-0.358659168658949585584833;
+     -0.357651647393710682556510;-0.356644851709687614097242;-0.355638778524435354114617;-0.354633424767988564942025;
+     -0.353628787382778257267641;-0.352624863323549089357157;-0.351621649557277299423852;-0.350619143063089265067382;
+     -0.349617340832180683769720;-0.348616239867736368502744;-0.347615837184850652567335;-0.346616129810448397848247;
+     -0.345617114783206600732730;-0.344618789153476590003681;-0.343621149983206811080181;-0.342624194345866191039560;
+     -0.341627919326368078915602;-0.340632322020994755827280;-0.339637399537322509551397;-0.338643148994147268210747;
+     -0.337649567521410787806985;-0.336656652260127388384144;-0.335664400362311233664909;-0.334672808990904149057107;
+     -0.333681875319703972982629;-0.332691596533293436535013;-0.331701969826969566525310;-0.330712992406673607028554;
+     -0.329724661488921454595204;-0.328736974300734602343382;-0.327749928079571588198463;-0.326763520073259942596795;
+     -0.325777747539928631019810;-0.324792607747940986773783;-0.323808097975828129478784;-0.322824215512222864778164;
+     -0.321840957655794060827042;-0.320858321715181497164882;-0.319876305008931181623254;-0.318894904865431130965352;
+     -0.317914118622847610998746;-0.316933943629061831947228;-0.315954377241607094911430;-0.314975416827606385291216;
+     -0.313997059763710409085611;-0.313019303436036068028315;-0.312042145240105369558602;-0.311065582580784767668659;
+     -0.310089612872224930709180;-0.309114233537800932275328;-0.308139442010052861334930;-0.307165235730626847800129;
+     -0.306191612150216499782543;-0.305218568728504748810367;-0.304246102934106099323804;-0.303274212244509278802667;
+     -0.302302894146020284917042;-0.301332146133705826128503;-0.300361965711337152205522;-0.299392350391334271152468;
+     -0.298423297694710549086888;-0.297454805151017689634687;-0.296486870298291089447284;-0.295519490682995566478939;
+     -0.294552663859971457696095;-0.293586387392381082923904;-0.292620658851655571567981;-0.291655475817442048981958;
+     -0.290690835877551179283553;-0.289726736627905061453614;-0.288763175672485475584010;-0.287800150623282476171267;
+     -0.286837659100243329383496;-0.285875698731221791258505;-0.284914267151927723820923;-0.283953362005877046135800;
+     -0.282992980944342017345393;-0.282033121626301848764835;-0.281073781718393642140920;-0.280114958894863651206581;
+     -0.279156650837518863691535;-0.278198855235678900977229;-0.277241569786128232611518;-0.276284792193068702925524;
+     -0.275328520168072367021790;-0.274372751430034633429268;-0.273417483705127710746720;-0.272462714726754355621958;
+     -0.271508442235501919439778;-0.270554663979096691116723;-0.269601377712358533425668;-0.268648581197155810297895;
+     -0.267696272202360602574682;-0.266744448503804209704503;-0.265793107884232934905753;-0.264842248133264151338477;
+     -0.263891867047342646851815;-0.262941962429697244896955;-0.261992532090297699218086;-0.261043573845811859956383;
+     -0.260095085519563108824305;-0.259147064941488061029464;-0.258199509948094531649116;-0.257252418382419764177804;
+     -0.256305788093988918991992;-0.255359616938773819496520;-0.254413902779151953738566;-0.253468643483865729295322;
+     -0.252523836927981979261970;-0.251579480992851717186640;-0.250635573566070138818935;-0.249692112541436868558283;
+     -0.248749095818916448507819;-0.247806521304599068058765;-0.246864386910661531949278;-0.245922690555328464760578;
+     -0.244981430162833749831766;-0.244040603663382200593156;-0.243100208993111462336147;-0.242160244094054142455665;
+     -0.241220706914100167219025;-0.240281595406959363132655;-0.239342907532124260995555;-0.238404641254833120745613;
+     -0.237466794546033175221892;-0.236529365382344090982915;-0.235592351746021644337605;-0.234655751624921610762049;
+     -0.233719563012463865891554;-0.232783783907596696293605;-0.231848412314761318243294;-0.230913446243856602738561;
+     -0.229978883710204005008221;-0.229044722734512696781181;-0.228110961342844899600545;-0.227177597566581417481393;
+     -0.226244629442387367225997;-0.225312055012178104725014;-0.224379872323085345587816;-0.223448079427423478459602;
+     -0.222516674382656069397254;-0.221585655251362555690039;-0.220655020101205127525302;-0.219724767004895795913128;
+     -0.218794894040163645297663;-0.217865399289722269296372;-0.216936280841237388021901;-0.216007536787294645454492;
+     -0.215079165225367585346059;-0.214151164257785804149977;-0.213223531991703279483542;-0.212296266539066872642729;
+     -0.211369366016585003701502;-0.210442828545696497740329;-0.209516652252539600760912;-0.208590835267921163856305;
+     -0.207665375727285994217637;-0.206740271770686371570609;-0.205815521542751728646714;-0.204891123192658494305804;
+     -0.203967074874100097938185;-0.203043374745257133785850;-0.202120020968767683833755;-0.201197011711697797933239;
+     -0.200274345145512129830755;-0.199352019446044727786039;-0.198430032793469978474653;-0.197508383372273702880600;
+     -0.196587069371224402895272;-0.195666088983344657349530;-0.194745440405882666216067;-0.193825121840283941729502;
+     -0.192905131492163145181806;-0.191985467571276068160731;-0.191066128291491757008872;-0.190147111870764779290812;
+     -0.189228416531107631065596;-0.188310040498563283771365;-0.187391982003177869538579;-0.186474239278973503757652;
+     -0.185556810563921243736198;-0.184639694099914182290314;-0.183722888132740675123476;-0.182806390912057700855673;
+     -0.181890200691364352574378;-0.180974315727975459787782;-0.180058734282995339669514;-0.179143454621291676492749;
+     -0.178228475011469528160158;-0.177313793725845458744699;-0.176399409040421795964603;-0.175485319234861012524278;
+     -0.174571522592460230261037;-0.173658017400125846045738;-0.172744801948348278393465;-0.171831874531176833748353;
+     -0.170919233446194691414565;-0.170006876994494006113221;-0.169094803480651127152813;-0.168183011212701933208286;
+     -0.167271498502117281711534;-0.166360263663778571863528;-0.165449305015953420285730;-0.164538620880271448335743;
+     -0.163628209581700180119424;-0.162718069448521050238859;-0.161808198812305520322682;-0.160898596007891303392273;
+     -0.159989259373358695124303;-0.159080187250007011076982;-0.158171377982331128954168;-0.157262829917998134988229;
+     -0.156354541407824073529228;-0.155446510805750798934573;-0.154538736468822928859817;-0.153631216757164898057737;
+     -0.152723950033958111799220;-0.151816934665418198035783;-0.150910169020772357429845;-0.150003651472236810385003;
+     -0.149097380394994340214742;-0.148191354167171931594021;-0.147285571169818503444195;-0.146380029786882735407659;
+     -0.145474728405190987074448;-0.144569665414425309128859;-0.143664839207101545589882;-0.142760248178547526324913;
+     -0.141855890726881349021849;-0.140951765252989749810222;-0.140047870160506561727510;-0.139144203855791260232257;
+     -0.138240764747907594970952;-0.137337551248602307010996;-0.136434561772283930757331;-0.135531794736001679775511;
+     -0.134629248559424415749181;-0.133726921664819699805000;-0.132824812477032925443096;-0.131922919423466532316157;
+     -0.131021240934059300105162;-0.130119775441265721744671;-0.129218521380035455255404;-0.128317477187792853446625;
+     -0.127416641304416570755579;-0.126516012172219246495889;-0.125615588235927263791463;-0.124715367942660583477010;
+     -0.123815349741912652250820;-0.122915532085530384369896;-0.122015913427694216181984;-0.121116492224898232793403;
+     -0.120217266935930366175908;-0.119318236021852664020104;-0.118419397945981628647141;-0.117520751173868625294633;
+     -0.116622294173280359096853;-0.115724025414179420083364;-0.114825943368704895524282;-0.113928046511153048954376;
+     -0.113030333317958065212148;-0.112132802267672860833969;-0.111235451840949959147184;-0.110338280520522429409946;
+     -0.109441286791184889349302;-0.108544469139774570452789;-0.107647826055152445372494;-0.106751356028184416804184;
+     -0.105855057551722567207705;-0.104958929120586468738427;-0.104062969231544552763038;-0.103167176383295538336451;
+     -0.102271549076449919020062;-0.101376085813511507424966;-0.100480785098859036867117;-0.0995856454387278195247507;
+     -0.0986906653411914604916251;-0.0977958433161436271229397;-0.0969011778752798730739420;-0.0960066675320795164344287;
+     -0.0951123108017875713654626;-0.0942181062013967326477174;-0.0933240522496294125539148;-0.0924301474669198294608308;
+     -0.0915363903753961476193256;-0.0906427794988626675037925;-0.0897493133627820661653203;-0.0888559904942576870157373;
+     -0.0879628094220158784725300;-0.0870697686763883808974295;-0.0861768667892947612642142;-0.0852841022942248949940076;
+     -0.0843914737262214943990342;-0.0834989796218626831784544;-0.0826066185192446164125211;-0.0817143889579641455038835;
+     -0.0808222894791015275174189;-0.0799303186252031783724914;-0.0790384749402644693440208;-0.0781467569697125663311978;
+     -0.0772551632603893113550998;-0.0763636923605341457488485;-0.0754723428197670745063028;-0.0745811131890716712576034;
+     -0.0736900020207781233421723;-0.0727990078685463164520295;-0.0719081292873489583205123;-0.0710173648334547409336790;
+     -0.0701267130644115407438380;-0.0692361725390296563667783;-0.0683457418173650832463727;-0.0674554194607028247723011;
+     -0.0665652040315402393386739;-0.0656750940935704228333480;-0.0647850882116656260497053;-0.0638951849518607065146125;
+     -0.0630053828813366142281973;-0.0621156805684039108129696;-0.0612260765824863215716701;-0.0603365694941043199550639;
+     -0.0594471578748587439426946;-0.0585578402974144438413873;-0.0576686153354839610080321;-0.0567794815638112370048940;
+     -0.0558904375581553526973823;-0.0550014818952742968058675;-0.0541126131529087634247656;-0.0532238299097659780237109;
+     -0.0523351307455035514472096;-0.0514465142407133614307148;-0.0505579789769054611525802;-0.0496695235364920143428400;
+     -0.0487811465027712564712261;-0.0478928464599114815382702;-0.0470046219929350539947460;-0.0461164716877024453160887;
+     -0.0452283941308962947597849;-0.0443403879100054938350531;-0.0434524516133092940154370;-0.0425645838298614372262090;
+     -0.0416767831494743086397292;-0.0407890481627031113131293;-0.0399013774608300622038846;-0.0390137696358486091000096;
+     -0.0381262232804476680027560;-0.0372387369879958805008091;-0.0363513093525258906760746;-0.0354639389687186410822098;
+     -0.0345766244318876873380993;-0.0336893643379635308794875;-0.0328021572834779694129740;-0.0319150018655484646175391;
+     -0.0310278966818625266397110;-0.0301408403306621149293976;-0.0292538314107280549642983;-0.0283668685213644704116734;
+     -0.0274799502623832302770910;-0.0265930752340884105905852;-0.0257062420372607701814493;-0.0248194492731422400936554;
+     -0.0239326955434204261946306;-0.0230459794502131245308374;-0.0221592995960528489842988;-0.0212726545838713707848753;
+     -0.0203860430169842694337440;-0.0194994634990754945941500;-0.0186129146341819385060947;-0.0177263950266780184821951;
+     -0.0168399032812602690424965;-0.0159534380029319432465412;-0.0150669977969876227814964;-0.0141805812689978363656155;
+     -0.0132941870247936860267611;-0.0124078136704514808161406;-0.0115214598122773775178119;-0.0106351240567920279148892;
+     -0.00974880501071523217374178;-0.00886250128095059790780139;-0.00797621147457020448290716;-0.00708993419879927212639732;
+     -0.00620366806100083540241769;-0.00531741166866042061615295;-0.00443116362937072670989908;-0.00354492255081630921408397;
+     -0.00265868704075826681650862;-0.00177245570701893011322260;-0.000886227157466552104565435;0.;
+     0.000886227157466552104565435;0.00177245570701893011322260;0.00265868704075826681650862;0.00354492255081630921408397;
+     0.00443116362937072670989908;0.00531741166866042061615295;0.00620366806100083540241769;0.00708993419879927212639732;
+     0.00797621147457020448290716;0.00886250128095059790780139;0.00974880501071523217374178;0.0106351240567920279148892;
+     0.0115214598122773775178119;0.0124078136704514808161406;0.0132941870247936860267611;0.0141805812689978363656155;
+     0.0150669977969876227814964;0.0159534380029319432465412;0.0168399032812602690424965;0.0177263950266780184821951;
+     0.0186129146341819385060947;0.0194994634990754945941500;0.0203860430169842694337440;0.0212726545838713707848753;
+     0.0221592995960528489842988;0.0230459794502131245308374;0.0239326955434204261946306;0.0248194492731422400936554;
+     0.0257062420372607701814493;0.0265930752340884105905852;0.0274799502623832302770910;0.0283668685213644704116734;
+     0.0292538314107280549642983;0.0301408403306621149293976;0.0310278966818625266397110;0.0319150018655484646175391;
+     0.0328021572834779694129740;0.0336893643379635308794875;0.0345766244318876873380993;0.0354639389687186410822098;
+     0.0363513093525258906760746;0.0372387369879958805008091;0.0381262232804476680027560;0.0390137696358486091000096;
+     0.0399013774608300622038846;0.0407890481627031113131293;0.0416767831494743086397292;0.0425645838298614372262090;
+     0.0434524516133092940154370;0.0443403879100054938350531;0.0452283941308962947597849;0.0461164716877024453160887;
+     0.0470046219929350539947460;0.0478928464599114815382702;0.0487811465027712564712261;0.0496695235364920143428400;
+     0.0505579789769054611525802;0.0514465142407133614307148;0.0523351307455035514472096;0.0532238299097659780237109;
+     0.0541126131529087634247656;0.0550014818952742968058675;0.0558904375581553526973823;0.0567794815638112370048940;
+     0.0576686153354839610080321;0.0585578402974144438413873;0.0594471578748587439426946;0.0603365694941043199550639;
+     0.0612260765824863215716701;0.0621156805684039108129696;0.0630053828813366142281973;0.0638951849518607065146125;
+     0.0647850882116656260497053;0.0656750940935704228333480;0.0665652040315402393386739;0.0674554194607028247723011;
+     0.0683457418173650832463727;0.0692361725390296563667783;0.0701267130644115407438380;0.0710173648334547409336790;
+     0.0719081292873489583205123;0.0727990078685463164520295;0.0736900020207781233421723;0.0745811131890716712576034;
+     0.0754723428197670745063028;0.0763636923605341457488485;0.0772551632603893113550998;0.0781467569697125663311978;
+     0.0790384749402644693440208;0.0799303186252031783724914;0.0808222894791015275174189;0.0817143889579641455038835;
+     0.0826066185192446164125211;0.0834989796218626831784544;0.0843914737262214943990342;0.0852841022942248949940076;
+     0.0861768667892947612642142;0.0870697686763883808974295;0.0879628094220158784725300;0.0888559904942576870157373;
+     0.0897493133627820661653203;0.0906427794988626675037925;0.0915363903753961476193256;0.0924301474669198294608308;
+     0.0933240522496294125539148;0.0942181062013967326477174;0.0951123108017875713654626;0.0960066675320795164344287;
+     0.0969011778752798730739420;0.0977958433161436271229397;0.0986906653411914604916251;0.0995856454387278195247507;
+     0.100480785098859036867117;0.101376085813511507424966;0.102271549076449919020062;0.103167176383295538336451;
+     0.104062969231544552763038;0.104958929120586468738427;0.105855057551722567207705;0.106751356028184416804184;
+     0.107647826055152445372494;0.108544469139774570452789;0.109441286791184889349302;0.110338280520522429409946;
+     0.111235451840949959147184;0.112132802267672860833969;0.113030333317958065212148;0.113928046511153048954376;
+     0.114825943368704895524282;0.115724025414179420083364;0.116622294173280359096853;0.117520751173868625294633;
+     0.118419397945981628647141;0.119318236021852664020104;0.120217266935930366175908;0.121116492224898232793403;
+     0.122015913427694216181984;0.122915532085530384369896;0.123815349741912652250820;0.124715367942660583477010;
+     0.125615588235927263791463;0.126516012172219246495889;0.127416641304416570755579;0.128317477187792853446625;
+     0.129218521380035455255404;0.130119775441265721744671;0.131021240934059300105162;0.131922919423466532316157;
+     0.132824812477032925443096;0.133726921664819699805000;0.134629248559424415749181;0.135531794736001679775511;
+     0.136434561772283930757331;0.137337551248602307010996;0.138240764747907594970952;0.139144203855791260232257;
+     0.140047870160506561727510;0.140951765252989749810222;0.141855890726881349021849;0.142760248178547526324913;
+     0.143664839207101545589882;0.144569665414425309128859;0.145474728405190987074448;0.146380029786882735407659;
+     0.147285571169818503444195;0.148191354167171931594021;0.149097380394994340214742;0.150003651472236810385003;
+     0.150910169020772357429845;0.151816934665418198035783;0.152723950033958111799220;0.153631216757164898057737;
+     0.154538736468822928859817;0.155446510805750798934573;0.156354541407824073529228;0.157262829917998134988229;
+     0.158171377982331128954168;0.159080187250007011076982;0.159989259373358695124303;0.160898596007891303392273;
+     0.161808198812305520322682;0.162718069448521050238859;0.163628209581700180119424;0.164538620880271448335743;
+     0.165449305015953420285730;0.166360263663778571863528;0.167271498502117281711534;0.168183011212701933208286;
+     0.169094803480651127152813;0.170006876994494006113221;0.170919233446194691414565;0.171831874531176833748353;
+     0.172744801948348278393465;0.173658017400125846045738;0.174571522592460230261037;0.175485319234861012524278;
+     0.176399409040421795964603;0.177313793725845458744699;0.178228475011469528160158;0.179143454621291676492749;
+     0.180058734282995339669514;0.180974315727975459787782;0.181890200691364352574378;0.182806390912057700855673;
+     0.183722888132740675123476;0.184639694099914182290314;0.185556810563921243736198;0.186474239278973503757652;
+     0.187391982003177869538579;0.188310040498563283771365;0.189228416531107631065596;0.190147111870764779290812;
+     0.191066128291491757008872;0.191985467571276068160731;0.192905131492163145181806;0.193825121840283941729502;
+     0.194745440405882666216067;0.195666088983344657349530;0.196587069371224402895272;0.197508383372273702880600;
+     0.198430032793469978474653;0.199352019446044727786039;0.200274345145512129830755;0.201197011711697797933239;
+     0.202120020968767683833755;0.203043374745257133785850;0.203967074874100097938185;0.204891123192658494305804;
+     0.205815521542751728646714;0.206740271770686371570609;0.207665375727285994217637;0.208590835267921163856305;
+     0.209516652252539600760912;0.210442828545696497740329;0.211369366016585003701502;0.212296266539066872642729;
+     0.213223531991703279483542;0.214151164257785804149977;0.215079165225367585346059;0.216007536787294645454492;
+     0.216936280841237388021901;0.217865399289722269296372;0.218794894040163645297663;0.219724767004895795913128;
+     0.220655020101205127525302;0.221585655251362555690039;0.222516674382656069397254;0.223448079427423478459602;
+     0.224379872323085345587816;0.225312055012178104725014;0.226244629442387367225997;0.227177597566581417481393;
+     0.228110961342844899600545;0.229044722734512696781181;0.229978883710204005008221;0.230913446243856602738561;
+     0.231848412314761318243294;0.232783783907596696293605;0.233719563012463865891554;0.234655751624921610762049;
+     0.235592351746021644337605;0.236529365382344090982915;0.237466794546033175221892;0.238404641254833120745613;
+     0.239342907532124260995555;0.240281595406959363132655;0.241220706914100167219025;0.242160244094054142455665;
+     0.243100208993111462336147;0.244040603663382200593156;0.244981430162833749831766;0.245922690555328464760578;
+     0.246864386910661531949278;0.247806521304599068058765;0.248749095818916448507819;0.249692112541436868558283;
+     0.250635573566070138818935;0.251579480992851717186640;0.252523836927981979261970;0.253468643483865729295322;
+     0.254413902779151953738566;0.255359616938773819496520;0.256305788093988918991992;0.257252418382419764177804;
+     0.258199509948094531649116;0.259147064941488061029464;0.260095085519563108824305;0.261043573845811859956383;
+     0.261992532090297699218086;0.262941962429697244896955;0.263891867047342646851815;0.264842248133264151338477;
+     0.265793107884232934905753;0.266744448503804209704503;0.267696272202360602574682;0.268648581197155810297895;
+     0.269601377712358533425668;0.270554663979096691116723;0.271508442235501919439778;0.272462714726754355621958;
+     0.273417483705127710746720;0.274372751430034633429268;0.275328520168072367021790;0.276284792193068702925524;
+     0.277241569786128232611518;0.278198855235678900977229;0.279156650837518863691535;0.280114958894863651206581;
+     0.281073781718393642140920;0.282033121626301848764835;0.282992980944342017345393;0.283953362005877046135800;
+     0.284914267151927723820923;0.285875698731221791258505;0.286837659100243329383496;0.287800150623282476171267;
+     0.288763175672485475584010;0.289726736627905061453614;0.290690835877551179283553;0.291655475817442048981958;
+     0.292620658851655571567981;0.293586387392381082923904;0.294552663859971457696095;0.295519490682995566478939;
+     0.296486870298291089447284;0.297454805151017689634687;0.298423297694710549086888;0.299392350391334271152468;
+     0.300361965711337152205522;0.301332146133705826128503;0.302302894146020284917042;0.303274212244509278802667;
+     0.304246102934106099323804;0.305218568728504748810367;0.306191612150216499782543;0.307165235730626847800129;
+     0.308139442010052861334930;0.309114233537800932275328;0.310089612872224930709180;0.311065582580784767668659;
+     0.312042145240105369558602;0.313019303436036068028315;0.313997059763710409085611;0.314975416827606385291216;
+     0.315954377241607094911430;0.316933943629061831947228;0.317914118622847610998746;0.318894904865431130965352;
+     0.319876305008931181623254;0.320858321715181497164882;0.321840957655794060827042;0.322824215512222864778164;
+     0.323808097975828129478784;0.324792607747940986773783;0.325777747539928631019810;0.326763520073259942596795;
+     0.327749928079571588198463;0.328736974300734602343382;0.329724661488921454595204;0.330712992406673607028554;
+     0.331701969826969566525310;0.332691596533293436535013;0.333681875319703972982629;0.334672808990904149057107;
+     0.335664400362311233664909;0.336656652260127388384144;0.337649567521410787806985;0.338643148994147268210747;
+     0.339637399537322509551397;0.340632322020994755827280;0.341627919326368078915602;0.342624194345866191039560;
+     0.343621149983206811080181;0.344618789153476590003681;0.345617114783206600732730;0.346616129810448397848247;
+     0.347615837184850652567335;0.348616239867736368502744;0.349617340832180683769720;0.350619143063089265067382;
+     0.351621649557277299423852;0.352624863323549089357157;0.353628787382778257267641;0.354633424767988564942025;
+     0.355638778524435354114617;0.356644851709687614097242;0.357651647393710682556510;0.358659168658949585584833;
+     0.359667418600413023280381;0.360676400325758007120738;0.361686116955375155485562;0.362696571622474653754979;
+     0.363707767473172885482797;0.364719707666579741216936;0.365732395374886611613717;0.366745833783455071567895;
+     0.367760026090906262156518;0.368774975509210977271926;0.369790685263780461897400;0.370807158593557929058249;
+     0.371824398751110802561401;0.372842409002723692717928;0.373861192628492112325372;0.374880752922416940270228;
+     0.375901093192499640195609;0.376922216760838241764837;0.377944126963724092138613;0.378966827151739385371447;
+     0.379990320689855477522284;0.381014610957531995364634;0.382039701348816746673194;0.383065595272446440156734;
+     0.384092296151948223201199;0.385119807425742045682276;0.386148132547243858203376;0.387177274984969653212906;
+     0.388207238222640357554016;0.389238025759287585100599;0.390269641109360258235345;0.391302087802832107029012;
+     0.392335369385310055084850;0.393369489418143501118366;0.394404451478534505450229;0.395440259159648890699302;
+     0.396476916070728266073398;0.397514425837202984767494;0.398552792100806044092875;0.399592018519687938075904;
+     0.400632108768532472381979;0.401673066538673551538719;0.402714895538212948552500;0.403757599492139067134282;
+     0.404801182142446706874082;0.405845647248257841828688;0.406890998585943423114100;0.407937239949246216222914;
+     0.408984375149404683917372;0.410032408015277925681117;0.411081342393471684846923;0.412131182148465434653718;
+     0.413181931162740554624251;0.414233593336909608794685;0.415286172589846737469336;0.416339672858819174317723;
+     0.417394098099619900777074;0.418449452286701449871496;0.419505739413310871709193;0.420562963491625873071421;
+     0.421621128552892143661369;0.422680238647561881737857;0.423740297845433532017696;0.424801310235792748891791;
+     0.425863279927554598163639;0.426926211049407010684775;0.427990107749955501430070;0.429054974197869167726505;
+     0.430120814582027980522297;0.431187633111671382759028;0.432255434016548209087678;0.433324221547067941350446;
+     0.434393999974453314433743;0.435464773590894287284015;0.436536546709703394067020;0.437609323665472490642935;
+     0.438683108814230911724245;0.439757906533605054280818;0.440833721222979402956917;0.441910557303659013468255;
+     0.442988419219033470153526;0.444067311434742334064264;0.445147238438842098189415;0.446228204741974666626704;
+     0.447310214877537374731816;0.448393273401854567498595;0.449477384894350753649022;0.450562553957725353140641;
+     0.451648785218129056031497;0.452736083325341810878516;0.453824452952952461084706;0.454913898798540047853642;
+     0.456004425583856798656448;0.457096038055012820367009;0.458188740982662516475480;0.459282539162192748048363;
+     0.460377437413912758365590;0.461473440583245881431235;0.462570553540923054824727;0.463668781183178157633870;
+     0.464768128431945194489628;0.465868600235057347005579;0.466970201566447914212286;0.468072937426353163868614;
+     0.469176812841517116828335;0.470281832865398286941320;0.471388002578378399274204;0.472495327087973109745849;
+     0.473603811529044749588163;0.474713461064017118363061;0.475824280883092349591560;0.476936276204469873381418;
+     0.478049452274567500775245;0.479163814368244654881929;0.480279367789027774200510;0.481396117869337913897394;
+     0.482514069970720571155200;0.483633229484077761074596;0.484753601829902369979366;0.485875192458514813349731;
+     0.486998006850302025989727;0.488122050515958812421373;0.489247328996731585891494;0.490373847864664524776547;
+     0.491501612722848175576772;0.492630629205670532103474;0.493760902979070620882506;0.494892439740794623223041;
+     0.496025245220654564833729;0.497159325180789604308388;0.498294685415929952250655;0.499431331753663453261643;
+     0.500569270054704863476697;0.501708506213167856807092;0.502849046156839793519900;0.503990895847459285274674;
+     0.505134061280996591228925;0.506278548487936880326002;0.507424363533566395388924;0.508571512518261555162144;
+     0.509720001577781030970368;0.510869836883560835199485;0.512021024643012459349639;0.513173571099824099964556;
+     0.514327482534265011304754;0.515482765263493024205197;0.516639425641865271140692;0.517797470061252158114871;
+     0.518956904951354624591287;0.520117736780024733298055;0.521279972053589632360901;0.522443617317178932853529;
+     0.523608679155055545499181;0.524775164190950020913310;0.525943079088398438444614;0.527112430551083889350577;
+     0.528283225323181600734249;0.529455470189707747371634;0.530629171976871999273835;0.531804337552433853555395;
+     0.532980973826062799920226;0.534159087749702369829439;0.535338686317938120181512;0.536519776568369603114798;
+     0.537702365581986374335740;0.538886460483548093183440;0.540072068441968768462873;0.541259196670705204915226;
+     0.542447852428149706044881;0.543638043018027089888778;0.544829775789796075195604;0.546023058139055096379661;
+     0.547217897507952606527904;0.548414301385601928668550;0.549612277308500716456468;0.550811832860955086394373;
+     0.552012975675508484690151;0.553215713433375352849737;0.554420053864879657122263;0.555626004749898347949980;
+     0.556833573918309816630205;0.558042769250447417470635;0.559253598677558124813090;0.560466070182266395414699;
+     0.561680191799043307809948;0.562895971614681051432447;0.564113417768772839452082;0.565332538454198320481909;
+     0.566553341917614565530137;0.567775836459952707816332;0.569000030436920314338033;0.570225932259509569364741;
+     0.571453550394511351351370;0.572682893365035286103023;0.573913969751035860388153;0.575146788189844681588133;
+     0.576381357376708970388645;0.577617686065336374962660;0.578855783068446196566639;0.580095657258327117971639;
+     0.581337317567401527679750;0.582580772988796534434417;0.583826032576921768121352;0.585073105448054064775488;
+     0.586322000780929135059570;0.587572727817340317262067;0.588825295862744517576917;0.590079714286875442175869;
+     0.591335992524364227366604;0.592594140075367575947128;0.593854166506203509719984;0.595116081449994850019300;
+     0.596379894607320540030513;0.597645615746874924647529;0.598913254706135105615973;0.600182821392036491754959;
+     0.601454325781656666134284;0.602727777922907694210153;0.604003187935236999091285;0.605280566010336932319637;
+     0.606559922412863170806831;0.607841267481162072869936;0.609124611628007128659254;0.610409965341344642667636;
+     0.611697339185048788456327;0.612986743799686178227838;0.614278189903290092422814;0.615571688292144517116603;
+     0.616867249841578139643369;0.618164885506768455582459;0.619464606323556143004491;0.620766423409269862694708;
+     0.622070347963561645949760;0.623376391269253034482742;0.624684564693192139971279;0.625994879687121793846348;
+     0.627307347788558961046672;0.628621980621685594656557;0.629938789898251111605494;0.631257787418486672937314;
+     0.632578985072031455556846;0.633902394838871105834534;0.635228028790288568996146;0.636555899089827491847212;
+     0.637886017994268400082143;0.639218397854617855207882;0.640553051117110799973398;0.641889990324226305141360;
+     0.643229228115716934468944;0.644570777229651948883019;0.645914650503474575043100;0.647260860875073567785696;
+     0.648609421383869300338197;0.649960345171914620681705;0.651313645485010717032493;0.652669335673838240103679;
+     0.654027429195103934604611;0.655387939612703037338166;0.656750880598897704268154;0.658116265935511734053255;
+     0.659484109515141860783048;0.660854425342385894008748;0.662227227535087989639255;0.663602530325601340875000;
+     0.664980348062068584081204;0.666360695209720220361604;0.667743586352191359587010;0.669129036192857099763483;
+     0.670517059556186860896197;0.671907671389117998920711;0.673300886762449031837214;0.674696720872252816899203;
+     0.676095189041310024580015;0.677496306720563262072639;0.678900089490592206274711;0.680306553063110113575629;
+     0.681715713282482081301025;0.683127586127265443385831;0.684542187711772690745695;0.685959534287657314902413;
+     0.687379642245522981697320;0.688802528116556450402505;0.690228208574184662218446;0.691656700435756431033936;
+     0.693088020664249178425344;0.694522186370001164193449;0.695959214812469673283114;0.697399123402015629710152;
+     0.698841929701715118137356;0.700287651429198304004255;0.701736306458516253629577;0.703187912822036166478639;
+     0.704642488712365542827071;0.706100052484305821364947;0.707560622656836032879092;0.709024217915127028034019;
+     0.710490857112586849451786;0.711960559272937830776354;0.713433343592326018207544;0.714909229441463523112322;
+     0.716388236367804427776112;0.717870384097754880153633;0.719355692538918027627315;0.720844181782374454291622;
+     0.722335872104998801164239;0.723830783971813263990758;0.725328938038378678969375;0.726830355153223922787785;
+     0.728335056360314369847696;0.729843062901560166465567;0.731354396219365099193895;0.732869077959216852218817;
+     0.734387129972319467070368;0.735908574318268836645561;0.737433433267772084805852;0.738961729305411702584492;
+     0.740493485132455332341197;0.742028723669712112047386;0.743567468060436513291432;0.745109741673280628576998;
+     0.746655568105295886066185;0.748204971184985193111121;0.749757974975406533741602;0.751314603777329069751936;
+     0.752874882132442820177428;0.754438834826623019790749;0.756006486893250283802416;0.757577863616587733240031;
+     0.759152990535216263530909;0.760731893445529167646167;0.762314598405287354806035;0.763901131737236436221627;
+     0.765491520032786980684311;0.767085790155759275037507;0.768683969246193957705777;0.770286084724229927541830;
+     0.771892164294050965314105;0.773502235947902541227630;0.775116327970180318981467;0.776734468941591905051428;
+     0.778356687743393431181906;0.779983013561702598512212;0.781613475891889853388483;0.783248104543049407761361;
+     0.784886929642551861182764;0.786529981640680226834519;0.788177291315351210791011;0.789828889776923641882848;
+     0.791484808473095999135830;0.793145079193895034858203;0.794809734076757544089917;0.796478805611707385363143;
+     0.798152326646629913608153;0.799830330392646043629732;0.801512850429588221935367;0.803199920711580645878829;
+     0.804891575572726132154983;0.806587849732902101709751;0.808288778303668215181754;0.809994396794288262140596;
+     0.811704741117868978705118;0.813419847597618541690289;0.815139752973227563323866;0.816864494407375488876719;
+     0.818594109492365380350606;0.820328636256890152754014;0.822068113172933415564498;0.823812579162808160822003;
+     0.825562073606336631023609;0.827316636348174794701583;0.829076307705284955373421;0.830841128474560120569394;
+     0.832611139940603861988961;0.834386383883669505636634;0.836166902587762602169693;0.837952738848910742789354;
+     0.839743935983604904963956;0.841540537837416635233851;0.843342588793795503465446;0.845150133783051394355273;
+     0.846963218291526337899961;0.848781888370960721117400;0.850606190648058868708733;0.852436172334259130778193;
+     0.854271881235713771374789;0.856113365763484112691231;0.857960674943956556465451;0.859813858429485276701938;
+     0.861672966509267556497199;0.863538050120457926759686;0.865409160859527456214132;0.867286350993874740539321;
+     0.869169673473695344084982;0.871059181944116660638341;0.872954930757605380467648;0.874856974986654979676502;
+     0.876765370436760885091340;0.878680173659691213822623;0.880601441967061241651857;0.882529233444220018882342;
+     0.884463606964457826649897;0.886404622203543451338165;0.888352339654600550118256;0.890306820643332687192334;
+     0.892268127343606938544748;0.894236322793406293395002;0.896211470911161423631094;0.898193636512472748831602;
+     0.900182885327234094639705;0.902179284017169626839285;0.904182900193796143139771;0.906193802436823220071163;
+     0.908212060313004144225512;0.910237744395451006093112;0.912270926283427801704433;0.914311678622635873015377;
+     0.916360075126006523317668;0.918416190595016169815512;0.920480100941539942827582;0.922551883210260210841589;
+     0.924631615601647103910482;0.926719377495528725731954;0.928815249475269388351157;0.930919313352574873985593;
+     0.933031652192944427270628;0.935152350341789909611329;0.937281493451243306720676;0.939419168507674572321134;
+     0.941565463859942616962252;0.943720469248403112623036;0.945884275834697682976616;0.948056976232349987744631;
+     0.950238664538195189410099;0.952429436364670313748928;0.954629388872994084354930;0.956838620807265927866568;
+     0.959057232529515013380354;0.961285326055731409118287;0.963523005092912714513298;0.965770375077160859351606;
+     0.968027543212865156493316;0.970294618513009154187884;0.972571711840640361501646;0.974858935951543519472376;
+     0.977156405538159765104812;0.979464237274795789251756;0.981782549864168927059146;0.984111464085336045518998;
+     0.986451102843056111574170;0.988801591218638441249043;0.991163056522330850848242;0.993535628347304261113059;
+     0.995919438625292750450167;0.998314621683950620428471;1.00072131430599073256377;1.00313965579017120729917;
+     1.00556978801420055183020;1.00801185549963441130450;1.01046600547884042676895;1.01293238796411114243671;
+     1.01541115581900854441199;1.01790246483202764360873;1.02040647379267054860440;1.02292334457002672171526;
+     1.02545324219395958661019;1.02799633493900437212031;1.03055279441108704931261;1.03312279563717946315342;
+     1.03570651715801129305399;1.03830414112396531630198;1.04091585339428861413581;1.04354184363975887267130;
+     1.04618230544895181216668;1.04883743643826305191789;1.05150743836584540982523;1.05419251724963077262264;
+     1.05689288348961428416044;1.05960875199458771638643;1.06234034231351854451039;1.06508787877178148051829;
+     1.06785159061246006671034;1.07063171214294743521462;1.07342848288708754562668;1.07624214774311116967833;
+     1.07907295714763465155357;1.08192116724600409166043;1.08478704006928314130699;1.08767084371819912164160;
+     1.09057285255437976350432;1.09349334739923158235067;1.09643261574083083726659;1.09939095194921926524473;
+     1.10236865750051942874064;1.10536604121030867162096;1.10838341947671646039726;1.11142111653373741828454;
+     1.11447946471528177395617;1.11755880473051639137566;1.12065948595108318104233;1.12378186671081768969028;
+     1.12692631461862921361960;1.13009320688524508602561;1.13328293066456607513313;1.13649588341042734231966;
+     1.13973247324961041494949;1.14299311937200641942412;1.14627825243888971549360;1.14958831501032442296701;
+     1.15292376199279451989056;1.15628506110822163735251;1.15967269338561384167854;1.16308715367667408672625;
+     1.16652895119678919456914;1.16999861009291979542268;1.17349667004001930287996;1.17702368686772646342933;
+     1.18058023321920212215387;1.18416689924411749800915;1.18778429332794946696112;1.19143304285989922198847;
+     1.19511379504192544842046;1.19882721774157318917829;1.20257400039148639653512;1.20635485493871746442485;
+     1.21017051684719268883733;1.21402174615696070959820;1.21790932860414388392794;1.22183407680583284397770;
+     1.22579683151451512054446;1.22979846294701294800680;1.23383987219332687396791;1.23792199271124470601816;
+     1.24204579191308427645764;1.24621227285149871837430;1.25042247601189031583104;1.25467748121966017307579;
+     1.25897840967127347229233;1.26332642609895248768275;1.26772274107973147711542;1.27216861350063008607206;
+     1.27666535319283750241650;1.28121432374906258624426;1.28581694553961189638512;1.29047469894432663626360;
+     1.29518912781926248355837;1.29996184321895768212866;1.30479452739733301983674;1.30968893811273510877283;
+     1.31464691326540954816980;1.31967037589881690640162;1.32476133959973289089009;1.32992191433606380401593;
+     1.33515431277582849530455;1.34046085713588930726969;1.34584398661485169413881;1.35130626547120660708750;
+     1.35685039181539261030877;1.36247920719316111112837;1.36819570704762188366770;1.37400305215884578004754;
+     1.37990458117316776469996;1.38590382434967794527797;1.39200451866918531861457;1.39821062447163920701438;
+     1.40452634381214160594658;1.41095614075394083115303;1.41750476384996789516378;1.42417727110354626495924;
+     1.43097905774508171768046;1.43791588721631257139615;1.44499392581891844119832;1.45221978156224685014342;
+     1.45960054783849765418572;1.46714385266652228379940;1.47485791438202173394033;1.48275160481818015318184;
+     1.49083452122409072994342;1.49911706841829737690302;1.50761055298380584861693;1.51632729169521713246473;
+     1.52528073684951304945905;1.53448562177771813914331;1.54395813058290834226583;1.55371609713170636373557;
+     1.56377923959071921153068;1.57416943844147646813481;1.58491106805948089387746;1.59603139479067878351924;
+     1.60756105826610791904552;1.61953465784330182137717;1.63199147310755983780432;1.64497635713318705017720;
+     1.65854085494046514649335;1.67274461920355980864179;1.68765722376985222676591;1.70336051775180604164031;
+     1.71995172674121274921475;1.73754760637571012247642;1.75629011016036985182043;1.77635428959895142455381;
+     1.79795957771209517163847;1.82138636771844967304021;1.84700119608962619338729;1.87529654534774725251058;
+     1.90695686467094561133550;1.94297485086632095708177;1.98487261261553223101737;2.03516768306608309650651;
+     2.09850761179367089535179;2.18512421913300426570596;2.32675376551352467056022]';
+endfunction
+[Xr,Yr]=erfinv_ref();
+y=erfinv(Xr);
+k=find(abs(Yr)>0);
+if norm((Yr(k)-y(k))./Yr(k),1) >1d-12 then bugmes();quit;end