ngiann/Azzalini2DDensity.jl

Azzalini2DDensity

using Azzalini2DDensity, PyPlot

α       = 5*2*(rand(6).-0.5)
f, logf = flexibleDensity2D()

range = -3.0:0.01:3.0

xgrid = zeros(length(range), length(range))
ygrid = zeros(length(range), length(range))
dens  = zeros(length(range), length(range))

for (i,x) in enumerate(range)
    for (j,y) in enumerate(range)
        xgrid[i,j] = x
        ygrid[i,j] = y
        dens[i,j]  = f([x;y], α)
    end
end

cla()
contourf(xgrid, ygrid, dens)