%-----------------------------------------
% Optim des Kräftegleichgewichtes
%-----------------------------------------
clc;

    x0=[-1,0.1] ;  %Startwert
    OPTIONS=optimset('Display','iter', 'TolX', 0, 'TolCon', 0, 'LargeScale','off','Diagnostics','on'...
                      ,'DiffMinChange', 0.01 ,'DiffMaxChange', 0.02 , 'MaxFunEvals',300);

    [X,FVAL,EXITFLAG,OUTPUT] = fsolve(@optim_fkt,x0,OPTIONS);

    OUT_alpha_s(i)=(X(2)*180/pi+50);
    OUT_spiel(i)=s;
    OUT_phi(i)=X(1)*180/pi;