function y=df4_newton(x)

y = 1/4*1/((cos(x/4))^2) + 1/4*1/((sin(x/4))^2);
