  The mass function of galaxy clusters is a powerful tool to constrain
cosmological parameters, e.g., the mass fluctuation on the scale of 8 h^{-1}
Mpc, sigma_8, and the abundance of total matter, Omega_m. We first determine
the scaling relations between cluster mass and cluster richness, summed r-band
luminosity and the global galaxy number within a cluster radius. These
relations are then used to two complete volume-limited rich cluster samples
which we obtained from the Sloan Digital Sky Survey (SDSS). We estimate the
masses of these clusters and determine the cluster mass function. Fitting the
data with a theoretical expression, we get the cosmological parameter
constraints in the form of sigma_8(Omega_m/0.3)^{alpha}=beta and find out the
parameters of alpha=0.40-0.50 and beta=0.8-0.9, so that sigma_8=0.8-0.9 if
Omega_m=0.3. Our sigma_8 value is slightly higher than recent estimates from
the mass function of X-ray clusters and the Wilkinson Microwave Anisotropy
Probe (WMAP) data, but consistent with the weak lensing statistics.
