Revert to previous version of pln for stability. Change the method in pln_solver() to BFGS for better results.