add phase_contrast_simu

this script is to simulate first order differential phase contrast
projection, in other word, gradient-hibert- transform backprojection.