Engel expansion: Difference between revisions
Added Quackery.
SqrtNegInf (talk | contribs) m (→{{header|Perl}}: count of terms in expansion) |
(Added Quackery.) |
||
Line 379:
Number of terms : 34, places : 6 (6 correct)
Back to rational: 25.628906
</pre>
=={{header|Quackery}}==
<syntaxhighlight lang="Quackery"> [ $ "bigrat.qky" loadfile ] now!
[ /mod 0 != + ] is ceiling ( n/d --> n )
[ [] unrot
[ 2dup 1/v ceiling
dip rot
dup dip
[ join unrot ]
1 v* 1 1 v-
2dup v0= until ]
2drop ] is v->engel ( n/d --> [ )
[ 0 1 rot
1 1 rot
witheach
[ n->v v/
2swap 2over v+
2swap ]
2drop ] is engel->v ( [ --> n/d )
$ "3.14159265358979 2.71828182845904 1.414213562373095"
nest$
witheach
[ $->v drop
2dup 200 point$ echo$ cr
v->engel dup echo cr
engel->v 200 point$ echo$ cr
cr ]
$ "3.1415926535897932384626433832795028841971693993751058"
$ "209749445923078164062862089986280348253421170679821480" join
$ "865132823066470938446095505822317253594081284811174502" join
$ "841027019385211 " join
nested
$ "2.7182818284590452353602874713526624977572470936999595"
$ "7496696762772407663035354759457138217852516642743 " join
nested join
$ "1.4142135623730950488016887242096980785696718753769480"
$ "731766797379907324784621070388503875343276415727350138" join
$ "462309122970249248360558507372126441214970999358314132" join
$ "226659275055927558" join
nested join
witheach
[ $->v drop
2dup 200 point$ echo$ cr
v->engel engel->v
200 point$ echo$ cr
cr ]</syntaxhighlight>
{{out}}
<pre>3.14159265358979
[ 1 1 1 8 8 17 19 300 1991 2768 4442 4830 10560 37132 107315 244141 651042 1953125 ]
3.14159265358979
2.71828182845904
[ 1 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 17 82 144 321 2289 9041 21083 474060 887785 976563 1953125 ]
2.71828182845904
1.414213562373095
[ 1 3 5 5 16 18 78 102 120 144 260 968 18531 46065 63005 65105 78125 ]
1.414213562373095
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211
2.71828182845904523536028747135266249775724709369995957496696762772407663035354759457138217852516642743
2.71828182845904523536028747135266249775724709369995957496696762772407663035354759457138217852516642743
1.4142135623730950488016887242096980785696718753769480731766797379907324784621070388503875343276415727350138462309122970249248360558507372126441214970999358314132226659275055927558
1.4142135623730950488016887242096980785696718753769480731766797379907324784621070388503875343276415727350138462309122970249248360558507372126441214970999358314132226659275055927558
</pre>
| |||