拜求迭代程序代码
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项按下式迭代
从第三次开始,以后公式一样,见下式
n20(X) = ((n20(X - 1) - n20(X - 2)) / (n66(X - 1) - n66(X - 2))) * (1 - n66(X - 2)) + n20(X - 2)
直到满意为止
也就是说20项是个变量,其他公式一样,不满足条件就按迭代式变化20项在带入以上的公式计算,直到满意为止
各位大侠,拜求这个程序迭代代码,谢谢