Disclaimer: This article is a blogger original article, follow the CC 4.0 BY-SA copyright agreement, reproduced, please attach the original source link and this statement.
#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;
}