%Gewichte als Parameter - globale Variable
global w;
w=1;

%Federkonstante als Parameter - globale Variable
global K;
K=[1;1];

% Anfangswert
x=[0.4;0];

% maximale Iterationsanzahl
maxit = 10000;

% Schrittweite 
eta=1;

% typische Werte
[f,g]=obj3a(x);
typf = f;
typx =[1;1];

j = 1;

%tolerance
tol =1e-4;
relgrad=tol+1;

while ((relgrad > tol) && (j < maxit))

 %Iterationsschritt
 x = x - eta*g

 % Auswertung von Funktion und Gradient
 [f,g]=obj3a(x);

 % Berechnung des relativen Gradienten
 relgrad = 0;
 for k=1:length(g)
  relgrad=max([relgrad; (abs(g(k))*max([x(k);typx(k)])/(max([f;typf])) ]);
 end;
 
 % Ausgabe des absoluten und relativen Gradienten
 [norm(g);relgrad]

j=j+1;

end;