(* Testing results of Segev - Heller paper *)
(* Printing input parameters and checking their validity *)

(* 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;
*)

(* Load a package for formatted printing of numbers *)
<<printing.txt;

fails=False;
Print["<PRE><B>Frequencies:</B>    <B style=\"color:black;background-color:#DAFCFE\">omx =",printF[omx,10,5],"</B>,   <B style=\"color:black;background-color:#DAFCFE\">omy =",printF[omy,10,5],"</B>"];
If[!(NumberQ[omx]&&NumberQ[omy]&&omx>0&&omy>0),
Print["<FONT COLOR=\"#FF0000\">Frequencies must be positive numbers!</FONT>\n"];
fails=True];
Print["<B>Displacements:</B>  <B style=\"color:black;background-color:#DAFCFE\">x0  =",printF[x0,10,5],"</B>,   <B style=\"color:black;background-color:#DAFCFE\">y0  =",printF[y0,10,5],"</B>"];
If[!(NumberQ[x0]&&NumberQ[y0]),
Print["<FONT COLOR=\"#FF0000\">Displacements must be numbers!</FONT>\n"];
fails=True];
Print["<B>Quantum number in the acceptor potential:</B> <B style=\"color:black;background-color:#DAFCFE\">n = ",nmax,"</B>"];
If[!(IntegerQ[nmax]&&nmax>=0),
Print["<FONT COLOR=\"#FF0000\">Quantum number in the acceptor potential must be a non-negative integer!</FONT>\n"];
fails=True];
If[!(nmax<=50),
Print["<FONT COLOR=\"#FF0000\">Quantum number in the acceptor potential is too large!</FONT>\n"];
fails=True];
If[fails,
   Print["Computations terminated because input parameters are illegitimate.</PRE>"];
   Exit[]
  ];
(* Testing ends *)
If[!(omx<100&&omy<100),
Print["<FONT COLOR=\"#0000FF\">Warning: one of frequencies is too large. Results may be unpredictable.</FONT>"];
];
If[!(omx>0.005&&omy>0.005),
Print["<FONT COLOR=\"#0000FF\">Warning: one of frequencies is too small. Results may be unpredictable.</FONT>"];
];
If[!(Abs[x0]<100&&Abs[y0]<100),
Print["<FONT COLOR=\"#0000FF\">Warning: one of displacements is too large. Results may be unpredictable.</FONT>"];
];
If[nmax<=20,
   {x0, y0, omx, omy}=N[{x0, y0, omx, omy}],
   {x0, y0, omx, omy}=SetPrecision[{x0, y0, omx, omy},128]
  ];
Print["</PRE>"];