[Kummi]

Hello Foamers,
I need to execute the following integrations in OpenFOAM as member functions.

H2Ox = @(EA) -H2O*(-8/(48000^8)).*(((EA)-E0).^7).*FE_H2O(EA).*A_PVM(EA);
H2Oz = integral(H2Ox, E0, Inf);

CO2x = @(EA) -CO2*(-4/(78000^4)).*(((EA)-E0).^3).*FE_CO2(EA).*A_PVM(EA);
CO2z = integral(CO2x, E0, Inf);

The above functions was executed in MATLAB. I don't know how to implement in OpenFOAM.
Kindly someone share your thoughts.
Thank you

[Kummi]

Hello all,
As per the following link,
can OpenFOAM calculate definite integral?
I tried this way using for-loop ~~

template<class thermo>inline Foam::scalar Foam::constSolidTransport<thermo>::C2H6(const scalar T) const
{
if (T>=623)
{
scalar dx = 1E-7, C2H6y =0, C2H6z = 0;

for (double x = 183000; x <= 1E100; x = x +dx)
{
if (*mC2H6_int < mC2H6_)
{
C2H6z = (SMALL + (-mC2H6_* ((x/x)*(-4/(VSMALL+SMALL+(Foam:ow(61000,4))))) * (Foam:ow((x-E0_),3)) * (Foam::exp(-(Foam:ow(((x-E0_)/61000),4)))) * ((1-(Foam::exp(-1.3E13*deltime_*Foam::exp(-x/(8.3144*T))))))));
*mC2H6_int += C2H6z * dx;
C2H6y = mC2H6_ - *mC2H6_int;
}
break;
}
std::cout<<"integral_out_C2H6 = " << C2H6y << std::endl;
return (C2H6y);
}
else { return mC2H6_;}
}

But that's not working properly..
Any thoughts Foamers ?
Thank you