Calculate pnls increment correctly.
Rename NLMM code to nlmm.jl from poppk.jl
Got the pnls function working, still having problems with fit function