function f=edo_f(x)
f = 4*(pi^2+1)*x.*sin(2*pi*x) - 4*pi*cos(2*pi*x);
return