function InfiniteHamil(E, EA, EAB, FA_list, FB_list, d_list, ϵ_list)
    #------------------------------
    num_steps = 200
    x_values  = collect(0:1/num_steps:1)
    y_values  = []

    #------------------------------
    for x in x_values
        k   = path(x)
        z1  = exp(im*k[1])
        z2  = exp(im*k[1]/2)*exp(im*sqrt(3)*k[2]/2)
        c1  = cos(k[1])
        c2  = cos(k[1]/2 + k[2]*sqrt(3)/2)
        c12 = cos(k[1]/2 - k[2]*sqrt(3)/2)

        secondHam = vcat( hcat(E+ 2*(c1*d_list[1] + c2*d_list[2]+ c12*d_list[5]),
                            FA_list[1]+ conj(z2)*(FA_list[2]+z1*FA_list[3]),
                            FB_list[1]+ conj(z2)*(FB_list[2]+z1*FB_list[3])),
                      hcat((FA_list[1]+ conj(z2)*(FA_list[2]+z1*FA_list[3]))',
                            EA+ 2*( c1*ϵ_list[1]+c2*ϵ_list[2]+c12*ϵ_list[5]),
                            EAB ),
                      hcat((FB_list[1]+ conj(z2)*(FB_list[2]+z1*FB_list[3]))',
                            EAB',
                            EA+ 2*( c1*ϵ_list[1]+c2*ϵ_list[2]+c12*ϵ_list[5]) ))


        println(typeof(E+ 2*(c1*d_list[1] + c2*d_list[2]+ c12*d_list[5])))
        println(typeof(FA_list[1]+ conj(z2)*(FA_list[2]+z1*FA_list[3])))
        println(typeof(FB_list[1]+ conj(z2)*(FB_list[2]+z1*FB_list[3])))
        println(typeof((FA_list[1]+ conj(z2)*(FA_list[2]+z1*FA_list[3]))'))
        println(typeof(EA+ 2*( c1*ϵ_list[1]+c2*ϵ_list[2]+c12*ϵ_list[5])))
        println(typeof(EAB ))
        println(typeof((FB_list[1]+ conj(z2)*(FB_list[2]+z1*FB_list[3]))'))
        println(typeof(EAB' ))
        println(typeof(EA+ 2*( c1*ϵ_list[1]+c2*ϵ_list[2]+c12*ϵ_list[5])))
        println(typeof(secondHam))

        push!(y_values ,  filter(x -> x > -10, sort(real(eigvals(secondHam))) ))

    end
    #------------------------------
    y_values = hcat(y_values...)
    plotData(x_values,y_values)
end