Describes a method for automated parameter estimation and testing of fit of nonstandard models for mean vectors and covariance matrices, allowing for nonlinear equality and inequality constraints on model parameters. Users need only provide subroutines to evaluate mean vector and covariance matrix according to the model and constraint functions. (SLD)