% Orr-Sommerfeld Re_c(k) function
   function [Rec,kc] = oserecrit(k,Re1,Re2,M)
   Rec = fzero('osefun2',[Re1;Re2],[],k,M);
   kc = k ;

