Using the MOPED algorithm we determine non-parametrically the Stellar Mass Function of 96,545 galaxies from the Sloan Digital Sky Survey data release one .
By using the reconstructed spectrum due to starlight we can eliminate contamination from either emission lines or AGN components .
Our results give excellent agreement with previous works , but extend their range by more than two decades in mass to 10 ^ { 7.5 } \mathrel { \raise 1.16 pt \hbox { $ < $ } \kern - 7.0 pt \lower 3.06 pt \hbox { { $% \scriptstyle \sim$ } } } M _ { s } / h ^ { -2 } M _ { \odot } \mathrel { \raise 1.16 pt \hbox { $ < $ } \kern% -7.0 pt \lower 3.06 pt \hbox { { $ \scriptstyle \sim$ } } } 10 ^ { 12 } .
We present both a standard Schechter fit and a fit modified to include an extra , high-mass contribution , possibly from cluster cD galaxies .
The Schechter fit parameters are \phi ^ { * } = ( 7.7 \pm 0.8 ) \times 10 ^ { -3 } h ^ { 3 } Mpc ^ { -3 } , M ^ { * } = ( 7.53 \pm 0.04 ) \times 10 ^ { 10 } h ^ { -2 } M _ { \odot } and \alpha = -1.167 \pm 0.004 .
Our sample also yields an estimate for the contribution from baryons in stars to the critical density of \Omega _ { b* } h = ( 2.40 \pm 0.04 ) \times 10 ^ { -3 } , in good agreement with other indicators .
Error bars are statistical and a Salpeter IMF is assumed throughout .
We find no evolution of the mass function in the redshift range 0.05 < z < 0.34 indicating that almost all stars were already formed at z \sim 0.34 with little or no star formation activity since then and that the evolution seen in the luminosity function must be largely due to stellar fading .