  New formulas are given for the grand partition function of paraboson systems
of order p with n orbitals and parafermion systems of order p with m orbitals.
These formulas allow the computation of statistical and thermodynamic functions
for such systems. We analyze and discuss the average number of particles on an
orbital, and the average number of particles in the system. For some special
cases (identical orbital energies, or equidistant orbital energies) we can
simplify the grand partition functions and describe thermodynamic properties in
more detail. Some specific properties are also illustrated in plots of
thermodynamic functions.
