function y=hardy(a,b,f)
%
%   function y=hardy(a,b,f)
%
% Esta función aplica el método de Hardy para poporcionar una aproximación 
% de la integral de la función f en el intervalo [a,b].
%
% La función f debe programarse vectorialmente de forma que si x = (x_i) es
% un vector, entonces f(x) es el vector "columna" cuyas componentes son los
% valores f(x_i).
%
x=linspace(a,b,7)';
w=[41 216 27 272 27 216 41];
fx=feval(f,x);
y=(b-a)/840*(w*fx);