public double Mu0(double T)
{
if (T == 0) return 0;
// Initialization of the integral.
double integralMu = 0;
int I = DVPLDOM.AdaptiveTimeDiscretization.DiscreteTime(T, this.CacheDates);
double delta = this.CacheDates[1] - this.CacheDates[0];
double[] b = B(0, I, delta);
for (int i = 0; i < I; i++)
{
integralMu += this.dDeltaT[I - i] * b[i];
}
integralMu *= this.sigma1SquaredTemp;
return -integralMu * delta;
}