Zaninetti Lorenzo

20 Reputation

3 Badges

15 years, 219 days

MaplePrimes Activity


These are questions asked by

The following spreadsheet contains an integral to be done

in the framework of the Maxwell–Jüttner distribution and therefore date back to 1911

integrale_juttner_primes.mw
 

restart;

with(IntegrationTools):

the following integrand represents the scaling of the Maxwell–Jüttner distribution (1911) in 3D

where x is the Lorentz factor which is equal or greater than 1

integrandx:= x*sqrt(x^2 - 1)*exp(-x);

x*(x^2-1)^(1/2)*exp(-x)

(1)

C:=1/int(integrandx,x=1..infinity);

1/(BesselK(0, 1)+2*BesselK(1, 1))

(2)

pdf:=C*integrandx;

x*(x^2-1)^(1/2)*exp(-x)/(BesselK(0, 1)+2*BesselK(1, 1))

(3)

the above pdf represents the probability density function of the Maxwell–Jüttner distribution

 

DF:=int(C*integrandx,x=1..xx) assuming xx>1;

int(x*(x^2-1)^(1/2)*exp(-x)/(BesselK(0, 1)+2*BesselK(1, 1)), x = 1 .. xx)

(4)

the above DF represents the inert form of the distribution function but the integral at the moment

does not exists

 

 


 

Download integrale_juttner_primes.mw

 

The command eval allows to simplify a complicated expression

in a more compact form for a later output in LATEX

In the example which follows I was able to insert A and B but non C in the expression.

There is already a post on this kind of topic but I failed to understand the details.

Perhaps a Maple worksheet of answer on this topic  would be useful !

bye Lorenzo
 

restart;

expression:=exp(-b*x^c/2)*((x^(-(3*c)/2 + a/2 + 1/2)*(c + a + 1)*b^(-(3*c + a + 1)/(2*c)) + c*x^(a/2 + 1/2 - c/2)*b^(-(c + a + 1)/(2*c)))*c*WhittakerM((-c + a + 1)/(2*c), (2*c + a + 1)/(2*c), b*x^c) + b^(-(3*c + a + 1)/(2*c))*x^(-(3*c)/2 + a/2 + 1/2)*WhittakerM((c + a + 1)/(2*c), (2*c + a + 1)/(2*c), b*x^c)*(c + a + 1)^2)/((a + 1)*(c + a + 1)*(2*c + a + 1));

exp(-(1/2)*b*x^c)*((x^(-(3/2)*c+(1/2)*a+1/2)*(c+a+1)*b^(-(1/2)*(3*c+a+1)/c)+c*x^((1/2)*a+1/2-(1/2)*c)*b^(-(1/2)*(c+a+1)/c))*c*WhittakerM((1/2)*(-c+a+1)/c, (1/2)*(2*c+a+1)/c, b*x^c)+b^(-(1/2)*(3*c+a+1)/c)*x^(-(3/2)*c+(1/2)*a+1/2)*WhittakerM((1/2)*(c+a+1)/c, (1/2)*(2*c+a+1)/c, b*x^c)*(c+a+1)^2)/((a+1)*(c+a+1)*(2*c+a+1))

(1)

``

(2)

expression_ABC:=eval(expression,[x^(-(3*c)/2 + a/2 + 1/2)*(c + a + 1)*b^(-(3*c + a + 1)/(2*c))=A,c*x^(a/2 + 1/2 - c/2)*b^(-(c + a + 1)/(2*c))=B,((a + 1)*(c + a + 1)*(2*c + a + 1))=C]);

exp(-(1/2)*b*x^c)*((A+B)*c*WhittakerM((1/2)*(-c+a+1)/c, (1/2)*(2*c+a+1)/c, b*x^c)+b^(-(1/2)*(3*c+a+1)/c)*x^(-(3/2)*c+(1/2)*a+1/2)*WhittakerM((1/2)*(c+a+1)/c, (1/2)*(2*c+a+1)/c, b*x^c)*(c+a+1)^2)/((a+1)*(c+a+1)*(2*c+a+1))

(3)

 

 

 

 


 

Download maple_primes_eval.mw

 

In dark matter cosmology the following integral in x  does not yet have an analytical solution. I report the integral in

Maple 2018 notation

int(1/sqrt((1+x)^3*a+(1-a)*(1+x)^(3+3*b)), x)

If someone has some smart idea is welcome

Page 1 of 1