Implementation of sampling formulas for the unified neutral