function y=sinc(t)
y=ones(size(t));
nn=find(t~=0);
y(nn)=sin(pi*t(nn))./(pi*t(nn));