clear
clear all
format longE
d = 10^(-3);

x = linspace(2-d, 2+d, 1000);

value1 = 0; %wartosci dla Algorytmu1 (x - 2)^4
value2 = 0; %wartosci dla Algorytmu2 (taki sam wzor, tylko po rozlozeniu na czynniki pierwsze)
                 % x*x*x*x - 4*x*x*x*2 + 6*x*x*4 - 4*x*8 + 16

diff = 0;
h=d/500;

for k=1:1000
   value1(k) = fun1(x(k));
   value2(k) = fun2(x(k));

   diff(k) = abs(value1(k)-value2(k));

   xk=2-d+k*h;
   max_k(k) = abs(fun1(xk)-fun2(xk));
end
plot(x,value1, x, value2, x, diff);
max(max_k)
%max(diff)
%    xk=2-d+k*h;
%    max_k = abs(fun1(xk)-fun2(xk));
%max(max_k)


function A1=fun1(x)
    a = (x - 2);
    A1 = a^4;
end

function A2=fun2(x)
    A2 = x*x*x*x-4*x*x*x*2+6*x*x*4-4*x*8+16;
end