Implement better fitting algorithms from Optim package. May also want to look into NLopt and JuMP libraries