function F = myfun(v)
R1 = evalin('base', 'R1');
x = evalin('base', 'x');
% adachstern = evalin('base', 'adachstern');
adachs = evalin('base', 'adachs');
a0 = evalin('base', 'a0');
adachr = evalin('base', 'adachr');
a = evalin('base', 'a');

V = [v^7,v^6,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)); 
if v <= 0.1
    v = 0.1;
    F = 0;
else
    F = [max(es,er)*xRx - 1];
end