function y = f(x)

y = log(x.^2+1)-exp(0.4*x).*cos(%pi*x)