def szrange(charge, nsingle):
    """
    Make a list giving :math:`S_{z}` values for given charge.

    Parameters
    ----------
    charge : int
        Value of the charge.
    nsingle : int
        Number of single particle states.

    Returns
    -------
    list
        List containing :math:`S_{z}` values for given charge.
    """
    szmax = min(charge, nsingle - charge)
    return list(range(-szmax, szmax + 1, +2))
