Public Function Ellipic1(d)'第1種完全楕円積分
  PID = 3.14159265358979 / 2
  C = 1#: dd = 1: S = 0: i1 = 1: i2 = 2
  If d >= 0.999999 Then
     Ellipic1 = 1E+308
     Exit Function
  End If
  For i = 1 To 100000
    SS = S: S = S + C * C * dd
    C = C * i1 / i2: dd = dd * d * d
    i1 = i1 + 2: i2 = i2 + 2
    If Abs(S - SS) < 0.000000001 Then Exit For
  Next
  Ellipic1 = S * PID
End Function
Public Function Ellipic2(d) '第2種完全楕円積分
  PID = 3.14159265358979 / 2
  C = 1#: dd = d * d: S = 1: i1 = 1: i2 = 2
  If d >= 0.999999 Then
     Ellipic2 = 1
     Exit Function
  ElseIf d < 0.000000001 Then
     Ellipic2 = PID
     Exit Function
  End If
  For i = 1 To 100000
    C = C * i1 / i2
    SS = S: S = S - C * C * dd / i1
    dd = dd * d * d
    i1 = i1 + 2: i2 = i2 + 2
    If Abs(S - SS) < 0.000000001 Then Exit For
  Next
  Ellipic2 = S * PID
End Function