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));
}
