double InertiaCalculations::zMomentOfInertia(SQParameters& sq) {
double result = 0.0;
result = 0.5 * sq.a1 * sq.a2 * sq.a3 * sq.e1 * sq.e2 * ((sq.a1 * sq.a1) + (sq.a2 * sq.a2)) *
betaFunction(1.5 * sq.e2, 0.5 * sq.e2) * betaFunction(0.5 * sq.e1, (2 * sq.e1 + 1));
return result;
}
double
BetaRV::getPDFvalue(double rvValue)
{
double result;
if ( a <= rvValue && rvValue <= b ) {
double par1 = pow(rvValue-a,q-1.0);
double par2 = pow(b-rvValue,r-1.0);
double par3 = betaFunction(q,r);
double par4 = pow(b-a,q+r-1.0);
result = par1*par2/(par3*par4);
}
else {
result = 0.0;
}
return result;
}
double BetaDistribution::density(const double x) {
if(x<=0.0 || x>=1.0) return 0.0;
return pow(x,a_-1.0)*pow(1-x,b_-1.0) / betaFunction(a_,b_);
}
