Index: tol/btol/bmodel/estim.cpp
===================================================================
--- tol/btol/bmodel/estim.cpp	(revision bab7ea70f5634f82d61a706f3493f50dff6560d7)
+++ tol/btol/bmodel/estim.cpp	(revision eccb6ceaa616e785eb426015152696e79dca7d63)
@@ -2498,8 +2498,26 @@
 
 //--------------------------------------------------------------------
+DeclareContensClass(BMat, BMatTemporary, BMatARMA_ACVF_ARFilter);
+DefExtOpr(1, BMatARMA_ACVF_ARFilter, "ARMA_ACVF_ARFilter", 
+  2, 2, "Matrix Polyn",
+  "(Matrix ACVF, Polyn phi)",
+  "",
+  BOperClassify::Sthocastic_);
+//--------------------------------------------------------------------
+void BMatARMA_ACVF_ARFilter::CalcContens()
+//--------------------------------------------------------------------
+{
+  BMatrix<BDat>& acvf = Mat(Arg(1));
+  BPolyn<BDat>& phi = Pol(Arg(2));
+  contens_.Alloc(acvf.Rows()-phi.Degree(),1);
+  BARIMA::ACVF_ARFilter(acvf.Data(), phi, contens_.GetData());
+}
+
+//--------------------------------------------------------------------
 DeclareContensClass(BSet, BSetTemporary, BSetARMAPreliminaryEstimation);
 DefExtOpr(1, BSetARMAPreliminaryEstimation, "ARMAPreliminaryEstimation", 
   4, 6, "Matrix Real Real Real Real Real",
-  "(Matrix ACVF, Real p, Real q, Real N [, Real maxIter=100*q, Real eps = 0.1/Sqrt(N))",
+  "(Matrix ACVF, Real p, Real q, Real N "
+  "[, Real maxIter=100*q, Real eps = 0.1/Sqrt(N)])",
   "Estimates a preliminary AR(p)MA(q) non seassonal model \n"
   "  \n"