public void RectangularInfititeTest() { var integral = new Integral(t => Pow(E, -t), 0, double.PositiveInfinity); var result = Integrals.RectangularInfinite(integral, 100000, 10); Assert.That(result, Is.EqualTo(1).Within(0.0001)); }
public virtual double HCoeffs(int s) { // TODO: use Gauss quadrature formula var coeffSym = mathlib.ScalarMul.LebesgueLaguerre(_right, Laguerre.Get(s)); return(Integrals.RectangularInfinite(coeffSym, 100000, 10)); }