function res = vepk(eps,v)

le = length(eps);

mat = [];

for k = 1:le
    mat(1,k) = v(k);
end;

for k = 1:le
    mat(2,k) = (2-eps(k))*v(k);
end;

for ke = 1:le
    for kp = 1:(le-2)
        mat(kp+2,ke) = mat(kp+1,ke)*(2-eps(ke))-mat(kp,ke);
    end;
end;

res = mat;

end
