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