Der Harmonische Oszillator

1) Definitionen

Ψq[n_, 
          x_] := ((HermiteH[n, x]^2 (E^(-(x^2)  )))/(n! (Sqrt[π]) (2^n)));
c[x_, A_] := 
        If[Abs[x]  <  Abs[A], (1/(π A (Sqrt[ 1 - ((((x/A)))^2)]))), 
          0];
\(\(p[n_] := 
        Block[{xm = \@\(n +  .5\), ym, 
            A}, 
				\ ym\ \ \ \  = \ \(-\(FindMinimum[\(-Ψq[
                      n, z]\), {z, 
                    2\ xm}]\)[\([1]\)]\); 
				A\ \ \ \ \ \ \ \  = \@\(\(\ \)\(2  n\  + \ 1\)\); \  (*gleiche\ Varianz*) 
				Plot[{Ψq[n, x], c[x, A]}, {x, \(-1.4\)\ A, 1.4\ A}, 
            PlotRange -> {{\(-1.4\)\ A, 1.4\ A}, {0, 2\ ym}}, 
            PlotPoints -> 
              1000, 
						PlotStyle -> {{Thickness[ .005]}, {Thickness[0.005], 
                  Dashing[{0.01, 
                      0.01}]}}, 
					AxesLabel -> {\*""x/(Subscript[x, 0])"", "|Ψ(x)|^2,  p(x)"}, 
            TextStyle -> {FontSize -> 
                  16}]
			];\)\)

In Einheiten von (Subscript[x, 0]) !

2) Plots für das Skriptum

SetDirectory["C:\Dateien\wvl\Tex\Skripten\QM-Skript\M&P"];
pp = p[0];
Display["Graph0.eps", pp, "EPS"];
pp = p[1];
Display["Graph1.eps", pp, "EPS"];
pp =  p[2];
Display["Graph2.eps", pp, "EPS"];
pp = p[20];
Display["Graph20.eps", pp, "EPS"];

3) Aufenthaltswahrscheinlichkeiten

p[100];