function [s]=prostokat(a,b,f)
    s=0;
    st=0;
    n=5;
    xd = (b-a)/n;
    X=[a:xd:b];
    for i = a:length(X)-1;
        x1=a+i*xd;
        st=st+f(x1-xd/2);
        if(i<n)
            s=s+f(x1);
        end
    end
    s=(xd/6)*(f(a)+f(b)+2*s+4*st);
    return;
end

% function [s]=prostokat(a,b,f)
%     s=0;
%     n=1000;
%     xd = (b-a)/n;
%     X=[a:xd:b];
%     for i = a:length(X)-1;
%         s=s+f(a+i*xd);
%
%     end
%     s=(s+(f(a)+f(b))/2)*xd;
%     return;
% end

% function s=prostokat(a,b,f)
%     s =0;
%     n=10000;
%     xd = (b-a)/n;
%     i=1;
%     while(i < n)
%         s=s+f(a+i*xd);
%         i=i+1;
%     end
%     s=(s+(f(a)+f(b))/2)*xd;
%     return;
% end

% function s=prostokat(a,b,f)
%     s =0;
%     n=10000;
%     xd = (b-a)/n;
%     i=1;
%     while(i <= n)
%         s=s+f(a+i*xd);
%         i=i+1;
%     end
%     s=s*xd;
%     return;
% end