function y = sinc_squared(x)
// The square of the sinc function, with special case for x = 0.
// We use the Scilab library's version of the sinc(x) function.
//
// MWR 4/9/2007

  if (x == 0)
    y = 1;
  else
    y = sinc(x)*sinc(x);
  end

endfunction
