I do realize that the normal distribution is defined by mean and variance, so this particular example is rather futile. I was deliberately trying to use a trivial example to help me figure out the calling conventions, etc. throughout the package. It's _possible_ that fit() is broken, but more likely that my understanding is.