function p(n, r, p_0, T)
    A = fill(T(0.0), n + 1)
    p = T(p_0)
    A[1] = p
    for i = 1 : n
        p = p + r * p * (T(1.0) - p)
        A[i + 1] = p
    end
    return p, A
end