MATLAB code describing the operation of the Wankel engine.

I would need the code to plot the geometry of the Wankel engine:

- the 2-lobed epithrochoid

- the rotor (Reauleaux's triangle)

- the eccentric shaft

(if possible also through an animation)