saf.m

来自「还是一些m文件集合」· M 代码 · 共 7 行

function y=saf(x)
n=length(x);
for k=1:n
	x(k)=x(k)+(x(k)==0)*eps;
	y(k)=sin(x(k))/x(k);
end