function r= par_soed(varargin)
% 琼圜屙桢 1/(1/x1+...+1/xn)
    r= 1./[varargin{:}];
    r= 1/sum(r);
end