C语言—定积分.docVIP

4.龙贝格求积公式，求解定积分 C/C++ code #includestdio.h #includemath.h #define f(x) (sin(x)/x) #define N 20 #define MAX 20 #define a 2 #define b 4 #define e 0.00001 float LBG(float p,float q,int n) { int i; float sum=0,h=(q-p)/n; for (i=1;in;i++) sum+=f(p+i*h); sum+=(f(p)+f(q))/2; return(h*sum); } void main() { int i; int n=N,m=0; float T[MAX+1][2]; T[0][1]=LBG(a,b,n); n*=2; for(m=1;mMaxims++) { for(i=0;im;i++) T[i][0]=T[i][1]; T[0][1]=LBG(a,b,n); n*=2; for(i=1;i=m;i++) T[i][1]=T[i-1][1]+(T[i-1][1]-T[i-1][0])/(pow(2,2*m)-1); if((T[m-1][1]T[m][1]+e)(T[m-1][1]T[m][1]-e)) { prints(Answer=%f\n,T[m][1]); getch(); return ; } } } 6. 牛顿-科特斯求积公式，求定积分 C/C++ code #includestdio.h #includemath.h int NC(a,h,n,r,f) float (*a)[]; float h; int n,f; float *r; { int nn,i; float ds; if(n1000||n2) { if (f) printf(\n Faild! Check if 1n1000!\n,n); return(-1); } if(n==2) { *r=0.5*((*a)[0]+(*a)[1])*(h); return(0); } if (n-4==0) { *r=0; *r=*r+0.375*(h)*((*a)[n-4]+3*(*a)[n-3]+3*(*a)[n-2]+(*a)[n-1]); return(0); } if(n/2-(n-1)/2=0) nn=n; else nn=n-3; ds=(*a)[0]-(*a)[nn-1]; for(i=2;i=nn;i=i+2) ds=ds+4*(*a)[i-1]+2*(*a)[i]; *r=ds*(h)/3; if(nnn) *r=*r+0.375*(h)*((*a)[n-4]+3*(*a)[n-3]+3*(*a)[n-2]+(*a)[n-1]); return(0); } main() { float h,r; int n,ntf,f; int i; float a[16]; printf(Input the x[i](16):\n); for(i=0;i=15;i++) scanf(%d,a[i]); h=0.2; f=0; ntf=NC(a,h,n,r,f); if(ntf==0) printf(\nR=%f\n,r); else printf(\n Wrong!Return code=%d\n,ntf); getch(); } C语言求函数定积分 2011-8-14 09:38 提问者： HYPERLINK /?businessaid=6un=%CC%EC%BF%D5%D6%AE%D6%F9 \l 2 \t _blank 天空之柱 | 浏览次数：185次 Auto.h中的内容： #includestdio.h #includemath.h #includestring.h #define g 9.81 #define e 2.71828182845904523536028747135266249775724709369995 #define F float #define SC scanf #define SCS scanf_s #define PR printf #define D %d\n #define PF %f\n #define D2 D D void swap( int * , int * ); void swap(int *ph1,int *ph2) { int temp; temp=*ph1; *ph1=*ph2; *ph2=tem