Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm