Re: integration problem???

*To*: mathgroup at smc.vnet.net*Subject*: [mg9306] Re: [mg9257] integration problem???*From*: David Withoff <withoff>*Date*: Mon, 27 Oct 1997 02:47:18 -0500*Sender*: owner-wri-mathgroup at wolfram.com

> I am trying to integrate the expression: ((2*r2-r1-r).^5).*(r.^(-8)) > with respect to r, between the limits r1=4.75e-4 to r2=2.25e-4. > > If I use Mathematica I get 1.9567e10 If I use quad8 in Matlab I get > 1.156e07 > > Does anyone know what is going wrong with this, please? > > Jason My guess is that there must be something fundamentally different between your inputs to Matlab and your inputs to Mathematica. The result from Mathematica In[8]:= r1=4.75*10^-4 Out[8]= 0.000475 In[9]:= r2 = 2.25*10^-4 Out[9]= 0.000225 In[10]:= NIntegrate[((2*r2-r1-r)^5)*(r^(-8)), {r, r1, r2}] 7 Out[10]= 1.15641 10 is correct. It would be surprising if Matlab would get the wrong answer for such an easy integral. Probably the Matlab result that you mentioned is somehow the result for a different integral. Dave Withoff Wolfram Research