Index: tol/bmath/bstat/logit.cpp
===================================================================
--- tol/bmath/bstat/logit.cpp	(revision af055cbb7760b746fd64e547daa52a59857caca1)
+++ tol/bmath/bstat/logit.cpp	(revision 664d02140293e4f91edd0c3ac572455c26ed2de3)
@@ -74,9 +74,20 @@
   clock_t tm_1, tm_0 = BTimer::Clocks();
   BReal tm_dif;
+
   for(iter=0; iter<BDat::MaxIter(); iter++)
   {
     H_.SetAllValuesTo(0);
     G.SetAllValuesTo(0);
-    p = X*B;
+    if(iter==0)
+    {
+      p.SetAllValuesTo(0);
+    }
+    else
+    {
+    //Std(BText("\nTRACE Logit B=\n")+B.Name());
+    //Std(BText("\nTRACE Logit X=\n")+X.Name());
+      p = X*B;
+    //Std(BText("\nTRACE Logit p=\n")+p.Name());
+    }
     oldLlh = lnLikelyhood;
     lnLikelyhood=0;