keskiarvojono.jl

function keskiarvojono(n,k)
# aluksi a= 1,2,..k
# sitten seuraava on k:n edellisen keskiarvo jne
# a[k+1] = 2,3,...k,ak jne
    a = rand(Float64,n) # arvotaan satunnaisjono, n lukua
    # alkuu lukujono 1,2,..,k
    for i = 1:k
        a[i] = i
    end
    # tulostetaan jonon alun k termiä
    for i= 1:k
        println(a[i])
    end
    for i = k+1 : n
        a[i] = sum(a[i-k:i-1])/k # edellisten k termin keskiarvo
        println(a[i]) # keskiarvotermi
    end
    println("jonon teoreettinen raja-arvo = (2k+1)/3 = ",(2*k+1)/3)
end