数值方法源代码 C 插值积分微分迭代.docVIP

高斯算法和追赶法求解线性方程组 高斯： #includestdio.h #includestdlib.h #includemath.h void guass(double a[],int n,double b[],double x[]) { int i,j,m; double s; for(i=1,j=0;in;i++,j++) { for(m=i;mn;m++){ s=a[j*n+j]/a[m*n+j]; //printf(%f ,s); //printf(\n); for(int k=j;kn;k++) { a[m*n+k]=a[m*n+k]*s-a[(i-1)*n+k]; } b[m]=b[m]*s-b[i-1]; } } for(i=0;in;i++) for(int j=0;jn;j++){ printf(%f ,a[n*i+j]); if(j==n-1){ printf(%f ,b[i]); printf(\n);} } for(i=n-1;i-1;i--) { if(i==n-1) x[i]=b[i]/a[i*n+i]; else { s=b[i]; for(j=1;jn-i;j++) { s=s