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
  17 plot    "tmp/gsl_test.stdout"   using 1:2,\
  18         "tmp/gsl_test.stdout"   using 1:(besj0($1)),\
  19         "tmp/gsl_test.stdout"   using 1:(sin($1)/$1)
  20 pause -1
  21 exit
  22 #b = sqrt(2*m*E)/hbar;
  23 f(x,A,delta,b)=A*sin(b*x-(l*pi/2)+delta);
  24 g(x,A,delta,b)=A*b*x*(cos(delta)*besj1(b*x)-sin(delta)*besy1(b*x));
  25 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0:3000:0 via b
  26 fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta,b
  27 #fit g(x,A,delta,1) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta
  28 plot    "tmp/noumerov.stdout"   using 1:2 every 2:::0::0,\
  29         ""                      using 1:2 every 2::2001:0::0,\
  30         f(x,A,delta,b),\
  31         g(x,A,delta,b) axis x1y2
  32 pause -1
  33 exit;
  34 plot    "tmp/noumerov.stdout"   using 0:1 every :::0::0,\
  35         ""                      using 0:2 every :::0::0 axis x1y2,\
  36         ""                      using 0:3 every :::0::0 axis x1y2
  37 pause -1
  38 plot    "tmp/noumerov.stdout"   using 1:4 every :::0::0,\
  39         ""                      using 1:2 every :::0::0 axis x1y2,\
  40         ""                      using 1:3 every :::0::0 axis x1y2
  41 pause -1
  42 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  43         "" using 0:7 every :::2::2 axis x1y2
  44 pause -1
  45 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  46         "" using 0:2 every :::2::2 axis x1y2
  47 pause -1
  48 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  49         "" using 0:4 every :::2::2 axis x1y2
  50 pause -1
  51 exit;
  52 plot [0:5E-16] "tmp/noumerov.stdout" using 1:2 every :::0::0,\
  53         "" using 1:3 every :::0::0,\
  54         "" using 1:4  every :::0::0 axis x1y2,\
  55         "" using 1:5  every :::0::0 axis x1y2
  56 pause -1
  57 plot "tmp/noumerov.stdout" using 1:2 every :::1::1,\
  58         "" using 1:($2+$4) every :::1::1,\
  59         "" using 1:($2+$5) every :::1::1,\
  60         "" using 1:($2+$6) every :::1::1
  61 pause -1
  62 plot "tmp/noumerov.stdout",\
  63         x*exp(-x**2/2)
  64 pause -1
  65 plot [0:1E-14] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
  66         "" using 1:3 every :::2::2 axis x1y2,\
  67         "" using 1:4 every 2::1:2::2
  68 pause -1
  69 x0 = 1E-11;
  70 plot [x0+0:x0+1E-10] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
  71         "" using 1:4 every :::2::2 axis x1y2,\
  72         "" using 1:5 every 2::1:2::2
  73 pause -1