L=cellfun(@(x) x.*0,coeff,'UniformOutput',fal) S=cellfun(@(x) x(find(x>r))==1,coeff,'UniformOutput',false)