In the papers attached, they considered the acceleration as a quadratic expression of time with some unknown coefficients. The trajectory of the lander can be obtained by integrating the acceleration twice. Also it is required to update the time-to-go and all coefficients in all guidance call. But using MatLab programming, I am facing some difficulties for the simulation.