integration by rectangles

I want to do the integral by rectangles of this function but I get an error, someone can help me please.

funcion:
function A = f_intRect(f,a,b,n)
% Devuelve la aproximacion de la integral por rectangulos centrados
% A: valor de la integral numerica
% f: funci ́on a evaluar
% a: extremo izquierdo
% b: extremo derecho
% n: cantidad de divisiones (numero de rectantulos)
h=(b-a)/n; %Base
k = 0:n-1; % ́ındices
centrosX = a+h/2+k*h; %Posici ́on de los centros en X
A = h*sum(feval(f,centrosX)); 
%Area = SUM base x altura
end

script
%función a evaluar f(x) = sin(100x/pi) + 9x^2 

A = f_intRect(f,a,b,n);
f=sin*(100*x/pi)+9*x^2;
a=-0.5;
b=1;
n=2;

disp(A)