astrophysics code snippets in matlab

simulate the poynting–robertson effect in matlab