function [y,x] = obszh_form(a,b,n,K,f)
y=1:n;
x=linspace(a,b,n);
h = x(2)-x(1);
y(1)=f(x(1));
A=1:n;
A(1)=h/2;
A(2)=h/2;

for i=2:n
f=f(x(i));
for j=1:(i-1)
if mod(i,2)==1
A = h.*[1, 2*(1+mod((1:i-2),2)), 1]./3;
%%A(i,i)=h/3;
else
A = h.*[1/3, 2/3*(1+mod((1:i-3),2)), 5/6, 1/2];
%%A(i)=h/2;
end
f=f+A(i,j)*K(x(i),x(j))*y(j);
end
y(i)=f+A(i,i)*K(x(i),x(j));
end
end