function fnonl = fnonl(u)
fnonl = .00000001*(300^4 - u^4);