Kompleksifunktionjuuret.jl

# Kompleksifunktion juuret
# Riemannin zeta-funktion nollakohdat
# Juhani Kaukoranta 5.2.2023
using RootsAndPoles,SpecialFunctions

function solverootspoles(func,xalku,xloppu,yalku,yloppu)
    # laskee kompelsifunktio nollakohdat halutulta väliltä
    # xalku,xloppu = etsintävälin realiosan rajat ala,ylä
    # yalku,yloppu = etsintävälin imag osan rajat ala,ylä
    r = 0.1  # etsintähilan askel step
    origcoords = rectangulardomain(complex(xalku, yalku), complex(xloppu, yloppu), r);
    zroots, zpoles = grpf(func, origcoords);
    return zroots # jätetään tulostuksesta zpoles pois
end

function zetazeros(alku,loppu) # zeta-funktion imaginaariosan nollakohtia
    zroots = solverootspoles(zeta,0.49,0.51,alku,loppu) # reaaliosa aina 0.5
    println(sort(imag(zroots))) # tulostetaan imaginaariosat
end
zetazeros(10,30) # zeta-funktion imaginaariosan nollakohdat välillä 10-30