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