Almkvist-Giullera formula for pi: Difference between revisions
Added Quackery.
m (added whitespace, aligned formulæ.) |
(Added Quackery.) |
||
Line 540:
π to 70 digits is 3.1415926535897932384626433832795028841971693993751058209749445923078164
mpmath π is 3.1415926535897932384626433832795028841971693993751058209749445923078164
</pre>
=={{header|Quackery}}==
<lang Quackery> [ $ "bigrat.qky" loadfile ] now!
[ 1 swap times [ i^ 1+ * ] ] is ! ( n --> n )
[ dup dup 2 ** 532 *
over 126 * + 9 +
swap 6 * ! * 32 *
swap ! 6 ** 3 * / ] is intterm ( n --> n )
[ dup intterm
10 rot 6 * 3 + **
reduce ] is vterm ( n --> n/d )
10 times [ i^ intterm echo cr ] cr
0 n->v
53 times [ i^ vterm v+ ]
1/v 70 vsqrt drop
70 point$ echo$ cr</lang>
{{out}}
<pre>96
5122560
190722470400
7574824857600000
312546150372456000000
13207874703225491420651520
567273919793089083292259942400
24650600248172987140112763715584000
1080657854354639453670407474439566400000
47701779391594966287470570490839978880000000
3.1415926535897932384626433832795028841971693993751058209749445923078164
</pre>
| |||