public double BSingle(int ti, int si, double delta) { // Initializes B integral. Vector b = new Vector(si); for (int i = ti; i < si; i++) b[i] = this.alphaT[i] / this.dDeltaT[si - i]; return 2 * this.dDeltaT[si - ti] * Fairmat.Math.Integrate.Trapezoid(b, delta); }