c--解線性方程組的幾種方法-在線瀏覽
2024-09-05 10:39本頁面
【正文】 1. 克拉默(Cramer)法則。 cout\n 3. Gauss全主元消去法。 cout\n 5. 退出。}void input(){ int i,j。 cinlenth。 return。 for(i=0。i++) x[i]=39。+i。 cout請(qǐng)?jiān)诿總€(gè)方程里輸入lenth系數(shù)和一個(gè)常數(shù):\n。 for(i=1。i++) { cout+i+1x[i]。 cout應(yīng)輸入:。ilenth。 cout10\n。 //輸入每個(gè)方程 for(i=0。i++) { cout輸入方程i+1:。jlenth。 cinb[i]。ilenth。jlenth。 for(i=0。i++) copy_b[i]=b[i]。}//輸入選擇int choose(){ int choice。 cinchoice。break。break。break。break。 default:cout輸入錯(cuò)誤,請(qǐng)重新輸入:。 break。 cinch。n39。N39。 recovery()。 return 1。double sum,sum_x。 //令第i行的列坐標(biāo)為i cout用克拉默(Cramer)法則結(jié)果如下:\n。ilenth。 sum=calculate_A(lenth,0)。 for(i=0。i++) { ch=39。+i。 for(j=0。j++) A_y[j]=j。 sum_x=calculate_A(lenth,0)。 exchange_lie(i)。 } cout\n。 calculate_A(int n,int m) //計(jì)算行列式{ int i。 } else{for(i=0。i++) { exchange(m,m+i)。 exchange(m,m+i)。}double quanpailie_A() //計(jì)算行列式中一種全排列的值{ int i,j,l。 for(i=0,l=0。i++) for(j=0。amp。j++) if(A_y[j]i) l++。ilenth。 sum+=p*((l%2==0)?(1):(1))。}//高斯列主元排列求解方程void gauss_row(){ int i,j。 //用高斯列主元消區(qū)法將系數(shù)矩陣變成一個(gè)上三角矩陣 for(i=0。i++) { for(j=0。j++) coutsetw(10)setprecision(5)a[i][j]。 } if(a[lenth1][lenth1]!=0) { cout系數(shù)行列式不為零,方程有唯一的解：\n。 for(i=0。i++) //輸出結(jié)果 { coutx[i]=b[i]\n。}void gauss_row_xiaoqu() //高斯列主元消去法{ int i,j,k,maxi。 cout用Gauss列主元消去法結(jié)果如下:\n。klenth1。 for(maxi=i=k。i++) if(a[i][j]a[maxi][j]) maxi=i。// for(i=k+1。i++) { lik=a[i][k]/a[k][k]。jlenth。 b[i]=b[i]b[k]*lik。 gauss_all_xiaoqu()。ilenth。jlenth。 coutsetw(10)b[i]endl。 gauss_calculate()。ilenth。x[j]!=39。+iamp。jlenth。 coutx[j]=b[j]endl。}void gauss_all_xiaoqu() //Gauss全主元消去法{ int i,j,k,maxi,maxj。 cout用Gauss全主元消去法結(jié)果如下:\n。klenth1。i
黨政相關(guān)相關(guān)推薦
鄂ICP備17016276號(hào)-1