@Klausklabauter 

It's not clear (at least for me) what is your final goal.
If you want the LDL decomposition then you can obtain it with 3 lines of code as shown. The complexity compared with a direct LDL algorithm is almost the same. You must have huge matrices to see a (marginal) difference.

 

@Annonymouse 

Basis returns the primitive parts of the polynomials, i.e. applies the primpart procedure (see ?primpart) to the "true" reduced basis. Note that this "basis" is also unique and has the advantage of not having non-integer coefficients.

It would be interesting to try this for t=1/10 instead of t=1.

@Mariusz Iwaniuk 

@weidade37211 

restart;
interface(version);
 Standard Worksheet Interface, Maple 2018.1, Windows 7, June 8  2018 Build ID 1321769
with(Statistics):
p1 := RandomVariable(BetaDistribution(1, 100)):
p2 := RandomVariable(BetaDistribution(1, 50)):
CodeTools:-Usage(
evalf(Int(z*PDF(0.3*p1+0.7*p2, z, inert), z = 0 .. 1))
);
memory used=0.56GiB, alloc change=142.00MiB, cpu time=6.33s, real time=6.18s, gc time=530.40ms
                         0.01669578722

@weidade37211 

As you see, this way it's computed fast too.

@aarjav 
It works for me.
BTW, note that Zeta is a built in function (Riemann Zeta). Probably you want zeta.

Here is a procedure for this.

Download IgnoreTerms.mw

@aarjav 

Your expressions are written incorrectly. You must use the argument t: so, theta(t) instead of theta etc.
Note that diff(theta, t) = 0, diff(beta, t, t) = 0 (without arguments).

[Actually it is possible to use alias(theta=theta(t)) but I'd recommend to avoid this at least for the moment].
Note also that zeta and Zeta are distinct objects.

It would be better to have either o formal description of the expression or a complete (not too long) example + the desired result.

@Rouben Rostamian  

Actually eq[2] is still undefined at (0,0,0). Maple does not see (without an extra simplify)  that the denominator is 0. 

@Carl Love 

The problem addressed by my answer was that strangely, fsolve is able to compute 4000 digits for z but not 5000.
Are you happy with this behavior?

It is strange that CubaCuhre works and other methods fail. After all, the function is well behaved so that any numerical integration should work. It seems that `evalf/MeijerG` is buggy.

@Carl Love 

Even when d=20 and we want 5000 digits for z (as above)?

@Mariusz Iwaniuk 

It seems that you have other integral in MMA.

@Markiyan Hirnyk 

Unfortunately you keep producing "anti-comments". It is a high school exercise to see that the equation 
arctan(x) = x/2 has 3 real roots and -tan(_Z) + 2 _Z = 0 has infinitely many real roots.

The important fact here is that RootOf uses a non-equivalent equation for the mentioned reason.

@Kitonum 

The name FrechetDistance for the procedure is misleading; this is actually a EuclideanDistance.
It would be interesting to compute the true FrechetDistance between the paths defined by P and Q.