step = 0.2
v = []
t = var("t")
for s in srange(0.01, 4.0*pi, step):
cycloid = parametric_plot((t - sin(t), 1 - cos(t)), (t, 0, s))
C = circle((s, 1), 1)
Cp = (s - sin(s), 1 - cos(s))
C_point = point(Cp, pointsize=25, rgbcolor="brown")
C_line = line([(s, 1), Cp], thickness=2, rgbcolor="brown")
v.append(cycloid + C + C_point + C_line)
animation = animate(
v, figsize=[10, 2], aspect_ratio=1,
xmin=-0.5, xmax=14, ymin=0, ymax=2.2
)
show(animation, delay=10)