function y = ftest(x)
    % Test function for Simpson's rule
    y = 1./(x.^3-2*x-5);