function [y] = ewert(v);
% ======================================================================
% Unterroutine zum CG-Verfahren fuer quadratische Funktionen. 
%	S.F.
% ----------------------------------------------------------------------
global nfc
global A
nfc = nfc+1;

%Z=[2 1 0 0 0 0 ;1 2 1 0 0 0;0 1 2 1 0 0;0 0 1 2 1 0; 0 0 0 1 2 1;0 0 0 0 1 2];
%z1=orth(Z)
%D=diag([1 1 1 1 2 2 ])     ;
% z1=diag([1 1 1 1 1 1])/Z   ;                 
% A=z1'*D*z1;
%A=D;

y=A*v;

