Keep It Simple, Stupid PerlMonks

### Approximation with chain fractions

by ambrus (Abbot)
 on Jan 24, 2005 at 07:55 UTC ( #424519=snippet: print w/ replies, xml ) Need Help??

 Description: Approximate any real number \$x with chain fractions. I'm not the first one offering code on this, see Continued Fractions. ```#use bignum "a", 32; # uncomment for greater precision \$x = 2**(1/12); # change this to any number you like print \$x, "\n"; (\$p, \$q, \$r, \$s) = (0, 1, 1, 0); { \$n = int((\$p - \$x*\$q) / (\$x*\$s - \$r)); (\$p, \$q) = (\$p + \$n*\$r, \$q + \$n +*\$s); print \$p, "/", \$q, " = ", \$y = \$p/\$q, "\n"; \$y == \$x and last; \$n = int((\$r - \$x*\$s) / (\$x*\$q - \$p)); (\$r, \$s) = (\$r + \$p*\$n, \$s + \$q +*\$n); print \$r, "/", \$s, " = ", \$y = \$r/\$s, "\n"; \$y == \$x and last; redo; } ```
Comment on Approximation with chain fractions
Replies are listed 'Best First'.
Re: Approximation with chain fractions
by Anonymous Monk on Jan 24, 2005 at 19:31 UTC

Back to Snippets Section

Create A New User
Node Status?
node history
Node Type: snippet [id://424519]
help
Chatterbox?
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others about the Monastery: (14)
As of 2015-11-24 22:52 GMT
Sections?
Information?
Find Nodes?
Leftovers?
Voting Booth?