GradConjugado

function [Xotimo CONT]=GradConjugado(alfa, x, f, grad, Hessiana, d0, x0)

CONT=1;

Xk=x0;
gradiK=subs(grad, x, x0);

dk=d0;
while vpa(norm(gradiK))>1e-4
    [alfaK Aux]=newton(alfa, x, f, dk, Xk);
    Xk=Xk+(alfaK*dk);

    Bk=((gradiK).' *(subs(Hessiana,x,Xk))*(dk))/((dk).' *(subs(Hessiana,x,Xk))*(dk))
    dk=-gradiK+Bk*dk
    CONT=CONT+1;
        if(CONT>200)
        break;
    end
end
Xotimo=vpa(Xk);