k= function (x) { x^4-8*x^3+16*x+1}
nlm(k,p=2)
nlm(k,p=-1)