xeulsum.c

来自「< C语言数值算法程序大全>>配套程序」· C语言 代码 · 共 39 行

/* Driver for routine eulsum */

#include <stdio.h>
#include <math.h>
#include "nr.h"
#include "nrutil.h"

#define NVAL 40

main()
{
	int i,j,mval;
	float sum,term,x,xpower,*wksp;

	wksp=vector(1,NVAL);
	/* evaluate ln(1+x)=x-x^2/2+x^3/3-x^4/4 ... for -1<x<1 */
	for (;;) {
		printf("\nHow many terms in polynomial?\n");
		printf("Enter n between 1 and %2d. (n=0 to end) ",NVAL);
		scanf("%d",&mval);
		printf("\n");
		if ((mval <= 0) || (mval > NVAL)) break;
		printf("%9s %14s %14s\n","x","actual","polynomial");
		for (i = -8;i<=8;i++) {
			x=i/10.0;
			sum=0.0;
			xpower = -1;
			for (j=1;j<=mval;j++) {
				xpower *= (-x);
				term=xpower/j;
				eulsum(&sum,term,j,wksp);
			}
			printf("%12.6f %12.6f %12.6f\n",x,log(1.0+x),sum);
		}
	}
	free_vector(wksp,1,NVAL);
	return 0;
}