Hmm...
Pomijając sprawy R minimalnego i R maksymalnego, oraz korzystając z programowania biegunowego...
zakładając, że Δω = 360/n, gdzie n - dowolna wartość uzależniona od dokładności
ΔR = ae/Δω
ω0 - kąt początkowy
CC x y
L PR Rmin AN ω0
LBL
L PR Rmin+ΔR AN ω0+Δω
CALL LBL dopuki Rmin+ΔR<Rmax
na Rmax wykonać jeszcze jedno pełne przejście

Oczywiście obliczenia trzeba wykonać na Q-parametrach wewnątrz pętli, a do zapętlenia użyć instrukcji skoku warunkowego.

Nie kumam, o co kaman... w (okrągłym) otworze trochoida staje się spiralą o skoku równym ae.