Question: Integral with Meijer Functions

Hi all,

 

I wonder if any of you guys can figure out why this integral is taking more than 2 hours to return a result??

I had similar problems in the past that were fixed creating a procedure, changing the solver and in case of trigonometric functions, changing the argument from float to rational number. 

Here it goes:Meijer_question.mw

restart

 

``

 

ms := (1.141448075+9.645873109*10^(-11)*I)*(MeijerG([[], []], [[1/4, 1/4], [1/2, 0]], .3956862293*(-.70*x+1)^4)+(2.148399968-2.148399963*I)*MeijerG([[], []], [[1/4], [1/2, 1/4, 0]], -.3956862293*(-.70*x+1)^4)+(-12.48809431-3.188863063*10^(-9)*I)*MeijerG([[], []], [[0], [1/2, 1/4, 1/4]], -.3956862293*(-.70*x+1)^4)+(4.061500400*10^(-10)-8.649913391*I)*MeijerG([[], []], [[1/2], [1/4, 1/4, 0]], -.3956862293*(-.70*x+1)^4));

b := .7;

``

NULL

"H3p:=proc(ms,b) local Hcub;  Hcub := evalf(Int(diff((diff(ms, x))*(int((-x*b+1)*(int((ms')^2, x = 0 .. x)), x = 1 .. x)), x)*ms, x=0..1,method = _d01ajc)): return Hcub; end proc:   "

``

 

st := time(); ``*H3p(ms, b); time()-st

``


Thanks!

Download Meijer_question.mw

Please Wait...