n20(1) = n(5) / n(19)

n22(1) = n20(1) / Sqr(1 + n20(1) ^ 2)

n24(1) = (n(5) - n(17) * n22(1)) / n(12)

n25(1) = n24(1) / Sqr(1 - n24(1) ^ 2)

n26(1) = n22(1) / n25(1)

n27(1) = 1 / Sqr(1 + n26(1) ^ 2)

n28(1) = n24(1) / n27(1)

n29(1) = Sqr(1 - n28(1) ^ 2)

n30(1) = (n(15) - n29(1)) / n28(1)

n31(1) = n28(1) * (n(9) - n30(1))

n32(1) = n(3) * n31(1)

n33(1) = n24(1) - (n22(1) * n32(1))

n34(1) = n33(1) / Sqr(1 - n33(1) ^ 2)

n35(1) = n22(1) / n34(1)

n36(1) = deg(Atn(n35(1)))

n37(1) = Cos(Atn(n35(1)))

n38(1) = n33(1) / n37(1)

n39(1) = deg(Atn((n38(1) / Sqr(1 - n38(1) ^ 2))))

n40(1) = Sqr(1 - n38(1) ^ 2)

n41(1) = (n(15) + n31(1) - n40(1)) / n38(1)

n42(1) = deg(Atn(n41(1)))

n43(1) = Cos(Atn(n41(1)))

n44(1) = n42(1) - n39(1)

n45(1) = Cos(rad(n44(1)))

n46(1) = Tan(rad(n44(1)))

n47(1) = n20(1) / n33(1)

n48(1) = deg(Atn(1 / n47(1)))

n49(1) = Sin((Atn(1 / n47(1))))

n50(1) = Cos((Atn(1 / n47(1))))

n51(1) = (n(17) + n(12) * n32(1)) / n37(1)

n52(1) = n(12) / n50(1)

n53(1) = n51(1) + n52(1)

n54(1) = n(12) * n45(1) / n49(1)

n55(1) = n43(1) * n51(1) / n35(1)

n56(1) = (n41(1) * n55(1) - n46(1) * n54(1)) / n53(1)

n57(1) = deg(Atn(n56(1)))

n58(1) = Cos(Atn(n56(1)))

n59(1) = (n41(1) * n56(1)) / n51(1)

n60(1) = (n46(1) * n56(1)) / n52(1)

n61(1) = n54(1) * n55(1)

n62(1) = (n54(1) - n55(1)) / n61(1)

n63(1) = n59(1) + n60(1) + n62(1)

n64(1) = (n41(1) - n46(1)) / n63(1)

n65(1) = n64(1) / n58(1)

n66(1) = n(7) / n65(1)
n20(1)项是个变量
如果1-n66(1)>0.001,那么n20(2)=0.9*n20(1)，然后按以上公式在计算一次，如果还不满足条件，20项按下式迭代

x=x+1

n20(X) = ((n20(X - 1) - n20(X - 2)) / (n66(X - 1) - n66(X - 2))) * (1 - n66(X - 2)) + n20(X - 2)

直到满意为止

也就是说20项是个变量，其他公式一样，不满足条件就按迭代式变化20项在带入以上的公式计算，直到满意为止

各位大侠，拜求这个程序迭代代码，谢谢

[此贴子已经被作者于2016-1-8 20:56编辑过]

各位大侠，请帮忙，谢谢