d5r1.cpp

visual c++查过能用数值算法 附带光盘里面的源代码

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

void main()
{
    //program d5r1
    //driver for routine eulsum
    int nval = 40;
	double x,sum,xpower,term;
    double wksp[40];
    //evaluate ln(1+x)=x-x^2/2+x^3/3-x^4/4... for -1<x<1
    cout<<endl;
    cout<<"how many terms in polynomial?"<<endl;
    //input mval , between 1 and 40, mval=0 to end
    int mval = 20;
    cout<<mval<<endl;
    if ((mval <= 0) || (mval > nval))
	{
		return;
	}
    cout<<"         x           actual      polynomial"<<endl;
    for (int i = -8; i<=8; i++)
	{
        x = i / 10.0;
        sum = 0.0;
        xpower = -1.0;
		cout.setf(ios::fixed|ios::right);
		cout.precision(6);
        for (int j = 1; j<=mval; j++)
		{
            xpower = -x * xpower;
            term = xpower / j;
            eulsum(sum, term, j, wksp);
        }
		cout.width(14);
        cout<<x;
		cout.width(14);
        cout<<log(1.0 + x);
		cout.width(14);
        cout<<sum<<endl;
    }
}