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