高斯消去法求解方程组

最新推荐文章于 2024-12-20 07:00:00 发布

Matrix_Designer

最新推荐文章于 2024-12-20 07:00:00 发布

阅读量926

点赞数 1

分类专栏：数值分析文章标签： float ini include file url

本文链接：https://blog.youkuaiyun.com/Matrix_Designer/article/details/5923924

版权

数值分析专栏收录该内容

22 篇文章

订阅专栏

#include "iostream.h"
void main()
{  
int i,j,n;
    float **a;
void guass(float **a,int n);
cout<<"请输入方程次数:";
cin>>n;
a=new float*[n];
for(i=0;i<n;i++)
  a=new float[n+1];
cout<<"请输入增广矩阵:"<<endl;
for(i=0;i<n;i++)
  for(j=0;j<n+1;j++)
   cin>>a[j];
guass(a,n);
cout<<"方程的解为："<<endl;
for(i=0;i<n;i++)
  cout<<"X"<<i+1<<"="<<a[n]<<endl; 
}
void guass(float **a,int n)
{
int i,j,k;
for(i=0;i<n-1;i++) [url=file://选/]file://选[/url]主元素
{
  for(flag=i,j=i+1;j<n;j++)
   if(a[j]>a[flag])
    flag=j;
  if(flag!=i)
   for(j=0;j<n+1;j++)
   {
    t=a[flag][j];
    a[flag][j]=a[j];
    a[j]=t;
   }
}
for(i=0;i<n;i++)
{
  for(j=n;j>=i;j--)
   a[j]=a[j]/a;
  for(k=i+1;k<n;k++)
   for(j=n;j>=i;j--)
    a[k][j]-=a[j]*a[k];
}
for(i=n-2;i>=0;i--)
  for(j=i+1;j<n;j++)
   a[n]-=a[j][n]*a[j]; 
}