function alpha=linesearch_secant(grad,x,d)
%Line search using secant method
%Note: I'm not checking for alpha > 0.

epsilon=10^(-5); %line search tolerance
max = 200; %maximum number of iterations
alpha_curr=0;
alpha=10^(-5);
dphi_zero=feval(grad,x)'*d;
dphi_curr=dphi_zero;
 
i=0;
while abs(dphi_curr)>epsilon*abs(dphi_zero),
  alpha_old=alpha_curr;
  alpha_curr=alpha;
  dphi_old=dphi_curr;
  dphi_curr=feval(grad,x+alpha_curr*d)'*d;
  alpha=(dphi_curr*alpha_old-dphi_old*alpha_curr)/(dphi_curr-dphi_old);
  i=i+1;
  if (i >= max) & (abs(dphi_curr)>epsilon*abs(dphi_zero)),
    disp('Line search terminating with number of iterations:');
    disp(i);
    break;
  end  
end %while