We report on the calculation of the inclusive semileptonic decay of the $D_s$ meson on the lattice. We simulate the $D_s \rightarrow X_s\ell\nu_\ell$ process with M\"obius domain-wall charm and strange quarks, whose masses are approximately tuned to their physical values. Our simulations cover the whole kinematical region. The focus of this work is to present updates on our strategies towards estimating the systematic uncertainties in the determination of the inclusive decay rate. We specifically focus on the systematic errors due to the choice of our approximation strategy and finite-volume effects.