Monte Carlo simulation for jet fragmentation in SUSY QCD