SUB acceleration(position(),accel(),GM)
    LET r2 = position(1)*position(1) + position(2)*position(2)
    LET r3 = r2*sqr(r2)
    FOR i = 1 to 2
        LET accel(i) = -GM*position(i)/r3
    NEXT i
END SUB