GADfit
Global nonlinear optimization with automatic differentiation
...The model curves (fitting functions) can be of essentially arbitrary complexity. This includes any nonlinear combination of elementary and special functions, single and/or double integrals, and any control flow statement allowed by the programming language. See the latest user guide under Files.