Almkvist-Giullera formula for pi: Difference between revisions
Content added Content deleted
Thundergnat (talk | contribs) m (→{{header|C#|Csharp}}: Regularize header markup to recommended on category page) |
(add PicoLisp) |
||
Line 1,813: | Line 1,813: | ||
Pi to 70 d.p.: 3.1415926535897932384626433832795028841971693993751058209749445923078164 |
Pi to 70 d.p.: 3.1415926535897932384626433832795028841971693993751058209749445923078164 |
||
Pi (builtin) : 3.1415926535897932384626433832795028841971693993751058209749445923078164 |
Pi (builtin) : 3.1415926535897932384626433832795028841971693993751058209749445923078164 |
||
</pre> |
|||
=={{header|PicoLisp}}== |
|||
<lang PicoLisp>(scl 70) |
|||
(de fact (N) |
|||
(if (=0 N) |
|||
1 |
|||
(* N (fact (dec N))) ) ) |
|||
(de almkvist (N) |
|||
(let |
|||
(A (* 32 (fact (* 6 N))) |
|||
B (+ (* 532 N N) (* 126 N) 9) |
|||
C (* (** (fact N) 6) 3) ) |
|||
(/ (* A B) C) ) ) |
|||
(de integral (N) |
|||
(*/ |
|||
1.0 |
|||
(almkvist N) |
|||
(** 10 (+ 3 (* 6 N))) ) ) |
|||
(let (S 0 N -1) |
|||
(do 10 |
|||
(println (inc 'N) (almkvist N)) ) |
|||
(prinl) |
|||
(setq N -1) |
|||
(while (gt0 (integral (inc 'N))) |
|||
(inc 'S @) ) |
|||
(setq S (sqrt (*/ 1.0 1.0 S) 1.0)) |
|||
(prinl "Pi to 70 decimal places is:") |
|||
(prinl (format S *Scl)) )</lang> |
|||
{{out}} |
|||
<pre> |
|||
0 96 |
|||
1 5122560 |
|||
2 190722470400 |
|||
3 7574824857600000 |
|||
4 312546150372456000000 |
|||
5 13207874703225491420651520 |
|||
6 567273919793089083292259942400 |
|||
7 24650600248172987140112763715584000 |
|||
8 1080657854354639453670407474439566400000 |
|||
9 47701779391594966287470570490839978880000000 |
|||
Pi to 70 decimal places is: |
|||
3.1415926535897932384626433832795028841971693993751058209749445923078152 |
|||
</pre> |
</pre> |
||