How would I program something like this in postscript?

a-radiusOfBigCircle . b - radiusOfSmallCircle d-lengthOfArm . r- numberOfTimesAround The hypotrochoid is produced by circling for t-0..360"r degrees using the following formulas. x(t)-(Ca -b)* cos(t)) (d *cos (Ca-b) * t)/b )) y(t)-(Ca -b) sin(t)) - (d * sin( (Ca-b) * t)/b )) For example Chypotrochoid 100 25 40 1) would produce a-radiusOfBigCircle . b - radiusOfSmallCircle d-lengthOfArm . r- numberOfTimesAround The hypotrochoid is produced by circling for t-0..360"r degrees using the following formulas. x(t)-(Ca -b)* cos(t)) (d *cos (Ca-b) * t)/b )) y(t)-(Ca -b) sin(t)) - (d * sin( (Ca-b) * t)/b )) For example Chypotrochoid 100 25 40 1) would produce