d8r6.cpp

vc++常用数值算法集合

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

void main()
{
    //program d8r6
    //driver for routine elmhes
    int i,j,np = 5;
    double  a[6][6],r[6],c[6];
    a[1][1]=1.0; a[1][2]=2.0; a[1][3]=300.0; a[1][4]=4.0; a[1][5]=5.0;
    a[2][1]=2.0; a[2][2]=3.0; a[2][3]=400.0; a[2][4]=5.0; a[2][5]=6.0;
    a[3][1]=3.0; a[3][2]=4.0; a[3][3]=5.0;   a[3][4]=6.0; a[3][5]=7.0;
    a[4][1]=4.0; a[4][2]=5.0; a[4][3]=600.0; a[4][4]=7.0; a[4][5]=8.0;
    a[5][1]=5.0; a[5][2]=6.0; a[5][3]=700.0; a[5][4]=8.0; a[5][5]=9.0;
    cout<< endl;
    cout<<"***** Original matrix *****"<< endl;
    cout<< endl;
    for (i = 1; i<=np; i++)
	{
        for (j = 1; j<=np; j++)
			cout<< setw(9)<<a[i][j];
			cout<<endl;
	}
    cout<< endl;
    cout<<"***** balance matrix *****";
    cout<< endl;
    balanc(a,np);
    for (i = 1; i<=np; i++)
	{
        for (j = 1; j<=np; j++)
        cout<< setw(9)<<a[i][j];
	    cout<<endl;
	}
    cout<< endl;
    cout<< "*****reduce to Hessenherg form *****";
    cout<< endl;
    elmhes(a,np);
    for (j = 1; j<=np - 2; j++)
	{
        for (i = j + 2; i<= np; i++)
		{
			a[i][j] = 0.0;
		}
	}
	cout<<setprecision(4);
    for (i = 1; i<=np; i++)
	{
        for (j = 1; j<=np; j++)
		{
			cout<< setw(13)<<a[i][j];
		}
        cout<<endl;
	}
}