A Kepleres feladathoz kapcsolódva az alábbi procedura egy bolygó Nap körüli mozgását szimulálja Newton-törvény és a gravitációs törvény alapján.
Az elméleti háttere Feynman: Mai fizika I. kötetében le van írva, és onnan származnak a példa kezdeti paraméterei értékei is. (Lehetne szebben is megírni, de most annyi idõt nem töltöttem vele.)
restart;with(plots):
bolygo:=proc(N,dt,x0,y0,vx0,vy0,M,G)
description "Bolygómozgás számítása Gravitációs és Newton törvénybõl Forrás: Feynmann : Mai fizika I. ";