d9r7.cpp

vc++常用数值算法集合

#include "iostream.h"
#include "math.h"
#include "stdlib.h"

void main()
{
    //program d9r7
    //driver for routine fleg
    int nval = 5;
    double x,dx = 0.2;
    int i,j,npoly = 5;
    double afunc[6];
    cout<<"Legendre polynomials"<<endl;
    cout<<"  n=1      n=2      n=3      n=4      n=5"<<endl;
	cout.setf(ios::fixed);
	cout.precision(4);
    for (i = 1; i<=nval; i++)
	{
        x = i * dx;
        fleg(x, afunc, npoly);
        cout<<" x = ";
		cout<<x<<endl;
        for (j = 1; j<=npoly; j++)
		{
			cout.width(11);
            cout<<afunc[j];
        }
        cout<<"   routine fleg"<<endl;
        for (j = 1; j<=npoly; j++)
		{
			cout.width(11);
            cout<<plgndr(j - 1, 0, x);
        }
        cout<<"   routine plgndr";
        cout<<endl;
    }
}