Numerical solution of the phase-space dependent backward master equation for the probability distribution of neutron number in a subcritical multiplying sample