gravity code snippets in matlab