trunk/noumerov/noumerov.plt

   1 #!/usr/bin/gnuplot
   2 set grid
   3 set y2tics
   4 set x2tics
   5 set isosample 40
   6 A =-2.36464e+31
   7 delta = 3.76155
   8 hbar=6.626068E-34
   9 E=5.0
  10 m=hbar*hbar/2;
  11 A=2.7354695938922176e+31
  12 delta=4.9303806576313238e-32
  13
  14 b = sqrt(2*m*E)/hbar;
  15 f(x,A,delta,b)=A*sin(b*x-delta);
  16 g(x,A,delta,b)=A*b*x*(cos(delta)*besj1(b*x)-sin(delta)*besy1(b*x));
  17 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0:3000:0 via b
  18 #fit f(x,A,delta,b) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta,b
  19 #fit g(x,A,delta,1) "tmp/noumerov.stdout" using 1:2 every ::2000:0::0 via A,delta
  20 plot    "tmp/noumerov.stdout"   using 1:2 every 2::2000:0::0,\
  21         ""                      using 1:2 every 2::2001:0::0,\
  22         f(x,A,delta,b) axis x1y2,\
  23         g(x,A,delta,1) axis x1y2
  24 pause -1
  25 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  26         "" using 0:2 every :::2::2 axis x1y2
  27 pause -1
  28 plot "tmp/noumerov.stdout" using 0:1 every :::2::2,\
  29         "" using 0:4 every :::2::2 axis x1y2
  30 pause -1
  31 exit;
  32 plot [0:5E-16] "tmp/noumerov.stdout" using 1:2 every :::0::0,\
  33         "" using 1:3 every :::0::0,\
  34         "" using 1:4  every :::0::0 axis x1y2,\
  35         "" using 1:5  every :::0::0 axis x1y2
  36 pause -1
  37 plot "tmp/noumerov.stdout" using 1:2 every :::1::1,\
  38         "" using 1:($2+$4) every :::1::1,\
  39         "" using 1:($2+$5) every :::1::1,\
  40         "" using 1:($2+$6) every :::1::1
  41 pause -1
  42 plot "tmp/noumerov.stdout",\
  43         x*exp(-x**2/2)
  44 pause -1
  45 plot [0:1E-14] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
  46         "" using 1:3 every :::2::2 axis x1y2,\
  47         "" using 1:4 every 2::1:2::2
  48 pause -1
  49 x0 = 1E-11;
  50 plot [x0+0:x0+1E-10] "tmp/noumerov.stdout" using 1:2 every 2:::2::2,\
  51         "" using 1:4 every :::2::2 axis x1y2,\
  52         "" using 1:5 every 2::1:2::2
  53 pause -1