BayesHistogram.jl main procedure: bayesianblocks( t::AbstractVector{T}; p0::T = T(5) / T(100), resolution::T = T(Inf), mincounts::Int = ceil(Int64, sqrt(length(t))/2) )