put "Plotting a trochoid. The fast way."

const pi: real := 3.141593
var x,y: real
var xn,yn: int
var u, r, a, du: real

setscreen("nocursor")
r  := 100
a  := 25
du := 0.01
u  := 0

x  := 0
y  := r

loop
	exit when u > 10*pi
	xn := round(x)
	yn := round(y)
	drawdot(xn,150+yn,1)
	u  += du
	x  := x + du*(a - y)
	y  := y + du*(x - a*u)
end loop