clear all,clc
f=inline('1./sqrt(x+0.0001)');
a=0;
b=1;
P=0;
T=0;
S=0;

for n=3:10
    c=2^n;
    hp=(b-a)/(2*c);
    m=a:hp:b;

    for i=1:c
    y=f(m(2*i));
    P=P+2*hp*y;
    end
bP=abs((((sqrt(10001)-1)/50)-P)/((sqrt(10001)-1)/50));
loglog(hp,bP,'*');
hold on;
P=0;

for i=1:c
    y=f(m(2*i+1));
    T=T+hp*(f(m(1))+f(m(2*c+1)))+2*hp*y-2*hp*f(m(2*c+1));
end

T
bT=abs((((sqrt(10001)-1)/50)-T)/((sqrt(10001)-1)/50));
loglog(hp,bT,'r*');
hold on;
T=0;

for i=1:c
    w1=f(m(1))+f(m(2*c+1));
    w2=f(m(2*i));
    w3=f(m(2*i+1))-f(m(2*c+1));
    S=S+(hp/3)*(w1+4*w2+2*w3);
end
S
bS=abs((((sqrt(10001)-1)/50)-S)/((sqrt(10001)-1)/50));
loglog(hp,bS,'g*');
xlabel('ilośc podprzedzialow'),ylabel('blad wzgledny'),title('Funkcja 2 ');
legend('prostokatow','trapezow','Simpsona',4);
hold on;
S=0;
end;