Program sfhilbert | Compute hilbert transform using different methods. 
[SYNOPSIS]
sfhilbert < in.rsf > out.rsf order=100 ref=1. phase=90. type=
[COMMENTS]
type=t -> time domain convolution
type=f -> frequency domain multiplication and use FFT
type=m -> closed-form design of maximally flat FIR hilbert transformer

[PARAMETERS]
int     | order | =100 |   | 	Hilbert transformer order if type=m 
float   | phase | =90. |   | 	phase shift (in degrees) 
float   | ref | =1. |   | 	Hilbert transformer reference (0.5 < ref <= 1) if type=m 
string  | type | = |   | 	Choosing hilbert transform method, type=t means time domain, 
       type=f means freqency domain, type=m means using more robust algorithm
       the default is type=m 
[DIRECTORY]
user/chenyk
