As all of you pointed that the result get converging with upper limit taken to be 100, And I am looking for 4-5 decimal places. I have a doubt what if in this the same case if it doesn't converge with these many points and we are asked to go for infinity or atleast 500 points. With 200 points it shows maple was unable allocate enough memory. Please see below code where a strange results are obtaining with Digits:=50
. Since I am a beginner with maple I am unaware of running the job in parallel. Also how is it possible to implement same in mathematica (where my approach takes 3-4 hours to complete the job) since I am much more comfortable with using mathematica.