function res=mysinc(x)
%res=sinc(x)
x=x+eps; % avoids errors at x=0

res= sin(pi*x)./(pi*x);



