function fnonlp = fnonlp(u)
fnonlp = .0000001*(-4)*u^3;