(* Testing results of Segev - Heller paper *)
(* Main program *)

(* Input:  x0, y0 - displacements
           omx, omy - frequencies
           nmax - quantum number of the acceptor state
Example of input line: {x0, y0, omx, omy} = {1.1519, 1.2649, 0.4689, 0.5555};nmax=10;
*)

(* Printing input parameters and checking their validity *)
<<harmchec.txt;

(* Quantum results *)
<<harmquan.txt;

(* Phase space results (constrained minimium of the Wigner function) *)
<<harmwign.txt;

(* Graphs of potentials in three dimensions *)
<<harmpots.txt;

(* Density diagram of the wave function of the final state *)
<<harmwfun.txt;

(* Distribution of the energy over two modes *)
<<harmener.txt;