function y=gu(u)
% partial derivative of g with respect to u
y=4*u^3-1;