Clear[E0]
e = 0.5
\[Lambda] = 3
w1 = 6
w2 = 4
a = 8
r = 3
E0 = (1 + e^2/2)*\[Lambda]^2 + w1^2/(4*\[Lambda]^2) + w2^2/(
8*\[Lambda]^2*e^2) -
Sqrt[2/\[Pi]]*(a*\[Lambda]*r*ArcSin[Sqrt[1 - e^2]])/Sqrt[1 - e^2] -
Sqrt[\[Pi]]/(2*\[Lambda])*F
Partition[
Flatten[Table[{w1, w2, e, \[Lambda], a, r, F, E0}, {F, 1, 3, 0.1}]],
8];