function [y]=sinc(x)
%function [y]=sinc(x)
%x may be a vector
pigreco=3.141592653589793;
y=zeros(1,length(x));
y(x==0)=1;
y(x~=0)=sin(pigreco*x(x~=0))./(pigreco*x(x~=0));