function [ value, time ] = integrate ( n )

  value = 0.0;

  tic;
  x = 1.0;
  for i = 1 : n
    fx = bessely ( 4.5, x );
    value = value + fx / n;
    x = x + 1.0 / ( n - 1);
  end
  time = toc;

  return
end
