function [y] = ewert2(v);
% ======================================================================
% Unterroutine zum CG-Verfahren fuer quadratische Funktionen. 
%	S.F.
% ----------------------------------------------------------------------
global nfc
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];
D=diag([1 1.0000000001 1.00001 1.0001 2 2.00001 ])     ;
 z1=diag([1 1 1 1 1 1])/Z   ;                 
 A=z1*D*Z;
A=D;

y=A*v;

