This module provides the fit() function for maximum likelihood estimation
and the femtofit class for representing fitted models. The interface
follows standard R conventions, implementing base R generics like
coef(), vcov(), confint(), logLik(), etc.