求大神帮忙看看程序错在哪，谢谢啦

md = 0.05*10^-3
L = 3
r0 = 0.5
ri = 0.1
\[Mu]r = 1
\[Epsilon]r = 1
d = 0.02
Z0 = Zc0 = Sqrt[\[Mu]r/\[Epsilon]r]*60*Log[r0/ri]
\[Epsilon]0 = 8.854*10^-12
\[Mu]0 = 4*Pi*10^-7
n = 10
f = 1.373*10^6
\[Beta] = 2*Pi*f*Sqrt[\[Epsilon]0*\[Epsilon]r*\[Mu]0*\[Mu]r]
L0 = md/n
For[i = 2, i < n + 1,
{p = Solve[{x1*x2 == r0^2, (x1 - md/n*i)*(x2 - md/n*i) == 0.01,
x1 < x2}, {x1, x2}],
x11 = p[[1, 1, 2]], x22 = p[[1, 2, 2]],
r = Abs[((r0 - x22)*(ri + md/n*i - x11))/((r0 - x11)*(ri + md/n*i -
x22))],
Zc = Sqrt[\[Mu]r/\[Epsilon]r]*60*Log[1/r],
L = md/n + 1/\[Beta] ArcTan[Z0/Zc*\[Beta]*L0], Z0 = Zc, L0 = L}, i++;
Print[L]]