Hi Sudarshan,
Plot your data in a grahp, and with the graph active, try the following code. You can copy-paste the code into script window, highlight-select all lines and then hit Enter to execute. You can then save it to OGS file etc - see help files for how to create script file.
Note that this just loops thru x values and computes area to chop up into one-thrids. This works best if there are enough points in the curve. If points are sparse, then area will not get divided into three equal halves etc.
Easwar
OriginLab
// Get size of active curve
limit %c;
nsize = limit.size;
// Find total area under curve
integ %c;
area=integ.area;
farea=area/3;
// Loop thru and find area for 1st one-third part
mks1=1;
for(i=mks1+1; i<=nsize; i++)
{
mks2 = i;
integ %c;
if( integ.area >= farea ) break;
}
xx2 = mks2;
area1=integ.area;
// Loop thru and find area for 2nd one-third part
mks1=xx2;
for(i=mks1+1; i<=nsize; i++)
{
mks2 = i;
integ %c;
if( integ.area >= farea ) break;
}
xx3 = mks2;
area2 = integ.area;
// Find area of 3rd part
mks1 = xx3;
mks2 = nsize;
integ %c;
area3 = integ.area;
// Get x dataset of active curve
%a=xof(%c);
// Report results
type "Total area \t\t\t=" $(area);
type "Area between x=$(%a[1]) to x=$(%a[xx2]) \t=" $(area1);
type "Area between x=$(%a[xx2]) to x=$(%a[xx3]) \t=" $(area2);
type "Area between x=$(%a[xx3]) to x=$(%a[nsize]) \t=" $(area3);
// Reset data markers and refresh page
mks1 = -1;
mks2 = -1;
doc -uw;
Divide Area under a curve into equal parts
Posted - 11/07/2005 : 10:27:53 AM
