d5r3.cpp

这个是常用的数值算法用VC编写的。相信对大家有用哦。

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

void poldiv(double u[], int n, double v[], int nv, double q[], double r[])
{
    for (int j = 1; j<=n; j++)
	{
        r[j] = u[j];
        q[j] = 0.0;
    }
    for (int k = n - nv; k>=0; k--)
	{
        q[k + 1] = r[nv + k] / v[nv];
        for (int j = nv + k - 1; j>=k + 1; j--)
		{
            r[j] = r[j] - q[k + 1] * v[j - k];
        }
    }
    r[nv] = 0.0;
}

void main()
{
    //program d5r3
    //driver for routine poldiv
    int i,n = 6;
    int nv = 4;
    double u[7], v[5], q[7], r[7];
    u[1] = -1.0;
    u[2] = 5.0;
    u[3] = -10.0;
    u[4] = 10.0;
    u[5] = -5.0;
    u[6] = 1.0;
    v[1] = 1.0;
    v[2] = 3.0;
    v[3] = 3.0;
    v[4] = 1.0;
    poldiv(u, n, v, nv, q, r);
	cout.setf(ios::fixed|ios::right);
	cout.precision(2);
    cout<<endl;
    cout<<"      x^0       x^1       x^2       x^3       x^4       x^5"<<endl;
    cout<<"Quotient polynomial coefficients:"<<endl;
    for (i = 1; i<=n; i++)
	{
		cout.width(10);
        cout<<q[i];
    }
	cout<< endl;
    cout<<"Expected quotient coefficients:"<<endl;
    cout<<"     31.00     -8.00      1.00      0.00      0.00      0.00"<<endl;
    cout<<"Remainder polynomial coefficients:"<<endl;
    for (i = 1; i<=nv; i++)
	{
		cout.width(10);
        cout<<r[i];
    }
	cout<<endl;
    cout<<"Expected quotient coefficients:"<<endl;
    cout<<"    -32.00    -80.00    -80.00      0.00"<<endl;;
}