% correl(x,y,N)

function xm = dosa(x,fa,M)

xf = iirtp(fa,0.45*fa/M,1,8,x);

xm = xf(1:M:lge(x));