步长.txt

利用变步长来计算积分

#include "math.h"
#include "stdio.h"

void main()
{
	double e=2.71828,z,f1,f2,r;
	float x,h;
	int i;
	printf("输入要计算的x的值:");
	scanf("%f",&x);
	z=pow(e,x);
	printf("\n*****@_@悄悄地告诉你它的标准值是f'(%lf)=%lf*****\n",x,z);
    printf("\n请您输入初始步长h的值:");
	scanf("%f",&h);
	for(i=0;;i++)
	{
		
        f1=(pow(e,x+h)-pow(e,x-h))/(2*h);
        f2=(pow(e,x+h/2)-pow(e,x-h/2))/h;
       
		r=f1-f2;
		  if(r<0) {r=-r;}
		
		printf("\nh=%f\t",h);
		printf("f'(%5.2f)=%lf\t\t",x,f1);
		printf("误差r=%f\t",r);
		h=h/2;
		
		if(r<0.00005)
		{
			printf("\n\n最佳步长是h=%f\t\t",2*h);
			printf("f'(%lf)=%lf\n",x,f1);
			break;
		}
			
	}
	printf("\n怎么样，跟标准值差不多吧，呵呵  @_@  \n");
	printf("\n");

}