Fitting the sunspot cycles 10–21 by a modified F-distribution density function