Add solver "lmdif.f"

Levenberg-Marquardt algorithm with finite difference.