%% f(y) for sorG.
function [s] = sorf(y)
%s = 1.0;
s = 1.0 + sinc(y);
