Based on the neutron induced fission fragment mass distribution data up to neutron energy 20 MeV measured with the double kinetic energy method (KEM) and the radio active method (RAM), the systematics of fission fragment mass distribution was investigated by using 5 Gaussian model and the systematics parameters were obtained by fitting the experimental data. With the systematics, the yields of any mass A and at any energy in the region from 0 to 20 MeV of neutron energy can be calculated. The calculated results could well reproduce the experimental data measured with KEM, but show some systematical deviation from the data measured by RAM, which reflects some systematical deviations between the two kinds of measured data.The error of systematics yield was calculated in an exact error transformation way, including from the error of the experimental yield data to the error of the discrete parameters, then to the systematics parameters,and at last to the yield calculated with systematics.