Continued fraction: Difference between revisions

{{trans|ALGOL_68}}

<lang Phix>function continued_fraction(integer steps, fa, fb)

atom res = 0

for n=steps to 1 by -1 do

res := fb(n) / (fa(n) + res)

end for

return fa(0) + res

end function

function sqr2_a(integer n) return iff(n=0?1:2) end function

function sqr2_b(integer n) return 1 end function

function nap_a(integer n) return iff(n=0?2:n) end function

function nap_b(integer n) return iff(n=1?1:n-1) end function

function pi_a(integer n) return iff(n=0?3:6) end function

function pi_b(integer n) return iff(n=1?1:power(2*n-1,2)) end function

constant precision = 10000

printf(1,"Sqr(2): %.10g\n", {continued_fraction(precision, sqr2_a, sqr2_b)})

printf(1,"Napier: %.10g\n", {continued_fraction(precision, nap_a, nap_b)})

printf(1,"Pi: %.10g\n", {continued_fraction(precision, pi_a, pi_b)})</lang>