(*H2: the results of RHF, RMP2, RMP3 and CISD within the  STO-3G atomic basis set*)

(*I. RHF*)

dataRHF={{0.5,-1.04300},{0.6,-1.10113},{0.7,-1.11735},{0.8,-1.11085},{0.9,-1.09191},{1.0,-1.06611},{1.2,-1.00511},{1.4,-0.94148},{1.6,-0.88173},{2.0,-0.78379},{2.4,-0.71591},{3.0,-0.65605}} (*H2: distances and  RHF_STO_3G energies*)

crveRHF=Fit[dataRHF, {1,Exp[-x],Exp[-2*x],Exp[-3*x],Exp[-4*x],Exp[-5*x]},x] (*creation of the curve crveRHF that approximates  dataRHF by using a four-exponent expansion*)

Show[ListPlot[dataRHF, PlotStyle->Red], Plot[crveRHF, {x, 0., 4.},PlotStyle->{Black,Thick}]] (*plot of  crveRHF with points*)

FindMinimum[crveRHF,{x,0.7}] (*glebokosc bezwzgledna i pozycja minimum*)

(*II. RMP2*)

dataRMP2={{0.5,-1.05149},{0.6,-1.11133},{0.7,-1.12958},{0.8,-1.12545},{0.9,-1.10927},{1.0,-1.08667},{1.2,-1.03366},{1.4,-0.98057},{1.6,-0.93424},{2.0,-0.87251},{2.4,-0.85285},{3.0,-0.88143}} (*H2: odleglosci i energie RMP2_STO_3G*)

crveRMP2=Fit[dataRMP2, {1,Exp[-x],Exp[-2*x],Exp[-3*x],Exp[-4*x],Exp[-5*x]},x] (*tworzenie funkcji crveRMP2 aproksymujacej dataRMP2 przez rozwiniecie czteroeksponencjalne*)

Show[ListPlot[dataRMP2, PlotStyle->Red], Plot[crveRMP2, {x, 0., 4.}]] (*wykres crveRMP2 z punktami*)

FindMinimum[crveRMP2,{x,0.7}] (*glebokosc bezwzgledna i pozycja minimum*)

(*III. RMP3*)

dataRMP3={{0.5,-1.05408},{0.6,-1.11472},{0.7,-1.13396},{0.8,-1.13105},{0.9,-1.11635},{1.0,-1.09555},{1.2,-1.04738},{1.4,-1.00128},{1.6,-0.96469},{2.0,-0.93281},{2.4,-0.95873},{3.0,-1.08132}} (*H2: odleglosci i energie RMP3_STO_3G*)

crveRMP3=Fit[dataRMP3, {1,Exp[-x],Exp[-2*x],Exp[-3*x],Exp[-4*x],Exp[-5*x]},x] (*tworzenie funkcji crveRMP3 aproksymujacej dataRMP3 przez rozwiniecie czteroeksponencjalne*)

Show[ListPlot[dataRMP3, PlotStyle->Red], Plot[crveRMP3, {x, 0., 4.}]] (*wykres crveRMP3 z punktami*)

FindMinimum[crveRMP3,{x,0.7}] (*glebokosc bezwzgledna i pozycja minimum*)

(*IV. CISD*)

dataCISD={{0.5,-1.05516},{0.6,-1.11629},{0.7,-1.13619},{0.8,-1.13415},{0.9,-1.12056},{1.0,-1.10115},{1.2,-1.05674},{1.4,-1.01547},{1.6,-0.98347},{2.0,-0.94864},{2.4,-0.93725},{3.0,-0.93363}} (*H2: odleglosci i energie CISD_STO_3G*)

crveCISD=Fit[dataCISD, {1,Exp[-x],Exp[-2*x],Exp[-3*x],Exp[-4*x],Exp[-5*x]},x] (*tworzenie funkcji crveCISD aproksymujacej dataCISD przez rozwiniecie czteroeksponencjalne*)

Show[ListPlot[dataCISD, PlotStyle->Red], Plot[crveCISD, {x, 0., 4.}]] (*wykres crveCISD z punktami*)

FindMinimum[crveCISD,{x,0.7}] (*glebokosc bezwzgledna i pozycja minimum*)

(*V. All methods*)

Plot[{crveRHF,crveRMP2, crveRMP3,crveCISD},{x, 0., 3.}, PlotStyle->{Red, Green, Blue,Black}] (*wszystkie wykresy razem*)