  We present results from a new Monte Carlo simulation for jet fragmentation in
QCD and SUSY QCD for large primary energies $\sqrt s$ up to $10^{16}$ GeV. In
the case of SUSY QCD the simulation takes into account not only gluons and
quarks as cascading particles, but also their supersymmetric partners. A new
model-independent hadronization scheme is developed, in which the hadronization
functions are found from LEP data. An interesting feature of SUSY QCD is the
prediction of a sizeable flux of the lightest supersymmetric particles (LSPs),
if R-parity is conserved. About 10% of the jet energy is transferred to LSPs
which, owing to their harder spectra, constitute an important part of the
spectra for large $x=E/E_{jet}$. Spectra of protons and of secondary particles,
photons and neutrinos, are also calculated. These results have implications for
the decay of superheavy particles with masses up to the GUT scale, which have
been suggested as a source of ultrahigh energy cosmic rays.
