function y=df1_newton(x)

y = 1/x + 3*(x+1)^2;
