(* Finding the minimum of the function W under the energy constraint H = E *)
(* Main program *)

(* Input:  x0, y0 - displacements
           omx, omy - frequencies
           h30,h21,h12,h03 - cubic anharmonicity for H
           h40,h31,h22,h13,h04 - 4th degree anharmonicity for H
           w30,w21,w12,w03 - cubic anharmonicity for W
           w40,w31,w22,w13,w04 - 4th degree anharmonicity for W
Example of input line: {x0, y0, omx, omy} = {0.1, 0.2, 0.8, 1.2};
                       {h30,h21,h12,h03}={-0.1,-0.3,0.2,0};
                       {h40,h31,h22,h13,h04}={0,0,0,0,0};
                       {w30,w21,w12,w03}={0.05,0.05,0.05,-0.05};
                       {w40,w31,w22,w13,w04}={0,0,0,0,0};
*)

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

(* Finding extrema and choosing the global minimum *)
<<contnumb.txt;

(* Density plot of functions with display of the minimum *)
<<contplot.txt;