VB中定义Arccos函数

游侠无极限

我是这样的：
Private Function Arccos(x As Double, n As Long)
# z! W3 I. N7 }; PDim i As Long, y As Double
+ L  I" X: i8 q) oFor i = 1 To n
    y = y + (-1 / Sqr(1 + (x / n * i) ^ 2)) * x / n
& h: I) ^! W( B' A% Y( [Next
Arccos = y
, I( c- k* w% ~& j( h/ wEnd Function

用的是Arccos的微分来算近似值，怎么不对？
; @* U7 O3 b# g# j) I# Y+ ?

煎饼

公式对就没有问题拉

游侠无极限

Const pi As Double = 3.14159265359
Private Function Arccos(x As Double, n As Long) ' n 为函数精度，建议为100
9 W- J8 m2 T$ m% L4 x  m6 B    Select Case x
3 g6 `, e" }* J  U) n" u    Case 1
        Arccos = 0
    Case -1
% ^; V: s) Z4 t; E& r5 V- I( z        Arccos = pi
    Case -1 To 1
! c9 c1 @1 K' ^$ j0 }, D+ T        Dim dx As Double, y As Double, i As Long
$ U' f3 P/ ?2 S  Q2 @7 l        dx = x / n
- O" c9 {$ N5 a/ a5 ~- m, U        For i = 1 To n
            y = y + (-1 / Sqr(1 - (dx * i + dx / 2) ^ 2)) * dx
        Next
6 Q5 }) Z1 z  I5 l. |        Arccos = y + pi / 2
' a' |8 M& f: e! E    Case Else
        '异常
    End Select
0 M& t+ h/ B- P+ S! \% u) ^# m; gEnd Function

搞定了