function [ ke ] = Crelib( r )
%
%  This function computes the reliability factor
%
%  r = reliability in percent

if (r < 50)
    ke = 1;
elseif (r > 99.9999)
    ke = 0.620;
else
    p = [50 90 95 99 99.9 99.99 99.999 99.9999]';
    k = [1.000 0.897 0.868 0.814 0.753 0.702 0.659 0.620]';

    ke = interp1(p, k, r);
end

end

