function NTF=clans6(order,OSR,Q,rmax,opt)%Version of clans for MATLAB >=6% Create the initial guessNTF = synthesizeNTF(order,OSR,opt,1+Q,0);Hz = NTF.z{1};x = zeros(1,order);odd = rem(order,2);if odd    z = NTF.p{1}(1)/rmax;    if any(abs(z))>1 %project poles outside rmax onto the circle	z = z./abs(z);    end    s = (z-1)./(z+1);    x(1)= sqrt(-s);endfor i=odd+1:2:order    z = NTF.p{1}(i:i+1)/rmax;    if any(abs(z))>1 %project poles outside rmax onto the circle	z = z./abs(z);    end    s = (z-1)./(z+1);    coeffs=poly(s);    wn = sqrt(coeffs(3));    zeta = coeffs(2)/(2*wn);    x(i) = sqrt(zeta);    x(i+1) = sqrt(wn);end% Run the optimizeroptions = optimset('TolX',1e-6, 'TolFun',1e-6, 'TolCon',1e-6, 'MaxIter',1000 );options = optimset(options,'Display','off');options = optimset(options,'Diagnostics','off');options = optimset(options,'LargeScale','off');[x,fval,exitflag,output] = fmincon(@(x)dsclansObj6a(x,order,OSR,Q,rmax,Hz),x, ...[],[],[],[],[],[],@(x)dsclansObj6b(x,order,OSR,Q,rmax,Hz), options );NTF = dsclansNTF(x,order,rmax,Hz);returnfunction f=dsclansObj6a(x,order,OSR,Q,rmax,Hz)% Objective function for clans; Optimization Toolbox version >= 6% f is the magnitude of H at the band-edgeH = dsclansNTF(x,order,rmax,Hz);f = abs(evalTF(H,exp(1i*pi/OSR)));returnfunction [g,g_eq]=dsclansObj6b(x,order,OSR,Q,rmax,Hz)% Constraint function for clans; Optimization Toolbox version >= 6% g =||h||_1 - QH = dsclansNTF(x,order,rmax,Hz);g = sum(abs(impulse(H,100))) -1 - Q;g_eq = [];return