% model order

function morder=order(finpre)

global m T

uni=ones(1,m);
mbar=1/(m);
for mod=1:m
ord(mod,m)=mbar;
end
  for i=m+1:T
    for mod=2:m
    ord(mod-1,i)=finpre(mod,i)*ord(mod-1,i-1);
    end
    
    ord(:,i)=ord(:,i)/(uni*ord(:,i));
  end
  
[pmax,morder]=max(ord(:,T));