unstable-commit

[applied-gslp/test1.git] / trunk / noumerov / noumerov.plt

diff --git a/trunk/noumerov/noumerov.plt b/trunk/noumerov/noumerov.plt
index 9e8e3b8..ea736a6 100644 (file)
--- a/trunk/noumerov/noumerov.plt
+++ b/trunk/noumerov/noumerov.plt
@@ -13,24 +13,61 @@ A=2.7354695938922176e+31
 delta=-(4.9303806576313238e-32+0.7+0.3+0.2+0.2)+(l*pi/2)
 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)/x-cos(x))/x;
+j1_(x) = (sin(x)/x-cos(x));
+y0(x) = -cos(x)/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:(besj0($1)),\
-       "tmp/gsl_test.stdout"   using 1:(sin($1)/$1)
-pause -1
-exit
 #b = sqrt(2*m*E)/hbar;
+#A               = 2.73547e+31      +/- 9.502e+27    (0.03474%)
+#delta           = 0.138346         +/- 0.0004599    (0.3324%)
+#b               = 2.22822          +/- 4.979e-05    (0.002234%)
+
+#A               = 2.73547e+31 - 0.6E+31
+#delta           = 0.138346
+#b               = 2.22822
+#A               = 2.73547e+31
+A=2.73547e+31
+delta=0.138474818354818
+b=2.22818982634668
+
+#A=1E31
+#delta=0.0
 f(x,A,delta,b)=A*sin(b*x-(l*pi/2)+delta);
-g(x,A,delta,b)=A*b*x*(cos(delta)*besj1(b*x)-sin(delta)*besy1(b*x));
+g(x,A,delta,b)=A*b*x*(cos(delta)*j1(b*x)-sin(delta)*y1(b*x));
+#g(x,A,delta,b)=A*(cos(delta)*j1_(b*x)-sin(delta)*y1_(b*x));
 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0:3000:0 via b
 fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta,b
-#fit g(x,A,delta,1) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A
+
+#A               = 2.13547e+31      +/- 4.767e+28    (0.2232%)
+#delta           = -0.79802         +/- 0.002965     (0.3716%)
+#b               = 2.24525          +/- 0.0003203    (0.01427%)
+#A               =  1.93924e+31#2.13547e+31
+#delta           = -0.79802
+#b               = 2.24525
+
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via delta
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via b
+#fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
+#print A;
+#print delta;
+#print b;
 plot   "tmp/noumerov.stdout"   using 1:2 every 2:::0::0,\
        ""                      using 1:2 every 2::2001:0::0,\
-       f(x,A,delta,b),\
-       g(x,A,delta,b) axis x1y2
+       f(x,A,delta,b)
+pause -1
+fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::6000:0::0 via A,delta,b
+plot   "tmp/noumerov.stdout"   using 1:2 every 2:::0::0,\
+       ""                      using 1:2 every 2::6001:0::0,\
+       ""                      using 1:(g($1,A,delta,b)) every ::6001:0::0
+       
 pause -1
-exit;
 plot   "tmp/noumerov.stdout"   using 0:1 every :::0::0,\
        ""                      using 0:2 every :::0::0 axis x1y2,\
        ""                      using 0:3 every :::0::0 axis x1y2
@@ -71,3 +108,14 @@ plot [x0+0:x0+1E-10] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
        "" using 1:4 every :::2::2 axis x1y2,\
        "" using 1:5 every 2::1:2::2
 pause -1
+plot   "tmp/gsl_test.stdout"   using 1:5,\
+       "tmp/gsl_test.stdout"   using 1:6,\
+       "tmp/gsl_test.stdout"   using 1:7,\
+       "tmp/gsl_test.stdout"   using 1:(y1($1))
+pause -1
+plot   "tmp/gsl_test.stdout"   using 1:2,\
+       "tmp/gsl_test.stdout"   using 1:3,\
+       "tmp/gsl_test.stdout"   using 1:4,\
+       "tmp/gsl_test.stdout"   using 1:(j1($1))
+pause -1
+exit;