simpson.cpp

数值分析的简单算法VC++实现

// Simpson.cpp: implementation of the Simpson class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "Simpson.h"

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

void Simpson::cal(double DATA[100][2])
{int i,len;
 double h;
 for(len=0;DATA[len][0]<999999;len++)
 {}
 h=2*(DATA[len-1][0]-DATA[0][0])/double(len-1);
 result=h*(DATA[len-1][1]+DATA[0][1])/6;
 for(i=1;i<len-1;i=i+2)
 {result=result+h*(2*DATA[i][1]+DATA[i+1][1])/3;
 }
 result=result-DATA[len-1][1]*h/3;
}