answered 13 years ago
I have a feeling that the answer to this ask.sagemath question may also be helpful. It mentions the built-in find_fit method.
find_fit