function H = hurst_rd(x)
%
% compute hurst parameter with relative dispersion method
%
H=hurst1(x);