Re: A problem with NIntegrate
- To: mathgroup at smc.vnet.net
- Subject: [mg119176] Re: A problem with NIntegrate
- From: Bob Hanlon <hanlonr at cox.net>
- Date: Mon, 23 May 2011 08:40:56 -0400 (EDT)
The magnitude of the imaginary part is less than 10^-25 on my machine, i.e., less than machine precision. These are arifacts of numerical noise. Look up Chop. Table[ NIntegrate[((1 - (u/40)^2.5)* ((1/((1 + 0.06*(u - t)))))/ Exp[0.01*u*(1.5 - 1)])^(1/1.5), {u, t, 40}], {t, 0, 40}] // Chop {19.7871, 19.2601, 18.7278, 18.1903, 17.6475, 17.0997, 16.5472, 15.99, \ 15.4286, 14.8633, 14.2945, 13.7225, 13.1479, 12.5712, 11.993, 11.4137, \ 10.8341, 10.2548, 9.6765, 9.1001, 8.52637, 7.95624, 7.39067, 6.83071, \ 6.27749, 5.73225, 5.19632, 4.67117, 4.15839, 3.65976, 3.17726, 2.71308, \ 2.26972, 1.85006, 1.45742, 1.09575, 0.769879, 0.485919, 0.252122, 0.0810213, 0} Or rationalize the powers Table[ NIntegrate[((1 - (u/40)^(5/2))* ((1/((1 + 0.06*(u - t)))))/ Exp[0.01*u*(1.5 - 1)])^(2/3), {u, t, 40}], {t, 0, 40}] {19.7871, 19.2601, 18.7278, 18.1903, 17.6475, 17.0997, 16.5472, 15.99, \ 15.4286, 14.8633, 14.2945, 13.7225, 13.1479, 12.5712, 11.993, 11.4137, \ 10.8341, 10.2548, 9.6765, 9.1001, 8.52637, 7.95624, 7.39067, 6.83071, \ 6.27749, 5.73225, 5.19632, 4.67117, 4.15839, 3.65976, 3.17726, 2.71308, \ 2.26972, 1.85006, 1.45742, 1.09575, 0.769879, 0.485919, 0.252122, 0.0810213, \ 0.} % == %% True Bob Hanlon ---- Anna Kaladze <anna.kaladze at gmail.com> wrote: ============= Hi all, I need to numerically integrate a function by utilizing =93Table=94 commandin the process (please see the code below), but for the parameter values I am interested in, I get complex results. I do not know why this is happening. Also, how do I get "normal" real values instead, and also ensuring that the answer is reasonably correct which one would expect from NIntegrate command. Thanks a lot. Table[NIntegrate[((1-(u/40)^2.5)*((1/((1+0.06*(u-t)))))/Exp[0.01*u*(1.5-1)])^(1/1.5),{u,t,40}],{t,0,40,1}] -- Bob Hanlon