N = 10;

A = zeros(N,N);

for i = 1:N;
	for j = 1:N;
    	if i == j
				A(i,j) = 1.5*cos(i+1.2*j);
        else
				A(i,j) = sin(0.5*i+0.2*j);
        end
    end
end

lum = 3.38304;

B = A - lum*eye(10);

x = [-0.43674
-0.906357
-1.96334
-1.08286
-1.2693
-1.07245
0.362511
1.68291
2.11274
1];
% 
% Q = hess(A);

% N = 5;
% 
% A = zeros(N,N);
% 
% for i = 1:N;
% 	for j = 1:N;
%     	A(i,j) = i+j;        
%     end
% end

%  B = A;
% 
% % B = hess(A);
% % Q = hess(A);
% ee = 10^-12;
% 
% for r = 1:N-2
%     temp = 0;
%     flag = 0;
%     for i = r+1:N
%         if abs(A(i,r)) > ee
%             flag = 1;
%             break;
%         end
%     end
%     
%     if flag
%         for i = (r+1):N
%             temp = temp + A(i,r)*A(i,r);
%         end
%          dr = (temp)^0.5;
% 
%         if A(r+1,r) <=0
%             cr = dr;
%         else
%             cr = -dr;
%         end
% 
%         hr = cr*cr - cr*A(r+1,r);    
%         ur = zeros(N,1);
% 
%         for i = 1:N
%             if i <= r
%                 ur(i) = 0;
%             end
%             if i == r+1
%                 ur(i) = A(i,r) - cr;
%             end
%               if i > r+1
%                 ur(i) = A(i,r);
%             end
%         end
%         
%         pr = A'*ur/hr;
%         qr = A*ur/hr;
%         tr = pr'*ur/hr;
%         wr = qr - tr*ur;
%         A = A-wr*ur'-ur*pr';
%         
%     end
%         
%  end
% 
% 
% % hr= 113.82085548649;
% % 
% % u=[0;12.2736;4;5;6];
% % 
% % pr = A'*u/hr;
% % 
% % qr = A*u/hr;
% % 
% % wr = qr - tr*u;
% % 
%