解方程组的c++代码下载_Word模板_20

-.z.Cramer算法解方程组Gauss列主元解方程组Gauss全主元解方程组用Doolittle算法解方程组//解线性方程组*include*include*include//----------------------------------------------全局变量定义区constintNumber=15;//方程最大个数doublea[Number][Number],b[Number],copy_a[Number][Number],copy_b[Number];//系数行列式intA_y[Number];//a[][]中随着横坐标增加列坐标的排列顺序,如a[0][0],a[1][2],a[2][1]...则A_y[]={0,2,1...};intlenth,copy_lenth;//方程的个数doublea_sum;//计算行列式的值char**;//未知量a,b,c的载体//----------------------------------------------函数声明区voidinput();//输入方程组voidprint_menu();//打印主菜单intchoose();//输入选择voidcramer();//Cramer算法解方程组voidgauss_row();//Gauss列主元解方程组voidguass_all();//Gauss全主元解方程组voidDoolittle();//用Doolittle算法解方程组intDoolittle_check(doublea[][Number],doubleb[Number]);//判断是否行列式>0,假设是,调整为顺序主子式全>0void*iaoqu_u_l();//将行列式Doolittle分解voidcalculate_u_l();//计算Doolittle结果double&calculate_A(intn,intm);//计算行列式doublequanpailie_A();//根据列坐标的排列计算的值,如A_y[]={0,2,1},得sum=a[0][A_y[0]]*a[1][A_y[1]]*a[2][A_y[2]]=a[0][0]*a[1][2]*a[2][1];voide*change(intm,inti);//交换A_y[m],A_y[i]voide*change_lie(intj);//交换a[][j]与b[];voide*change_hang(intm,intn);//分别交换a[][]和b[]中的m与n两行voidgauss_row_*iaoqu();//Gauss列主元消去法voidgauss_all_*iaoqu();//Gauss全主元消去法voidgauss_calculate();//根据Gauss消去法结果计算未知量的值voide*change_a_lie(intm,intn);//交换a[][]中的m和n列voide*change_*(intm,intn);//交换*[]中的*[m]和*[n]voidrecovery();//恢复数据//主函数voidmain(){intflag=1;input();//输入方程while(flag){print_menu();//打印主菜单flag=choose();//选择解答方式}}//函数定义区voidprint_menu(){system("cls");cout<<"------------方程系数和常数矩阵表示如下:\n";for(intj=0;j>lenth;if(lenth>Number){cout<<"Itistoobig.\n";return;}*=newchar[lenth];for(i=0;i>a[i][j];cin>>b[i];}//备份数据for(i=0;i>choice;switch(choice){case1:cramer();break;case2:gauss_row();break;case3:guass_all();break;case4:Doolittle();break;case5:return0;default:cout<<"输入错误,请重新输入:";choose();break;}cout<<"\n是否换种方法求解(Y/N):";cin>>ch;if(ch=='n'||ch=='N')return0;recovery();cout<<"\n\n\n";return1;}//用克拉默法则求解方程.voidcramer(){inti,j;doublesum,sum_*;charch;//令第i行的列坐标为icout<<"用克拉默(Cramer)法则结果如下:\n";for(i=0;ii)l++;for(p=1,i=0;ia[ma*i][j])ma*i=i;if(ma*i!=k)e*change_hang(k,ma*i);//for(i=k+1;ia[ma*i][ma*j]){ma*i=i;ma*j=j;}}if(ma*i!=k)e*change_hang(k,ma*i);if(ma*j!=k){e*change_a_lie(ma*j,k);//交换两列e*change_*(ma*j,k);}for(i=k+1;i=0;i--){for(j=i+1,sum_a*=0;j=0;i--){for(j=i+1,sum_a*=0;jtemp_a[ma*i][j])ma*i=i;if(ma*i!=k){e*change_hang(k,ma*i);for(j=0;j

                    本文档为【解方程组的c++代码】，请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑，
                    图片更改请在作品中右键图片并更换，文字修改请直接点击文字进行修改，也可以新增和删除文档中的内容。 
 该文档来自用户分享，如有侵权行为请发邮件ishare@vip.sina.com联系网站客服，我们会及时删除。

                    [版权声明] 本站所有资料为用户分享产生，若发现您的权利被侵害，请联系客服邮件isharekefu@iask.cn，我们尽快处理。

                    本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权，请谨慎使用。

                    网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传，仅限个人学习分享使用，禁止用于任何广告和商用目的。
                

下载需要：免费已有0 人下载

立即下载

解方程组的c++代码

你可能还喜欢