Hey
I understand that you wish to fit two equations describing two different parts of the same dataset. I suggest that you break the curve into two parts and fit them separately. Check out this link: Fit two parts of a curve with two equations. Hope that solves your problem. As for the factorial part of the equation, since factorial applies to only integers, I think your approach is right.
Cheers
