数值计算编程作业.docVIP

数值计算作业 \ 一. 为求方程X3-X2-1=0在X0=1.5附近的一个根，设将方程写 成X=1+1/X2，迭代格式为Xi+1=1+1/Xi2。 解：#includestdio.h main() {long float x,x0; X0=1.5; While(fab(x-x0)1e-3) {x=1+1/x0*x0; x0=x;} printf(“x=%lf\n”,x);} 运行结果：x=1.465 二.用牛顿法求f（x）=x3-3x-1=0在x0=2附近的根，要求有 四位有效数字。 解：#includestdio.h #includemath.h main() {long float x,x0,f,f1; X=2.0; do{x0=x; f=x0*x0*x0-3*x0-1; f1=3*x0*x0-3; x=x0-f/f1;} while(fabs(x-x0)=1e-4); printf(“The root of the EQU is%6.3lf\n”,x);} 运行结果：The root of the EQU is1.879 三. 用消去法解方程组 2x1+2x2+3x3=3 ① 4x1+7x2+7x3=1 ② -2x1+4x2+5x3=-7 ③ 解：#includeiostream #includecmath #define N 20 using namespace std; void load(); float a[N][N]; int m; m=3; int main() {int i,j,m; int c,k,n,p,r; float x[N],l[N][N],s,d; coutendl; cout请按顺序输入增广矩阵a:endl; load(); for(i=0;im;i++) {for(j=i;jm;j++) c=(fabs(a[j][i])fabs(a[i][i]))?j:i; for(n=0;nm+1;n++) {s=a[i][n]; a[i][n]=a[c][n]; a[c][n]=s;} for(p=0;pm+1;p++) couta[i][p]\t; coutendl; for(k=i+1;km;k++) {l[k][i]=a[k][i]/a[i][i]; for(r=i;rm+1;r++) a[k][r]=a[k][r]-l[k][i]*a[i][r];} } x[m-1]=a[m-1][m]/a[m-1][m-1]; for(i=m-2;i=0;i--) {d=0; for(j=i+1;jm;j++) d=d+a[i][j]*x[j]; x[i]=(a[i][m]-d)/a[i][i];} cout该方程组的解为:endl; for(i=0;im;i++) coutx[i]=x[i]\t; //system(pause); return 0;} void load() {int i,j; for(i=0;im;i++) for(j=0;jm+1;j++) cina[i][j];} 请按顺序输入增广矩阵a: 2 2 3 3 4 7 7 1 -2 4 5 -7 4 7 7 1 0 -1.5 -0.5 3 0 7.5 8.5 -7 运行结果：该方程组的解为：x[0]=2 x[1]=-2 x[2]=1 四.用梯形公式求∫19x1/2dx的积分。 解：#includestdio.h #includemath.h doublef(double x) {return (sqrt x);} double T(double a,double b) {double m,n; h=(b-a)/2,n=f(a)+f(b); Return(m*n);} main() {double a,b; a=1，b=9； double T(double a,double b); printf(“%lf”,T(a,b));} 运行结果：16 五.取步长h=0.2用龙格-库塔格式解初值问题。 y’=x+y 0≤x≤0.6 {y (0)=1 解：#includestdio.h #includemath.h #define f(x,y),x+y; main() {int n,k=1; long float a,b,h,x,y,y0,t[4]; scanf(“%lf,%lf,%lf,%lf”,a,b,y0,h); n=(b-a)/h; x=a; y=y0; printf(“%l