Beefy Boxes and Bandwidth Generously Provided by pair Networks
Welcome to the Monastery
 
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
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?
Username:
Password:

What's my password?
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 lurking in the Monastery: (7)
As of 2015-07-02 04:02 GMT
Sections?
Information?
Find Nodes?
Leftovers?
    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