Re: Numerical integration and list of points
- To: mathgroup at smc.vnet.net
- Subject: [mg87654] Re: [mg87622] Numerical integration and list of points
- From: Curtis Osterhoudt <cfo at lanl.gov>
- Date: Tue, 15 Apr 2008 05:49:52 -0400 (EDT)
- Organization: LANL
- References: <200804141053.GAA13528@smc.vnet.net>
- Reply-to: cfo at lanl.gov
Hi, Gueromo,
It seems to work for me:
In[1]:= data = ({#1, Sin[13.4*#1]} & ) /@ Range[0, 4, 0.01];
f = Interpolation[data];
NIntegrate[f[x], {x, 0, 4}]
Integrate[f[x], {x, 0, 4}]
Integrate[Sin[a*x], x]
% /. a -> 13.4
(% /. {x -> 4}) - (% /. {x -> 0})
During evaluation of In[1]:= NIntegrate::ncvb: NIntegrate failed to \
converge to prescribed accuracy after 9 recursive bisections in x \
near {x} = {0.126211}. NIntegrate obtained 0.14786850709355207` and \
1.9590212682839143`*^-7 for the integral and error estimates. >>
Out[3]= 0.147869
Out[4]= 0.147869
Out[5]= -(Cos[a x]/a)
Out[6]= -0.0746269 Cos[13.4 x]
Out[7]= 0.147869
On Monday 14 April 2008 04:53:33 guerom00 wrote:
> Hello everyone,
>
> I have a function which I read as a list of points. I want then to
> estimate its integral. I do more or less this :
>
> data={{x1,y1},{x2,y2},...,{xN,yN}}
> f=Interpolation[data]
> NIntegrate[f[x],{x,x1,xN}]
>
> Is it the correct way ? Because Mathematica hangs without giving me an
> answer although it seems a pretty straightforward thing to do...
>
> Thanks for any suggestions.
--
==========================================================
Curtis Osterhoudt
cfo at remove_this.lanl.and_this.gov
PGP Key ID: 0x4DCA2A10
Please avoid sending me Word or PowerPoint attachments
See http://www.gnu.org/philosophy/no-word-attachments.html
==========================================================
- References:
- Numerical integration and list of points
- From: guerom00 <guerom00@gmail.com>
- Numerical integration and list of points