d8r7.cpp

Visual C++ 常用数值算法集 源代码

#include <math.h>
#include <iomanip.h>
#include <iostream.h>

void main()
{
    //program d8r7
    //driver for routine hqr
    int i,j,np = 5;
    double  a[6][6],wr[6],wi[6];
    a[1][1]=1.0;  a[1][2]=2.0; a[1][3]=0.0;   a[1][4]=0.0; a[1][5]=0.0;
    a[2][1]=-2.0; a[2][2]=3.0; a[2][3]=0.0;   a[2][4]=0.0; a[2][5]=0.0;
    a[3][1]=3.0;  a[3][2]=4.0; a[3][3]=50.0;  a[3][4]=0.0; a[3][5]=0.0;
    a[4][1]=-4.0; a[4][2]=5.0; a[4][3]=-60.0; a[4][4]=7.0; a[4][5]=0.0;
    a[5][1]=-5.0; a[5][2]=6.0; a[5][3]=-70.0; a[5][4]=8.0; a[5][5]=-9.0;
    cout<< "Matrix:"<<endl;
    cout<<endl;
    for (i = 1; i<=np; i++)
	{
		cout<<setprecision(3)<<setiosflags(ios::fixed);
        for (j = 1; j<= np; j++)
            cout<<setw(12)<<a[i][j];
         cout<<endl;
    }
    balanc(a,np);
    elmhes(a,np);
    hqr(a,np,wr,wi);
    cout<<endl;
    cout<<"Eigenvalues:"<<endl;
    cout<<endl;
    cout<<"        Real          Imag."<<endl;
	cout<<setprecision(6)<<setiosflags(ios::fixed);
    for (i = 1; i<=np; i++)
	{
        cout<<setw(14)<<wr[i];
        cout<<setw(14)<<wi[i];
    	cout<<endl;
	}
}