public double D(double deltaT) { return (2 * this.gamma * Math.Exp(this.gamma * (deltaT))) / ((this.alpha1Temp + this.gamma) * Math.Exp(2 * this.gamma * (deltaT)) + (this.gamma - this.alpha1Temp)); }