function p = Estep(X, compParams)
  pix = mixturePiX(X, compParams) ;
  p = op(pix,'./',sum(pix,2)) ;
  
    