> with(orthopoly,H);

                                       [H]
> b := sqrt(2*n+1) ;

                                              1/2
                                b := (2 n + 1)
> k := sqrt(b / (2^n * n! * sqrt(Pi)));\


                                  /         1/2\1/2
                                  |(2 n + 1)   |
                             k := |------------|
                                  |  n      1/2|
                                  \ 2  n! Pi   /
> psi:= k * H(n, b*x) * exp(-b^2 * x^2 / 2);

             /         1/2\1/2
             |(2 n + 1)   |                  1/2                         2
      psi := |------------|    H(n, (2 n + 1)    x) exp(- 1/2 (2 n + 1) x )
             |  n      1/2|
             \ 2  n! Pi   /
> pqu := psi * psi;\


                       1/2               1/2   2                      2 2
              (2 n + 1)    H(n, (2 n + 1)    x)  exp(- 1/2 (2 n + 1) x )
       pqu := -----------------------------------------------------------
                                       n      1/2
                                      2  n! Pi
> n := 13;

                                     n := 13
> pcl := 1/(Pi * sqrt(1-x^2));

                                            1
                              pcl := --------------
                                              2 1/2
                                     Pi (1 - x )
> interface(plotdevice=hpgl, plotoutput=HPGLherm); 
> plot({pcl,pqu},x=-2..2,title='Christian');
> interface(plotdevice=x11);
