function [W] = taylor(xX,yY)

pochodne = zeros(length(xX))
pochodneD = zeros(4)
pochodneT = zeros(3)

    for i = length(xX) : -1: 2
            pochodne( i) = yY(i) - yY(i-1) ;
            disp(pochodne)
    end

    for i = 5 : -1: 3
            pochodneD( i - 1) = pochodne(i) - pochodne(i-1) ;
            disp(pochodneD)
    end

     for i = 4 : -1: 3
            pochodneT( i - 1) = pochodneD(i) - pochodneD(i-1) ;
            disp(pochodneT)
    end

    pochodneC = pochodneT(3) - pochodneT(2);
     disp(pochodneC)

     W = 1/(xX(2)/xX(1)) * (pochodne(5) + 1.0/2.2 * pochodneD(4) + 1.0/3.0  * pochodneT(3) + 1.0/4.0 * pochodneC);