gauss_xuanzy.h

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

#include"fangcheng.h"
class Guass_XuanZY:public fangcheng
{
      public:
             void jie(void);
};
void Guass_XuanZY::jie(void)
{
     int i,m,k,max_h=1;
     float beishu=0,he;
     for(i=1;i<=(N-1);i++)
     {
           for(m=i+1;m<=N;m++)
           {
                if(a[m][i]>a[m-1][i])max_h=i;              
           }
           if(max_h>i)change(i,max_h);
           if(a[i][i]==0){cout<<"Wrong!";break;}
           for(m=i+1;m<=N;m++)
           {
                beishu=a[m][i]/a[i][i];
                b[m]=b[m]-beishu*b[i];
                for(k=i+1;k<=N;k++)
                {
                     a[m][k]=a[m][k]-beishu*a[m][i];              
                }                  
           }
     }

     x[3]=b[3]/a[3][3];
     for(i=N-1;i>=1;i--)
     { 
          for(k=i+1;k<=N;k++)he+=b[i]*a[i][k];
          x[i]=(b[i]-he)/a[i][i];
     }
     cout<<"*********************测试到了这里**********************"; 
}