repeat_s.cpp

从出版社求得的经典数值算法

//////////////////////////////////////////////////////
//	程序3.2 — 复化Simpson公式

#include <stdlib.h>
#include <stdio.h>
#include <conio.h>
#include <string.h>
//#include <graphics.h>
#include <math.h>

#include "expressi.cpp"

float RepeatS(char FxString[],float a,float b,int n)
{
	int i;
	float h,Value1,Value2;

	if(n<=0) return 0;
	h=(b-a)/n;
	if(CreateFx(FxString)) return 0;

	Value1=0;
	Value2=0;
	for(i=0;i<n;++i)
	   Value1+=f(a+(i+0.5)*h);
	for(i=1;i<n;++i)
	   Value2+=f(a+i*h);
	return((f(a)+f(b)+4*Value1+2*Value2)*h/6);
}

void main()
{
	float i1,a,b;
	int n;
	char FxString[200];

	printf("\nInput function,a,b,n: ");
	scanf("%s %f %f %d",FxString,&a,&b,&n);

	i1=RepeatS(FxString,a,b,n);

	printf("\n%f",i1);

	getch();
}

/*
	运行实例：

Input function,a,b,n: 1/(1+x) 0 1 10

0.693147
*/