素の Emacs で Calc の多倍長整数を使って円周率を1000桁計算する方法

calc-pi.el

;; 出力結果は文字列だけど、3の後ろに小数点が無いのと最後の4桁は未収束の為除いたものが正しい値

(defmath mpiarccot (x)
  (let ((base (expt 10 (+ 1000 4)))
        (a x)
        (b 1)
        pos
        (sum 0))
    (while (> (setq pos (idiv base (* a b))) 1)
      (setq sum (+ sum (if (/= (logand 2 b) 0) (- pos) pos))
            a (* a x x)
            b (+ b 2)))
    sum))

(calc-eval "4 * (4 * mpiarccot(5) - mpiarccot(239))") ; => "3141592653...42019893644"