function [y] = trapezowa(a,b,n,f)
h = (b - a)/n;
y = (feval(f,a) + feval(f,b))/2;
for i = 1 : n-1
y = y + feval(f,a+i*h);
end
y = h*y;

endfunction 




function [y] = funkcja1(x)
y=sin(4*x)/x; // napiszac
endfunction 



clc
clear all
disp('Zad. 4')
y=trapezowa(2,10,12,'funkcja1')