fangcheng.h

高斯选主元分解法 用于解方程

#include<iostream>
using namespace std;
//所有均以3元3次方程为例 
const int N=3;
/****************以下为类定义***************/ 


class fangcheng
{
      public:
             float a[N+1][N+1];
             float b[N+1];
             float x[N+1];
             fangcheng(void)
             {     
                   a[1][1]=2;a[1][2]=-1;a[1][3]=3;
                   a[2][1]=4;a[2][2]=2;a[2][3]=5;
                   a[3][1]=0;a[3][2]=2;a[3][3]=2;
                   b[1]=1;b[2]=4;b[3]=6;
                   x[1]=1;x[2]=1;x[3]=1;
                   
             }                 
             void get_a_b(void);       //用于获取方程系数 
             inline void put_x(void);         //给出计算结果 
             void shiyan(void);        //用于检验结果 
             void change(int n,int m);
}; 



/*************以下为函数定义******************/ 



void fangcheng::get_a_b(void)
{
    int i,k;
    for(i=1;i<=N;i++)
    {
        cout<<"\n输入第"<<i<<"个方程的系数：";             
        k=1;
        for(;k<=N;k++)
        cin>>a[i][k];
     }
    cout<<"\n输入b的值：" ;
    for(i=1;i<=N;i++)
    {
        cin>>b[i];
     }              
}
void fangcheng::put_x(void)
{
     cout<<"\n计算结果是：x="<<x[1]<<";y="<<x[2]<<";z="<<x[3];
}
void fangcheng::shiyan(void)
{
     cout<<"\n这个方程是：\n";
     int i,k;
     for(i=1;i<=N;i++)
     {
         cout<<a[i][1]<<"x+"<<a[i][2]<<"y+"<<a[i][3]<<"z=" <<b[i]<<"\n";            
     }
}
void fangcheng::change(int n,int m)
{
     int i;
     float vol;
     vol=b[n];b[n]=b[m];b[m]=vol;
     for(i=1;i<=N;i++)
     {
          vol=a[n][i];
          a[n][i]=a[m][i];
          a[m][i]=vol;
     }
}