h = 0.3;
x = 0.5;
p1 = '2*exp(x^2+1)*x';
p2 = '2*exp(x^2+1)*(2*x^2+1)';
p3 = '4*exp(x^2+1)*x*(2*x^2+3)';
w1 = eval(p1);
w2 = eval(p2);
w3 = eval(p3);

for i=1:70
    H(i) = h*2^(-(i-1));
    RoznicaCentralna(i) = rcentralna(x,H(i));
    RoznicaProgresywna(i) = rprogresywna(x, H(i));
    BladWzglednyCentralna(i) = (RoznicaCentralna(i)-w1)/w1;
    BladWzglednyProgresywna(i) = (RoznicaProgresywna(i)-w1)/w1;
    BladCentralna(i) = ((-1/6)*(H(i))^2*w3)/w1;
    BladProgresywna(i) = ((1/2)*H(i)*w2)/w1;
end

loglog(H, abs(BladWzglednyCentralna), H, abs(BladWzglednyProgresywna), H, abs(BladCentralna), H, abs(BladProgresywna))
xlabel('log(h)');
ylabel('log(B)');