function [] = testgrad(fktname,x0,gradname)
% ======================================================================
% Funktion, dient zum Testen der Ableitung an der Stelle x0
% es werden folgende Angaben ben"otigt:
%
% fktname  : Name der Funktion;
% x0       : Startwert;
% gradname : Name des Gradienten;
% ======================================================================

[F0]=feval(fktname,x0);		%	f(x0)
[g]=feval(gradname,x0);		%	g(x0)

%	Numerische Approximation von g(x0)
h=1E-8;
[F]=feval(fktname ,x0*ones(size(x0))'+h*eye(length(x0)));%  "f(x0+h)"
dF=(F-F0)/h;		%	g~(x0)

%	Bestimmung des maximalen Fehlers
[w,i]=max(abs(g-dF'));

if w>sqrt(h) 
   error(['Gradient "uberpr"ufen! Tip: ',int2str(i),'. Element falsch.'])
end;

return;






