function prior(w0::Float64, w1::Float64)::Float64
    μ = [0., 0.]
    return pdf(MultivariateNormal(μ, I(length(μ))), [w0, w1])
end