Question: is it possible to calculate this

restart;with(Statistics):

MartinPoisson := (p/((p+2*q*T)/(1+p*T+q*T^2)-(p*T+q*T^2)*(p+2*q*T)/(1+p*T+q*T^2)^2))^r*r^2*T^n/factorial(n);

Dist := subs(T=t,MartinPoisson) assuming t > 0;

MartinDist := Distribution(CDF = unapply(piecewise(t>1,Dist,0),t));

MartinDist := Distribution(CDF = unapply(Dist, t)) assuming t > 0; #Change 2

X:=RandomVariable(MartinDist);

MartinDensity := PDF(X,t);

MartinDensity := subs(t=x, MartinDensity);

MartinDensity := subs(n=a, MartinDensity);

MartinDensity := subs(r=b, MartinDensity);

MartinDensity := subs(a=t/v, MartinDensity);

MartinDensity := subs(x=g, MartinDensity);

MartinDensity := subs(b=1/v, MartinDensity);

MartinDensity := simplify(MartinDensity);

tm := simplify(int(expand(1/(rho*sqrt(2*Pi*g))*exp(-((X-theta*g)^2)/(2*g*rho^2))*MartinDensity),g=0..infinity));

a := int(exp(i*X*u)*tm, X = -infinity .. infinity);Martin := simplify(a, power) assuming rho > 0, theta > 0, v > 0;

Martin2 := subs(i^2=-1,subs(u=-i,Martin));

Please Wait...