c语言编程- 微积分程序

原题在这个链接里 http://hi.baidu.com/170418878/blog/item/08691054bdc70652574e0002.html 现在要使用trapezoidal rule 求下列函数的积分值， trapezoidal rule 如下，
（不能贴图片，贴个链接）请用c语言编出，求下列积分值得程序，并且用 n=10 ， 100 ， 1000 时的值和 exact solution 做比较

举报该问题

推荐答案推荐于2016-08-12

辛普森法

#include <stdio.h>
#include <math.h>
#define F(X) (4.0/(1+X*X))
static float EPS = 1.0E-14;
static int COUNT=1;
static double a=0.0, b=1.0;
static double M[32],T[32],S[32],C[32],R[32],E[32];

double GETM(int K)
{ unsigned long j,n=1;
double x=0,y=0,step=0;
for(j=0;j<K;j++) n*=2;
step = (b-a)/n;
x = a+step/2;
for(j=0;j<n;j++){y+=F(x);x+=step;}
return (y*step);
}

Simpson(float EPS)
{ int k=0;
T[0]= (b-a)*(F(a)+F(b))/2.0;
S[0]=T[0];
E[0]=1.0E10;
for(k=0;k<20;k++)
{ M[k] = GETM(k);
T[k+1]=(M[k]+T[k])/2.0;
S[k+1]=(4.0*T[k+1]-T[k])/3.0;
COUNT++;
E[k+1]=fabs(S[k+1]-S[k]);
if(k<3)continue;
if(E[k+1]<EPS)break;
}
return;
}

ShawResult()
{ int k;
system("cls");
printf("\n K M[K] T[K] S[K]");
printf(" E[K]");
printf("\n-------------------------------------------");
printf("-----------------------------------");
for(k=0;k<COUNT;k++)
printf("\n%2d %20.15lf%20.15lf%20.15lf%12.4e",k,M[k],T[k],S[k],E[k]);
printf("\n--------------------------------------------");
printf("-----------------------------------");
getch();
system("cls");
return;
}

SaveResult()
{ int k;
FILE * fp;
fp=fopen("Simpson.htm","w");
if(!fp) return;
fprintf(fp,"<html><head></head>");
fprintf(fp,"<body bgcolor = #006699 text = #FFff00>");
fprintf(fp,"<pre><font size=\"6\">");
fprintf(fp,"\n K M[K] T[K] S[K]");
fprintf(fp," E[K]");
fprintf(fp,"\n----------------------------------------------");
fprintf(fp,"----------------------------------");
for(k=0;k<COUNT;k++)
fprintf(fp,"\n%2d %20.15f%20.15f%20.15f%12.4le",k,M[k],T[k],S[k],E[k]);
fprintf(fp,"\n----------------------------------------------");
fprintf(fp,"----------------------------------");
fprintf(fp,"</pre></font></body>");
fclose(fp);
return;
}

main()
{ Simpson(EPS);
ShawResult();
SaveResult();
}

温馨提示：答案为网友推荐，仅供参考

当前网址：http://99.wendadaohang.com/zd/XWBjzjtX.html

其他回答

第1个回答 2019-04-08

#include<stdio.h>
#include<math.h>
void
main()
{
int
i,n;
float
x[10001],y[10001];
float
x0,xn,h,JiFeng;
/*f=x^3*/
printf("input
x0,xn,n:");
scanf("%f",&x0);
scanf("%f",&xn);
scanf("%d",&n);
h=(xn-x0)/n;
x[0]=x0;
for(i=0;i<=n;i++)
{
x[i]=x[0]+i*h;
y[i]=pow(x[i],3);
}
JiFeng=0.0;
for(i=0;i<n;i++)
{
JiFeng=JiFeng+y[i]+y[i+1];
}
JiFeng=JiFeng*h/2.0;
printf("\n");
printf("n=%d
JiFeng
is
%f",n,JiFeng);
}

相似回答

C语言编写一个矩形求微积分的程序答：/*被积函数由宏定义给出*/ define f(x) ((x)*(x)+2*(x))main(){ int i,n;float a,b,t,m,s=0;printf("请输入积分上下限");scanf("%f%f",&a,&b);printf("请输入精确度m");scanf("%f",&m);for(n=1;;n++){ t=(b-a)/n;if(t<=m)break;} for(i=0;i<n;i++...

c语言编程求助- 微积分程序答：include<math.h> void main(){ int i,n;float x[10001],y[10001];float x0,xn,h,JiFeng;/*f=x^3*/ printf("input x0,xn,n:");scanf("%f",&x0);scanf("%f",&xn);scanf("%d",&n);h=(xn-x0)/n;x[0]=x0;for(i=0;i<=n;i++){ x[i]=x[0]+i*h;y[i]=pow...

请微积分和C语言高手帮忙编下程序答：include <stdio.h> include <math.h> int dhrt(double a,double b,double h,double eps,double x[],int m,double(*f)()){ int n,js;double z,y,z1,y1,z0,y0;n=0;z=a;y=(*f)(z);while((z<=b+h/2.0)&&(n!=m)){ if(fabs(y)<eps){ n=n+1;x[n-1]=z;z=z+h...

微积分用C语言怎么表达编程答：c不是专门用来运算的，所以说你也不用期待有几个符号加上函数就能求出微分和积分。如果实在要算微积分只能用微元法，把间隔取得足够小，用数值方法算出微分和积分

C语言编程求函数y=sinx在(0,1)区间的定积分?答：include<math.h> define N 10000000 void main(){ double sum=0;int i;for(i=0; i<N; i++){ sum += sin((double)(i)/N)/N;} printf("%lf\n%lf",sum,1-cos(1));} N后面的0有点多了，不过这个数刚好能精确到小数点后6位。当然，这个程序的算法实在是太垃圾了，效率非常的低...

怎样编写一个C语言编程求Y=sin(√ax)+ln(a+x)。要求a和x从键盘输入。当...答：include <math.h> { float a,x,y;printf("enter the number a,x:\n");scanf("%f,%f"＆a,＆x);y=sin(sqrt(a*x))+log(a+x);printf("y=%f\n",y);} 例如：include <stdio.h> include <math.h> int main(){ baidouble a, x, y;scanf("%lf%lf", &a, &x);y = pow(...

C语言如何实现微积分运算答：计算微积分有很多数值逼近的算法，任何可以计算微积分的语言都是用这种方法比如插值多项式，构造数值积分来计算的。只有近似值，没有准确值。你需要自己编程，编运算方法来计算。具体的lz去参看相关的书籍，比如计算方法

c++的编程语言中有微积分 导数距阵等高等数学方面的函数吗答：语言本身不含任何函数，其中包括著名的printf函数（这都不属于c语言范畴）.c语言的一些函数，是由库的形式提供，比如math.h提供的数学库,stdio.h提供的基本io库。目前倒是没有比较权威的高等数学函数。如果楼主有兴趣，可以自己补充这方面。或者找相对较完善的库。参考资料：百度 ...

求解C语言编程y=1*1/2*1/3*1/4*1/5*1/6*...前十项之积,用while语句编写...答：include<stdio.h> main(){ int n=10,i;float s=1;for(i=1;i<=10;i++)s=(1/i)*s;printf("%f",s);} 这个程序应该不是很难吧，楼主要多思考下哦

大家正在搜

用高级程序设计语言编写的程序 c语言程序编程用c语言编写的代码程序 c语言编译程序是什么软件 c语言编写的程序被称为 C语言程序设计教程汇编语言程序 c语言程序设计基础 c语言程序100例