Re: Area under Curve
- To: mathgroup at smc.vnet.net
- Subject: [mg7961] Re: [mg7943] Area under Curve
- From: Andrei Constantinescu <constant at athena.polytechnique.fr>
- Date: Wed, 30 Jul 1997 02:37:25 -0400
- Sender: owner-wri-mathgroup at wolfram.com
On Fri, 25 Jul 1997, Sankar Baruah wrote:
> I am new to Mathematica and would really appreciate some pointers from
> experienced users. I have a 2 variable data set of about 50
> observations. I want to plot each obseravtion and calculate the area
> under the curve. How do I do this in Mathematica ?
> Are there examples in the Mathematica Book and/or are there separate
> application notes etc. available ?
> Thanks.
> Sankar.
>
>
>
Hi Sankar,
there is no need to be experienced to solve this problem.
You have many solutions to do this (take also a look at the Mma books):
1. your data is really: {..., {xi, f[xi]}, ... }
- try ListIntegrate
- try an Interpolation + Integration
- write a simple Trapez system as the one below.
2. you data is {..., {xi,yi}, ... } - a closed curve.
- I do not know the function but, the following
simple function functioned pretty well:
Cint[ data_ ] := Module[{mylist,intg},
imax = Length[data];
mylist = Append[ data, data[[1]] ];
intg = 0;
Do[ intg += 0.25 (mylist[[i,2]] + mylist[[i+1,2]])
(mylist[[i+1,1]] - mylist[[i,1]]);
intg += - 0.25 (mylist[[i,1]] + mylist[[i+1,1]])
(mylist[[i+1,2]] - mylist[[i,2]])
, {i,imax} ];
Print[ intg / 4524., " [MJ /m^3]"]
]
a + andrei