\newcommand{\vt}[1]{\boldsymbol{#1}}
\newcommand{\uv}[1]{\hat{\boldsymbol{#1}}}
\newcommand{\tv}[1]{\tilde{\boldsymbol{#1}}}
\newcommand{\dvg}{\operatorname{div}}
\newcommand{\grd}{\operatorname{grad}}
using WiltonInts84
using FixedSizeArrays

using PyPlot

p1 = Vec(0.0, 0.0, 0.0)
p2 = Vec(1.0, 0.0, 0.0)
p3 = Vec(0.0, 1.0, 0.0)
h = 0.1

x = -0.5 : 0.01 : 1.5
y = -0.5 : 0.01 : 1.5

t, T = 0.15, 1
r, R = max(0,t-T), max(0,t)

M = zeros(typeof(t), length(x), length(y))
for i in eachindex(x)
    for j in eachindex(y)
        c = Vec(x[i],y[j],h)
        try
            I, K = wiltonints(p1,p2,p3,c,r,R,Val{0})
            M[i,j] = I[2]
        catch
            #@show i j
            #@show c
        end
    end
end

imshow(M); colorbar()