public double Evaluate(double a, double b)
{
var gamma = new GammaFunction();
var result = (gamma.EvaluateWithLanczos(a) * gamma.EvaluateWithLanczos(b)) / (gamma.EvaluateWithLanczos(a + b));
return(result);
}
public void When_Solve_Gamma_Function_With_Lanczos_Approximation()
{
var gammaFunction = new GammaFunction();
var firstResult = gammaFunction.EvaluateWithLanczos(2.5);
var secondResult = gammaFunction.EvaluateWithLanczos(-0.5);
Assert.Equal(1.3293403881791388, firstResult);
Assert.Equal(-3.5449077018110287, secondResult);
}
public IActionResult ComputeGammaFunction([FromBody] GetGammaRequest request)
{
var gamma = new GammaFunction();
var result = gamma.EvaluateWithLanczos(request.A);
return(new OkObjectResult(new FunctionResponse {
Value = result
}));
}
