Add matlab example code for orbital motion