(* Combination of several programs -
    calculates diabatic/adiabatic exact/primitive potentials/
        wave functions along with "uniform" convolutions
      ALL examples *)
ClearAll["@"];
Get["C:\\SERGEEV\\Tulane\\math\\exactqua\\examples\\examples.m"];
mexamp = Length[examples];
names = Table[
      examples[[nexamp]];
      name, {nexamp, mexamp}];
Print["Examples: ", names // TableForm];
Do[
    examples[[nexamp]];
    x0 = x1 = 0;
    Get["C:\\sergeev\\tulane\\math\\exactqua\\exactwf1.m"];
    gr = GraphicsArray[{{pltv, pltw},
          {pltdc1, pltac1}, {pltdc2, pltac2}}];
    Show[gr, ImageSize -> 600,
      DisplayFunction -> $DisplayFunction];
    (* Calculating convolutions *)
    gr = GraphicsArray[{plotgamma[1], plotgamma[10], plotgamma[100], 
          plotgamma[1000]}];
    Show[gr, ImageSize -> 600,
      DisplayFunction -> $DisplayFunction];
    (* Calculating P12(E) *)
    denmin = 0.01;
    den = 10.;
    Plot[p12sc[ene], {ene, enmin + denmin, enmin + den}, ImageSize -> 600, 
      PlotRange -> All,
      DisplayFunction -> $DisplayFunction];
    Plot[p12exact[ene] - p12sc[ene], {ene, enmin + denmin, enmin + den}, 
      ImageSize -> 600, PlotRange -> All,
      DisplayFunction -> $DisplayFunction];
    , {nexamp, mexamp}];