


DOUSTAFOD Discrete Oustaloup approximation of a s^r operator
USAGE: Z = DOUSTAFOD(R,N,WB,WH,TS)
where R is the operator order,
N relates to the number of zeros and poles: N_ZP = 2N+1,
WB and WH are lower and higher frequency bounds of approximation,
TS is the sampling time.