Clear[alfa,gama]
alfa=1./137.036
gama=(alfa/8)^2*13.6

E1=Function[z,-gama +2 z gama]
E2=Function[z,-gama -2 z gama]
E3=Function[z,-3 gama + z gama/2 + Sqrt[4 gama^2 + 2/3*gama^2*z 
					+ gama^2*z^2/4]]
E4=Function[z,-3 gama + z gama/2 - Sqrt[4 gama^2 + 2/3*gama^2*z 
					+ gama^2*z^2/4]]
E5=Function[z,-3 gama - z gama/2 + Sqrt[4 gama^2 - 2/3*gama^2*z 
					+ gama^2*z^2/4]]
E6=Function[z,-3 gama - z gama/2 - Sqrt[4 gama^2 - 2/3*gama^2*z 
					+ gama^2*z^2/4]]
ans = Table[{z, E1[z], E2[z], E3[z], E4[z], E5[z], E6[z]}, {z, 0, 10, 0.01}];

stmp = OpenWrite["prob7exacto.dat",PageWidth -> 200]
Do[Write[stmp, FortranForm[ans[[i, 1]]], "  ", FortranForm[ ans[[i, 2]]], "  ",  FortranForm[ans[[i, 3]]], "  ", FortranForm[ ans[[i, 4]]], "  ",  FortranForm[ans[[i, 5]]], "  ",  FortranForm[ans[[i, 6]]], "  ", FortranForm[ans[[i, 7]]]], {i, 1, Length[ans]}]
Close[stmp]


(* Campo Fraco *)

E32BFraco1=Function[z,-gama + z*gama*4/3*3/2 ]
E32BFraco2=Function[z,-gama + z*gama*4/3*1/2 ]
E32BFraco3=Function[z,-gama - z*gama*4/3*1/2 ]
E32BFraco4=Function[z,-gama - z*gama*4/3*3/2 ]

E12BFraco1=Function[z,-5*gama + z*gama*2/3*1/2 ]
E12BFraco2=Function[z,-5*gama - z*gama*2/3*1/2 ]


ans = Table[{z, E32BFraco1[z], E32BFraco2[z], E32BFraco3[z], E32BFraco4[z], E12BFraco1[z], E12BFraco2[z]}, {z, 0, 4, 0.1}];


stmp = OpenWrite["prob7BFraco.dat",PageWidth -> 200]
Do[Write[stmp, FortranForm[ans[[i, 1]]],"  ", FortranForm[ ans[[i, 2]]], "  ",  FortranForm[ans[[i, 3]]], "  ", FortranForm[ ans[[i, 4]]], "  ",  FortranForm[ans[[i, 5]]], "  ",  FortranForm[ans[[i, 6]]], "  ", FortranForm[ans[[i, 7]]]], {i, 1, Length[ans]}]
Close[stmp]


(* Campo Forte *)

EBForte=Function[{z,ml,ms},gama/3*(-7 + 8*ml*ms) + z*gama*(ml+2*ms)]


ans = Table[{z, EBForte[z,1,1/2], EBForte[z,-1,-1/2], EBForte[z,0,1/2], EBForte[z,1,-1/2], EBForte[z,-1,1/2], EBForte[z,0,-1/2]}, {z, 1, 10, 0.1}];


stmp = OpenWrite["prob7BForte.dat",PageWidth -> 200]
Do[Write[stmp, FortranForm[ans[[i, 1]]],"  ", FortranForm[ ans[[i, 2]]], "  ",  FortranForm[ans[[i, 3]]], "  ", FortranForm[ ans[[i, 4]]], "  ",  FortranForm[ans[[i, 5]]], "  ",  FortranForm[ans[[i, 6]]], "  ", FortranForm[ans[[i, 7]]]], {i, 1, Length[ans]}]
Close[stmp]


