trunk/noumerov/noumerov.plt

   1 #!/usr/bin/gnuplot
   2 set grid
   3 set y2tics
   4 set x2tics
   5 set isosample 40
   6 l=1;
   7 A =-2.36464e+31
   8 delta = 3.76155
   9 hbar=6.626068E-34
  10 E=5.0
  11 m=hbar*hbar/2;
  12 A=2.7354695938922176e+31
  13 delta=-(4.9303806576313238e-32+0.7+0.3+0.2+0.2)+(l*pi/2)
  14 b=2.236068e+00
  15
  16 # besj0 = \e$B1_Cl%Y%C%;%k!#1_%Y%C%;%k$G$J$$!#\e(B
  17 j0(x) = sin(x)/x;
  18 j1(x) = (sin(x)-cos(x))/x;
  19 j1_(x) = (sin(x)-cos(x));
  20 y0(x) = -cos(x)/x;
  21 y1(x) = -(cos(x)+sin(x))/x;
  22 y1_(x) = -(cos(x)+sin(x));
  23
  24 plot    "tmp/gsl_test.stdout"   using 1:2,\
  25         "tmp/gsl_test.stdout"   using 1:3,\
  26         "tmp/gsl_test.stdout"   using 1:(j0($1)),\
  27         "tmp/gsl_test.stdout"   using 1:(sin($1)/$1)
  28 pause -1
  29 #b = sqrt(2*m*E)/hbar;
  30 #A               = 2.73547e+31      +/- 9.502e+27    (0.03474%)
  31 #delta           = 0.138346         +/- 0.0004599    (0.3324%)
  32 #b               = 2.22822          +/- 4.979e-05    (0.002234%)
  33 A               = 2.73547e+31 - 0.6E+31
  34 delta           = 0.138346
  35 b               = 2.22822
  36
  37 f(x,A,delta,b)=A*sin(b*x-(l*pi/2)+delta);
  38 g(x,A,delta,b)=A*b*x*(cos(delta)*j1(b*x)-sin(delta)*y1(b*x));
  39 g(x,A,delta,b)=A*(cos(delta)*j1_(b*x)-sin(delta)*y1_(b*x));
  40 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0:3000:0 via b
  41 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta,b
  42 #fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A
  43
  44 #A               = 2.13547e+31      +/- 4.767e+28    (0.2232%)
  45 #delta           = -0.79802         +/- 0.002965     (0.3716%)
  46 #b               = 2.24525          +/- 0.0003203    (0.01427%)
  47 A               =  1.93924e+31#2.13547e+31
  48 delta           = -0.79802
  49 b               = 2.24525
  50
  51 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
  52 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A
  53 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via delta
  54 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via b
  55 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
  56 print A;
  57 print delta;
  58 print b;
  59 plot    "tmp/noumerov.stdout"   using 1:2 every 2:::0::0,\
  60         ""                      using 1:2 every 2::4001:0::0,\
  61         f(x,A,delta,b),\
  62         g(x,A,delta,b)
  63 #       ""                      using 1:(f($1,A,delta,b)),\
  64 #       ""                      using 1:(g($1,A,delta,b))
  65 pause -1
  66 exit;
  67 plot    "tmp/noumerov.stdout"   using 0:1 every :::0::0,\
  68         ""                      using 0:2 every :::0::0 axis x1y2,\
  69         ""                      using 0:3 every :::0::0 axis x1y2
  70 pause -1
  71 plot    "tmp/noumerov.stdout"   using 1:4 every :::0::0,\
  72         ""                      using 1:2 every :::0::0 axis x1y2,\
  73         ""                      using 1:3 every :::0::0 axis x1y2
  74 pause -1
  75 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  76         "" using 0:7 every :::2::2 axis x1y2
  77 pause -1
  78 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  79         "" using 0:2 every :::2::2 axis x1y2
  80 pause -1
  81 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  82         "" using 0:4 every :::2::2 axis x1y2
  83 pause -1
  84 exit;
  85 plot [0:5E-16] "tmp/noumerov.stdout" using 1:2 every :::0::0,\
  86         "" using 1:3 every :::0::0,\
  87         "" using 1:4  every :::0::0 axis x1y2,\
  88         "" using 1:5  every :::0::0 axis x1y2
  89 pause -1
  90 plot "tmp/noumerov.stdout" using 1:2 every :::1::1,\
  91         "" using 1:($2+$4) every :::1::1,\
  92         "" using 1:($2+$5) every :::1::1,\
  93         "" using 1:($2+$6) every :::1::1
  94 pause -1
  95 plot "tmp/noumerov.stdout",\
  96         x*exp(-x**2/2)
  97 pause -1
  98 plot [0:1E-14] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
  99         "" using 1:3 every :::2::2 axis x1y2,\
 100         "" using 1:4 every 2::1:2::2
 101 pause -1
 102 x0 = 1E-11;
 103 plot [x0+0:x0+1E-10] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
 104         "" using 1:4 every :::2::2 axis x1y2,\
 105         "" using 1:5 every 2::1:2::2
 106 pause -1