function foo=f(x)
	foo = 5 * x+(sin(9*x-5)/x);
endfunction

steps = 10;
x_start = 1;
x_end = 10;
x_delta = (x_end - x_start) / steps;

total = 0; 
for i=0:(steps-1),
	x = x_start + i * x_delta;
	total=total+f(x);
end
total = total * x_delta;
printf("Rectangular integration value %f\n", total);