Add solver "lmdif.f"

Add "lmdif.f",
  which is a Least-Squares solver
  based on Levenberg-Marquardt Algorithm.