计算物理作业

原创于 2019-10-24 16:39:33 发布 · 752 阅读

CC 4.0 BY-SA版权

本文提供了一个使用C语言实现的数值积分程序，通过复合梯形法则和辛普森法则计算从1到7区间内函数1/(2*sqrt(x))的积分值，并比较了不同分割数量下积分结果的准确性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#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;
}