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:4,\
  27         "tmp/gsl_test.stdout"   using 1:(j1($1))
  28 pause -1
  29 exit;
  30 #b = sqrt(2*m*E)/hbar;
  31 #A               = 2.73547e+31      +/- 9.502e+27    (0.03474%)
  32 #delta           = 0.138346         +/- 0.0004599    (0.3324%)
  33 #b               = 2.22822          +/- 4.979e-05    (0.002234%)
  34 A               = 2.73547e+31 - 0.6E+31
  35 delta           = 0.138346
  36 b               = 2.22822
  37
  38 f(x,A,delta,b)=A*sin(b*x-(l*pi/2)+delta);
  39 g(x,A,delta,b)=A*b*x*(cos(delta)*j1(b*x)-sin(delta)*y1(b*x));
  40 g(x,A,delta,b)=A*(cos(delta)*j1_(b*x)-sin(delta)*y1_(b*x));
  41 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0:3000:0 via b
  42 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta,b
  43 #fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A
  44
  45 #A               = 2.13547e+31      +/- 4.767e+28    (0.2232%)
  46 #delta           = -0.79802         +/- 0.002965     (0.3716%)
  47 #b               = 2.24525          +/- 0.0003203    (0.01427%)
  48 A               =  1.93924e+31#2.13547e+31
  49 delta           = -0.79802
  50 b               = 2.24525
  51
  52 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
  53 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A
  54 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via delta
  55 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via b
  56 fit g(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::4000:0::0 via A,delta,b
  57 print A;
  58 print delta;
  59 print b;
  60 plot    "tmp/noumerov.stdout"   using 1:2 every 2:::0::0,\
  61         ""                      using 1:2 every 2::4001:0::0,\
  62         f(x,A,delta,b),\
  63         g(x,A,delta,b)
  64 #       ""                      using 1:(f($1,A,delta,b)),\
  65 #       ""                      using 1:(g($1,A,delta,b))
  66 pause -1
  67 exit;
  68 plot    "tmp/noumerov.stdout"   using 0:1 every :::0::0,\
  69         ""                      using 0:2 every :::0::0 axis x1y2,\
  70         ""                      using 0:3 every :::0::0 axis x1y2
  71 pause -1
  72 plot    "tmp/noumerov.stdout"   using 1:4 every :::0::0,\
  73         ""                      using 1:2 every :::0::0 axis x1y2,\
  74         ""                      using 1:3 every :::0::0 axis x1y2
  75 pause -1
  76 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  77         "" using 0:7 every :::2::2 axis x1y2
  78 pause -1
  79 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  80         "" using 0:2 every :::2::2 axis x1y2
  81 pause -1
  82 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  83         "" using 0:4 every :::2::2 axis x1y2
  84 pause -1
  85 exit;
  86 plot [0:5E-16] "tmp/noumerov.stdout" using 1:2 every :::0::0,\
  87         "" using 1:3 every :::0::0,\
  88         "" using 1:4  every :::0::0 axis x1y2,\
  89         "" using 1:5  every :::0::0 axis x1y2
  90 pause -1
  91 plot "tmp/noumerov.stdout" using 1:2 every :::1::1,\
  92         "" using 1:($2+$4) every :::1::1,\
  93         "" using 1:($2+$5) every :::1::1,\
  94         "" using 1:($2+$6) every :::1::1
  95 pause -1
  96 plot "tmp/noumerov.stdout",\
  97         x*exp(-x**2/2)
  98 pause -1
  99 plot [0:1E-14] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
 100         "" using 1:3 every :::2::2 axis x1y2,\
 101         "" using 1:4 every 2::1:2::2
 102 pause -1
 103 x0 = 1E-11;
 104 plot [x0+0:x0+1E-10] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
 105         "" using 1:4 every :::2::2 axis x1y2,\
 106         "" using 1:5 every 2::1:2::2
 107 pause -1