function y=f4_newton(x)

y = tan(x/4)-1/tan(x/4);
