function x = rect(f, a, b, n)
	x = (b-a)/n * sum(f(a+(b-a)/(2*n):(b-a)/n:b-(b-a)/(2*n)));