d8r5.cpp

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

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

void main()
{
    //program d8r5
    //driver for routine balanc
    int i,j,np = 5;
    double a[6][6],r[6],c[6];
    a[1][1]=1.0; a[1][2]=100.0; a[1][3]=1.0; a[1][4]=100.0; a[1][5]=1.0;
    a[2][1]=1.0; a[2][2]=1.0;   a[2][3]=1.0; a[2][4]=1.0;   a[2][5]=1.0;
    a[3][1]=1.0; a[3][2]=100.0; a[3][3]=1.0; a[3][4]=100.0; a[3][5]=1.0;
    a[4][1]=1.0; a[4][2]=1.0;   a[4][3]=1.0; a[4][4]=1.0;   a[4][5]=1.0;
    a[5][1]=1.0; a[5][2]=100.0; a[5][3]=1.0; a[5][4]=100.0; a[5][5]=1.0;
    //cout<< norms
    for (i = 1 ;i<=np; i++)
	{
        r[i] = 0.0;
        c[i] = 0.0;
        for (j = 1; j<=np; j++)
		{
            r[i] = r[i] + fabs(a[i][j]);
            c[i] = c[i] + fabs(a[j][i]);
        }
    }
    cout<<endl;
    cout<<setw(5)<<"Rows:"<<endl;
    cout<<endl;
    for (i = 1; i<=np; i++)
        cout<< setw(10)<<r[i];
    cout<<endl;
	cout<<endl;
	cout<<"Columns:"<<endl;
	cout<<endl;
    for (i = 1; i<=np; i++)
	{
        cout<<setw(10)<<c[i];
    }
    cout<<endl;
	cout<<endl;
    cout<< " ********** balancing matrix **********"<<endl;
    balanc(a, np);
    //cout<< norms
    for (i = 1; i<=np; i++)
	{
        r[i] = 0.0;
        c[i] = 0.0;
        for (j = 1; j<=np; j++)
		{
            r[i] = r[i] + fabs(a[i][j]);
            c[i] = c[i] + fabs(a[j][i]);
        }
    }
    cout<<endl;
    cout<<"Rows:"<<endl;
	cout<<endl;
    for (i = 1; i<=np; i++)
        cout<<setw(10)<<r[i];
    cout<<endl;
	cout<<endl;
    cout<<"Columns:"<<endl;
    cout<<endl;
    for (i = 1; i<=np; i++)
	{
        cout<< setw(10)<<c[i];
    }
    cout<< endl;
}