Get["helium_po.m"]
(* Plot some trail trajectory *)
Table[PlotHelium[x, 5], {x, 0.5, 1.8, .1}]

(* Starting from some Guess for x and t find improved solution *)
l1 = CleanRoot[1.73, 1.0, 2, 0.5, 5]
l2 = CleanRoot[1.35, 3.2, 3.5, 0.2, 5]
l3 = CleanRoot[1.06, 4.8, 6, 0.05, 10]
l4 = CleanRoot[0.95, 7, 8, 0.02, 5]

(* Plot the improved solutions *)
ph1 = PlotHelium[l1]
ph2 = PlotHelium[l2]
ph3 = PlotHelium[l3]
ph4 = PlotHelium[l4]

(* Save graphs (note that you need to use psfix to convert the files *
 * to real postscript                                                *)
Display["gr1.mps", ph1]
Display["gr2.mps", ph2]
Display["gr3.mps", ph3]
Display["gr4.mps", ph4]

(* Calculate the action for the solutions *)
{Action[l1], Action[l2], Action[l3], Action[l4]}

