library("mcPAFit")
T <- 10000
t <- 1:T
p_t <- function(p_max,p_min,t) {p_max + (t-1)/T * (p_min - p_max)}
#
plot(t, p_t(0.5,0.01,t))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.