sinc.m

来自「A set of C++ and Matlab routines impleme」· M 代码 · 共 8 行

function y = sinc(x)

y = pi * x;

sel = abs(y) > 1e-15;

y(sel) = sin(y(sel)) ./ y(sel);
y(~sel) = 1;