Mean-field treatment of nonlinear susceptibilities for a ferromagnet of arbitrary spin