function F = YY(v)
init
optimierungsergebnis
R1 = [    9.1882   28.4887   23.2505    8.1342    1.5389;
   28.4887  123.5995  110.7081   39.7564    7.8915;
   23.2505  110.7081  177.5209   67.8702   16.7624;
    8.1342   39.7564   67.8702   28.9017    7.0587;
    1.5389    7.8915   16.7624    7.0587    2.2227]
x = [-0.0392;   0.0085 ;  -0.0358   ; 0.0167    ;0.2226]

V = [v^5,v^4,v^3,v^2,v]
D = diag(V)

R=D\R1/D
xRx = x'*R*x

ks = D\adachstern-a0
hs = D\adachs-a0
hr = D\adachr-a

es = hs'/R*hs
er = hr'/R*hr

hsx = sqrt(sum((hs'*x).^2))
% F = [hsx-1; max(es,er)*xRx-1];
F = [max(es,er)*xRx-1];