The Maxey-Riley equation (MRE) models the motion of a finite-sized, spherical particle in a fluid. It is a second-order integro-differential equation with a kernel with a singularity at initial time. Because solving the integral term is numerically challenging, it is often neglected despite its often non-negligible impact. Recently, Prasath et al. showed that the MRE can be rewritten as a time-dependent heat equation on a semi-infinite domain with a nonlinear, Robin-type boundary condition. This approach avoids the need to deal with the integral term. They also describe a numerical approach for solving the transformed MRE based on Fokas method. We provide a Python toolbox implementing their approach, verify it against some of their numerical examples and demonstrate its flexibility by computing the trajectory of a particle in a velocity field given by experimental data.