jisuanqi.cpp

这个是数字积分器 用ＶＣ＋＋６.０环境开发的 大家有用的着的尽量下载 共同学习

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

#include "stdafx.h"
#include "数值积分计算器.h"
#include "jisuanqi.h"
#include "math.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif
#define N 100
//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

jisuanqi::jisuanqi()
{

}

jisuanqi::~jisuanqi()
{

}

double jisuanqi::TiXing(double a, double b)
{

    double x,h,T;
    int i;

    h=(b-a)/N;
    T=0.0;
    x=a;
    for(i=1;i<N;i++)
	{  
		x=x+h;
        T=T+exp(x*x);
	}
    T=h/2*(exp(a*a)+2*T+exp(b*b));
    return T;
}

double jisuanqi::XPS(double a, double b)
{
    double T;
	double s1,s2,s3;
	double x,h;
  

	h=(b-a)/(2*N);
	s1=exp(a*a)+exp(b*b);
	s2=0;s3=0;
	for(int i=1;i<=2*N-1;i++)
	{
		x=a+i*h;
		if(i%2==0)
			s3=s3+exp(x*x);
		else
			s2=s2+exp(x*x);
	}
	T=h*(s1+4*s2+2*s3)/3;
	return T;
}


double jisuanqi::Guass(double a, double b)
{
    double T;
	double s1;
	double x,h;
 
	h=(b-a)/N;
	s1=0;
	for(int i=0;i<=N-1;i++)
	{
		x=a+(i+0.5)*h;
		s1=s1+exp(x*x);
	}
	T=h*s1;
	return T;
}