function y=F(t)if(t~=0) y=integral(@(x) sin(x)/x,t,1e-10);else y=0;endendx=-1:.01:1;y=x;for i=1:length(x) y(i)=F(x(i));endplot(x,y)