int sum=0;
for(j=0;j<b;j++)
   sum+=eab[j]*p[j][i];
eab[i]=sum;

这怎么行呢?
在计算过程中eab[]已发生改变
增加一个数组吧
for ( int j = 0; j < b; j++ )
    sum += eab[j] * p[j][i];
ans_new[i] = sum;