clear all
clc

%Wartosci poczatkowe wspolrzendych przedzialu
a = 20;
b = 30;

%Wartosci aktualne wspolrzendych przedzialu
an = a;
bn = b;

%Wartosc srodka przedzialu
tn = (an + bn) /2;

%Wartosc dlugosci przedzialu
dl = bn - an;

%Wartosci chwilowe funkcje od wartosci przedzialow oraz srodka przedzialu
yan = (an.^6) + (5*an.^3) +2;
ybn = (bn.^6) + (5*bn.^3) +2;
ytn = (tn.^6) + (5*tn.^3) +2;

%Pomocniczy iterator
n = 0;

while((yan < ytn) || (ybn < ytn) )
    if (yan > ybn)
        an = an + 0.5 * dl;
        bn = bn + 0.5 * dl;


    else
        an = an - 0.5 * dl;
        bn = bn - 0.5 * dl;

    end
    tn = (an + bn) /2;
    n = n +1;
    yan = (an.^6) + (5*an.^3) +2;
    ybn = (bn.^6) + (5*bn.^3) +2;
    ytn = (tn.^6) + (5*tn.^3) +2;
end

%Wyznaczenie wartości X minimum
X_coordin_min = fminbnd(@insertX,an,bn);
%Wyznaczenie wartości Y minimum
Y_coordin_min = (X_coordin_min.^6) + (5*X_coordin_min.^3) +2;