新机械原理第2版 教学课件 廖汉元 孔建益 GUASS.pptVIP

* 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： GUASS /**********/ #define pi 3.1415926 #include stdio.h #include math.h #include stdlib.h main() { int i,j,n=3,k,z,s; float x[10],b,c; static float a[10][10]= {{0,0,0,0,0},{0,3,2,1,14},{0,1,1,1,10},{0,2,3,-1,1}}; GUASS主元消去法程序 /**********/ #define pi 3.1415926 #include stdio.h #include math.h #include stdlib.h main() { int i,j,n=3,k,z,s; float x[10],b,c; static float a[10][10]= {{0,0,0,0,0},{0,3,2,1,14},{0,1,1,1,10},{0,2,3,-1,1}}; a[0][j]未用 GUASS主元消去法程序 /**********/ #define pi 3.1415926 #include stdio.h #include math.h #include stdlib.h main() { int i,j,n=3,k,z,s; float x[10],b,c; static float a[10][10]= {{0,0,0,0,0},{0,3,2,1,14},{0,1,1,1,10},{0,2,3,-1,1}}; GUASS主元消去法程序 = n=4 for(i=1;in+1;i++) for(j=1;jn+2;j++) printf(a[%d][%d]=%f\n,i,j,a[i][j]); for(j=1;jn+1;j++) { if(a[j][j]!=0) goto loop1; for(i=j;in+1;i++){if(a[i][j]!=0) goto loop2; } printf(no solution); goto loop4; loop2:for(k=j;kn+2;k++){b=a[j][k];a[j][k]=a[i][k];a[i][k]=b; } loop1:c=1/a[j][j]; for(z=1;zn+2;z++) a[j][z]=a[j][z]*c; for(i=1;in+1;i++) if(i!=j) {c=-a[i][j]; for(s=j;sn+2;s++) a[i][s]=a[i][s]+c*a[j][s]; } } loop3: printf(\n******solution******\n); for(i=1;in+1;i++) printf(x[%d][%d]=%f\n,i,n+1,a[i][n+1]); goto loop5; loop5:printf(finish\n); loop4:getch(); } X[1][4]?x1 X[2][4]?x2 X[3][4]?x3 /**********/ #define pi 3.1415926 #include stdio.h #include math.h #include stdlib.h main() { int i,j,n=4,k,z,s,numb=0; float x[10],b,c; static float a[10][10],f1,f2,f3,f4; float l1=20,l2=50,r3=50,l3=40,l4=30,l6=70.7,e=10, sda3=68*pi/180,fai1=30*pi/180,sda6=135*pi/180, fai2=100*pi/180,fai3=185*pi/180,fai4=45*pi/180, sy=55; 已知值 初始值 四元方程 迭代次数 例3-1 loop10: a[1][1]=l2*sin(fai2); a[1][2]=l3*sin(fai3);a[1][3]=l4*sin(fai4);a[1][4]=0;a[1][5]=-l1*cos(fai1)+l2*cos(fai2)+l3*cos(fai3)+l4*cos(fai4)-l6*cos(sda6); f1=a[1][5]; a[2][1]=-l2*