ketjumurto_to_rationals.jl

# ketjumurto_to_rationals.jl
# Juhani Kaukoranta 4.1.2023
A = [3, 7, 15, 1, 292, 1, 1, 1, 2, 1, 3, 1, 14, 2, 1, 1, 2, 2, 2, 2] # pi
#A = [1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2] # sqrt(2)
function rationals(A) # esitys rationaalilukuina
    println("ketjumurtoluku = ",A) # tulostaa ketjun
    println("ketjun pituus = ",length(A))
    h0,k0 = A[1],1 # ekan osoittaja ja nimittäjä
    h1,k1 = A[1]*A[2]+1,A[2] # tokan osoittaja ja nimittäjä
    println(h0,"/",k0," = ",h0/k0) # eka rational
    println(h1,"/",k1," = ",h1/k1) # toka rational
    for i in A[3:end] # ketjumurtoluku kolmosalkiosta lähtien
        h = i*h1 + h0
        k = i*k1 + k0
        println(h," / ",k," = ",h/k)
        h1,h0 = h,h1
        k1,k0 = k,k1
    end
end