sharonyue 
April 7, 2014 23:17 
Quote:
Originally Posted by chrisb2244
(Post 484535)
In terms of a for loop, what I mean is something conceptually like
Code:
for (double x = 0.5; x <= 1; x = x+0.01)
{
sum += ( pow(x, 0.55) * exp(a * x));
}
where a is a volScalarField equal to the product of your

This method is feasible. But I think if OpenFOAM can deal with definite integral directly, that would be perfect. And the Gamma function you'd mentioned Im not familar with that, maybe it will take me a while to dig into it.
Actually Im tring to implant PBM into Foam, the RHS of the equation is luo breakage kernel~. Right now Im stuck here. I will try your loop method and keep this thread update. However, I dont know if it would be accurate, but its a good start I think. Thanks very much.
