Clear[E0, Eb, \[Lambda], t]
h = 6.58*10^(-13)
m = 0.657
a = 0.068
W1 = 5.497*10^(13)
Z = 1.6*10^(-22)
wc = 1*10^13
E0 = 3/2*h*((v0*Z)/(3*m*r^2) + wc^2/6)^(1/2) - v0 - a*h*W1
Partition[
Flatten[Table[{wc, r, v0, E0, Eb, \[Lambda], t}, {r, 0.5, 1.5,
0.1}, {v0, 1, 2, 0.1}]], 7]；
h = 6.58*10^(-13)
m = 0.657
a = 0.068
W1 = 5.497*10^(13)
Z = 1.6*10^(-22)
wc = 1*10^13
E0 = 3/2*h*((v0*Z)/(3*m*r^2) + wc^2/6)^(1/2) - v0 - a*h*W1
Partition[
Flatten[Table[{wc, r, v0, E0, Eb, \[Lambda], t}, {r, 0.5, 1.5,
0.1}, {v0, 1, 2, 0.1}]], 7]；
