nrc=65; xymax=3; freq=1; m=sinx_div_x(xymax,nrc,freq); %figure(3);imagesc(m) %figure(4);surf(m) xv=repmat(-xymax:2*xymax/(nrc-1):xymax,nrc,1); yv=xv'; mexp=zeros(nrc); for r=1:nrc for c=1:nrc Rv=sqrt(xv(r,c)^2+yv(r,c)^2); if Rv>eps mexp(r,c)=sin(Rv*pi*freq)/(Rv*pi*freq); else mexp(r,c)=1; end end % c end %r %figure(1);imagesc(mexp) %figure(2);surf(mexp) assert(~any(any(isnan(m)))) assert(max(max(abs(m-mexp)))<.01)

nrc=127; xymax=3; freq=4; m=sinx_div_x(xymax,nrc,freq); %figure(3);imagesc(m) %figure(4);surf(m) xv=repmat(-xymax:2*xymax/(nrc-1):xymax,nrc,1); yv=xv'; mexp=zeros(nrc); for r=1:nrc for c=1:nrc Rv=sqrt(xv(r,c)^2+yv(r,c)^2); if Rv>eps mexp(r,c)=sin(Rv*pi*freq)/(Rv*pi*freq); else mexp(r,c)=1; end end % c end %r %figure(1);imagesc(mexp) %figure(2);surf(mexp) assert(~any(any(isnan(m)))) assert(max(max(abs(m-mexp)))<.01)

nrc=96; xymax=16; freq=0.5; m=sinx_div_x(xymax,nrc,freq); %figure(3);imagesc(m) %figure(4);surf(m) xv=repmat(-xymax:2*xymax/(nrc-1):xymax,nrc,1); yv=xv'; mexp=zeros(nrc); for r=1:nrc for c=1:nrc Rv=sqrt(xv(r,c)^2+yv(r,c)^2); if Rv>eps mexp(r,c)=sin(Rv*pi*freq)/(Rv*pi*freq); else mexp(r,c)=1; end end % c end %r %figure(1);imagesc(mexp) %figure(2);surf(mexp) assert(~any(any(isnan(m)))) assert(max(max(abs(m-mexp)))<.01)