Python codes that numerically compute inflaton trajectories for ultra slow-roll inflation to evaluate the power spectrum of curvature pertrubations by solving the Mukhanov-Sasaki equation.