Beefy Boxes and Bandwidth Generously Provided by pair Networks
Welcome to the Monastery

Approximation with chain fractions

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


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
Download Code
Re: Approximation with chain fractions
by Anonymous Monk on Jan 24, 2005 at 19:31 UTC
    Can you add code tags please?

Back to Snippets Section

Log In?

What's my password?
Create A New User
Node Status?
node history
Node Type: snippet [id://424519]
and the web crawler heard nothing...

How do I use this? | Other CB clients
Other Users?
Others lurking in the Monastery: (7)
As of 2015-07-02 04:02 GMT
Find Nodes?
    Voting Booth?

    The top three priorities of my open tasks are (in descending order of likelihood to be worked on) ...

    Results (27 votes), past polls