@testset "nodahilbert.jl" begin
    @test N(2) ≈ [0.0                0.3183098861837907;
                 -0.3183098861837907 0.0               ]
    @test N(3) ≈ [0.0                  0.3183098861837907 0.15915494309189535;
                 -0.3183098861837907   0.0                0.3183098861837907 ;
                 -0.15915494309189535 -0.3183098861837907 0.0                ]
    @test N(4) ≈ [0    1    1/2 1/3;
                 -1    0    1   1/2;
                 -1/2 -1    0   1  ;
                 -1/3 -1/2 -1   0  ] / π

    function f(i, j)
        if i == j
            return 0.0
        end
        return 1 / (π * (j - i))
    end
    for n in 5:10
        @test N(n) ≈ [f(i, j) for i in 1:n, j in 1:n]
    end

    @test N * [1, 2] ≈ [0.6366197723675814, -0.3183098861837907]
    @test N * [1, 2, 3] ≈ [1.1140846016432673, 0.6366197723675813, -0.7957747154594768]
end
