clear all, clc
x=2;
h0=.2;
f=@(x)  x^2*cos(x);
p1=@(x) 2*x*cos(x)-x^2*sin(x);
p2=@(x) 2*cos(x)-4*x*sin(x)-x^2*cos(x);
p3=@(x) -6*sin(x)-6*x*cos(x)+x^2*sin(x);

for i=1:40

    h(i)=h0*2^(-(i-1)); %%krok

    dp(i)   = (f(x+h(i))-f(x))/h(i);    %%progresywna
    wzdp(i) = abs((dp(i)-p1(x))/p1(x)); %%blad względny roznicy prog.
    bzdp(i)  = (h(i)/2*p2(x))/p1(x);     %%błąd bezwzględny

    dc(i)   = (f(x+h(i))-(f(x-h(i))))/(2*h(i))    %%różnica centralna
    wzdc(i)  = abs((dc(i)-p1(x))/p1(x));         %%błąd względny różnicy centralnej
    bzdc(i)  = (((h(i)^2)*p3(x))/6)/(p1(x));            %%błąd bezwzgędny
end

loglog(h,wzdp,'r.');
hold on
grid on
loglog(h,wzdc,'blue.');
loglog(h,abs(bzdp),'black');
loglog(h,abs(bzdc));
xlabel('log(h)');
ylabel('log(błąd)');
legend('różnica progresywna','różnica centralna','błąd względny odcięcia r. prog','błąd względny odcięcia r.cent.');
legned('location','South');
hold off