计算物理作业

#include<stdio.h> 
#include<math.h>

int main(){
	int i,n;
	float sum,a=1,b=7,h;
	for(n=2;n<101;n+=5){
		sum=0;
		h=(b-a)/n;
		for(i=0;i<=n-1;i++){
			sum += 1/(2*sqrt(a+i*h));
		}
		sum = h*(1/(2*sqrt(a))+1/(2*sqrt(b))+2*sum)/2;
		printf("%d %lf\n",n,fabs(sum-sqrt(7)+1));
	}
	return 0;
}

#include<stdio.h> 
#include<math.h>

int main(){
	int i,n;
	double sum,a=1,b=7,h,s1,s2;
	for(n=2;n<101;n+=5){
		sum=0;
		h=(b-a)/n;
		s1=1/(2*sqrt(a+h/2));
		s2=0;
		s1 += 1/(2*sqrt(a+i*h+h/2));
		s2 += 1/(2*sqrt(a+i*h));
		for(i=0;i<=n-1;i++){
			sum += 1/(2*sqrt(a+i*h));
		}
		sum = h*(1/(2*sqrt(a)))+4*s1+2*s2+1/(2*sqrt(b))/6;
		printf("%d %lf\n",n,fabs(sum-sqrt(7)));
	}
	return 0;
}