void SystLinSurResolu::PushDifferentialEquation
     (
          Fonc_Num      aFonc,
          const PtsKD & aPts,
          REAL          aPds
     )
{
    INT anIndEq = mNbEqCur;
    SetNbEquation(mNbEqCur+1);

    mDataB[anIndEq] =  - aFonc.ValFonc(aPts);

    for  (INT iVar=0; iVar<mNbVarCur ; iVar++)
    {
          mDataA[anIndEq][iVar] = aFonc.ValDeriv(aPts,iVar); 
          mDataB[anIndEq] += mDataA[anIndEq][iVar] * aPts(iVar);
    }
    mDataPds[anIndEq] = aPds;
}