unstable-commit

git-svn-id: svn+ssh://svn.sourceforge.jp/svnroot/applied-gslp@176 0e905424-1750-4201-adbd-bf076348bc1c

diff --git a/trunk/noumerov/noumerov.plt b/trunk/noumerov/noumerov.plt
index d1da956..b7db3cc 100644 (file)
--- a/trunk/noumerov/noumerov.plt
+++ b/trunk/noumerov/noumerov.plt
@@ -15,11 +15,11 @@ b=2.236068e+00
 
 # besj0 = \e$B1_Cl%Y%C%;%k!#1_%Y%C%;%k$G$J$$!#\e(B
 j0(x) = sin(x)/x;
-j1(x) = (sin(x)-cos(x))/x;
-j1_(x) = (sin(x)-cos(x));
+j1(x) = (sin(x)/x-cos(x))/x;
+j1_(x) = (sin(x)/x-cos(x));
 y0(x) = -cos(x)/x;
-y1(x) = -(cos(x)+sin(x))/x;
-y1_(x) = -(cos(x)+sin(x));
+y1(x) = -(cos(x)/x+sin(x))/x;
+y1_(x) = -(cos(x)/x+sin(x));
 
 plot   "tmp/gsl_test.stdout"   using 1:2,\
        "tmp/gsl_test.stdout"   using 1:3,\

diff --git a/trunk/test/gsl_test b/trunk/test/gsl_test
index accf854..5d5e865 100755 (executable)
Binary files a/trunk/test/gsl_test and b/trunk/test/gsl_test differ

diff --git a/trunk/test/gsl_test.c b/trunk/test/gsl_test.c
index 438c128..5463374 100644 (file)
--- a/trunk/test/gsl_test.c
+++ b/trunk/test/gsl_test.c
@@ -30,13 +30,15 @@ void test()
        int l = 1;
        double x = 0.0;
        for(i=0;i<1000;i++){
+               x += 0.01;
                printf("% e\t", x);
                printf("% e\t", gsl_sf_bessel_jl(l, x));
                printf("% e\t", gsl_sf_bessel_j1(x));
-               printf("% e\t", (sin(x)-cos(x))/x);
+               printf("% e\t", (sin(x)/x-cos(x))/x);
                printf("% e\t", gsl_sf_bessel_yl(l, x));
+               printf("% e\t", gsl_sf_bessel_y1(x));
+               printf("% e\t", -(cos(x)/x+sin(x))/x);
                printf("% d\n", i);
-               x += 0.1;
        }
        return;
 }

diff --git a/trunk/test/gsl_test.o b/trunk/test/gsl_test.o
index 874a037..6f11832 100644 (file)
Binary files a/trunk/test/gsl_test.o and b/trunk/test/gsl_test.o differ