y=stud_arofmvc1(x)

function y=stud_arofmvc1(x)
alfa=x(1);
nu=x(2);
g1 = 2.30753;
g2 = 0.27061;
g3 = 0.99229;
g4 = 0.04481;
p = alfa;
k = nu;
p0 = 0.5 - sign(0.5 - p) * (0.5 - p);
z1 = sqrt(-2 * log(p0));
z = z1 - (g1 + g2 * z1) / (1 + z1 * (g3 + g4 * z1));
z = z * sign(0.5 - p);
t = sign(z) * sqrt(k * (exp((z * (8 * k + 3) / (8 * k + 1)) ^ 2 / k) - 1));
if abs(t) >= 3
    t = t / (1 + 0.05 * log(abs(t)) ^ 2 / (2 * k - 1));
end
y=t;