open Impact
clear()
t <- 0.
p <- 1.0
m <- 80.0
while t < 200. do 
    dot (t, p) blue   
    t <- t + 1.
    s <- 0.001*p*(m-p)
    p <- p + s