dá sa to rýchlejšie so zaokrúhľovaním: pýtajme sa na jeden ťah a zaokrúhlime: round(13/14+15/16+1/12)=20:10 😉 prečo som ich sčítal? lebo chcem dosiahnuť nulu vo výsledku, á čo logicky musím z nuly odrátať alebo pridať? no tých 20:10 🙂 misia na konci, odchádzam o5 spať...

using Pkg
Pkg.add("Plots")
Pkg.add("SymPy")

using SymPy
using Plots

x = symbols("x")
ikses = Sym[]  # na uloženie riešení
speculs = Sym[]  # na uloženie špekulatívnych odhadov

for k in 1:10
    eq = (k*x - 13)/(k*x - 14) - (k*x - 15)/(k*x - 16) + 1//12
    sol = solve(eq, x)
    specul = rationalize(k*1/(13/14 + 15/16 + 1/12))

    for s in sol
        push!(ikses, s)
        push!(speculs, specul)
    end
end

# Vypočítať pomery susedných riešení
ratios = Sym[]
for k in 1:length(ikses)-1
    r = rationalize(round(N(ikses[k+1])) / round(N(ikses[k])))
    push!(ratios, r)
end
# Výpisy
println("speculs: ", speculs)
println("reals: ", ratios)
# Vykreslenie grafu
scatter(
    1:length(speculs), speculs, 
    orientation = :vertical, 
    title = "Speculation vs Real guess comparison", 
    label = "speculs", 
    xlabel = "iteration", 
    ylabel = "value"
)
scatter!(
    1:length(ratios), ratios,
    orientation = :vertical,
    label = "reals"
)
hline!([4], color=:green, linestyle=:dash, label="75% speculs")
vline!([14], color=:green, linestyle=:dash, label="75% correct")
display(current())
# Pre interaktívne prostredie
if !Base.isinteractive()
    println("Press enter to quit:")
    readline()
end

---