%[H,k,r,I,h]=fht(h,R;r,K;k,n;I,M,m)
%----------------------------------
%
%Quasi fast Hankel transform of order n.
%
%Input:
% h Function h(r)
% or
% h Signal h(r) ?
% R Maximum radius [m]
% or
% r Radial positions [m] * {default}
% K Maximum frequency [rad/m]
% or
% k Spatial frequencies [rad/m] * {default}
% n Transform order {0}
% or
% I Integration kernel * {default}
% M Maximum points per cycle {7}
% m Minimum points per cycle {5}
%
%Output:
% H Spectrum H(k)
% k Spatial frequencies [rad/m]
% r Radial positions [m]
% I Integration kernel
% h Signal h(r)
%
%