Perl Monk, Perl Meditation PerlMonks

### Approximation with chain fractions

by ambrus (Abbot)
 on Jan 24, 2005 at 07:55 UTC 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; }
```
Replies are listed 'Best First'.
Re: Approximation with chain fractions
by Anonymous Monk on Jan 24, 2005 at 19:31 UTC

Create A New User
Node Status?
node history
Node Type: snippet [id://424519]
help
Chatterbox?
and God said, "Let Newton be!"...

How do I use this? | Other CB clients
Other Users?
Others surveying the Monastery: (7)
As of 2018-03-20 08:43 GMT
Sections?
Information?
Find Nodes?
Leftovers?
Voting Booth?
When I think of a mole I think of:

Results (248 votes). Check out past polls.

Notices?