public void initTKF91(){
  int iNumNodes=iParent.length;
  double iBeta[]=new double[iNumNodes];
  iB=new double[iNumNodes];
  iE=new double[iNumNodes];
  iH=new double[iNumNodes];
  iN=new double[iNumNodes];
  iInitial=1.0;
  for (int i=0; i < iNumNodes; i++) {
    if (i == iNumNodes - 1) {
      iBeta[i]=1.0 / iMu;
      iH[i]=0.0;
    }
 else {
      iBeta[i]=Math.exp((iLambda - iMu) * iTau[i]);
      iBeta[i]=(1.0 - iBeta[i]) / (iMu - iLambda * iBeta[i]);
      iH[i]=Math.exp(-iMu * iTau[i]) * (1.0 - iLambda * iBeta[i]);
    }
    iB[i]=iLambda * iBeta[i];
    iE[i]=iMu * iBeta[i];
    iN[i]=(1.0 - iMu * iBeta[i]) * (1.0 - iB[i]) - iH[i];
    iInitial*=(1.0 - iB[i]);
  }
}

public void initTKF91 ( ) { int iNumNodes = iParent . length ; double iBeta [ ] = new double [ iNumNodes ] ; iB = new double [ iNumNodes ] ; iE = new double [ iNumNodes ] ; iH = new double [ iNumNodes ] ; iN = new double [ iNumNodes ] ; iInitial = 1.0 ; for ( int i = 0 ; i < iNumNodes ; i ++ ) { if ( i == iNumNodes - 1 ) { iBeta [ i ] = 1.0 / iMu ; iH [ i ] = 0.0 ; } else { iBeta [ i ] = Math . exp ( ( iLambda - iMu ) * iTau [ i ] ) ; iBeta [ i ] = ( 1.0 - iBeta [ i ] ) / ( iMu - iLambda * iBeta [ i ] ) ; iH [ i ] = Math . exp ( - iMu * iTau [ i ] ) * ( 1.0 - iLambda * iBeta [ i ] ) ; } iB [ i ] = iLambda * iBeta [ i ] ; iE [ i ] = iMu * iBeta [ i ] ; iN [ i ] = ( 1.0 - iMu * iBeta [ i ] ) * ( 1.0 - iB [ i ] ) - iH [ i ] ; iInitial *= ( 1.0 - iB [ i ] ) ; } }
populates the ib , ie , ih and in arrays , and iinitial value
