[BACK]Return to taylor CVS log [TXT][DIR] Up to [local] / OpenXM / src / asir-contrib / packages / sample
File: [local] / OpenXM / src / asir-contrib / packages / sample / taylor (download)

Revision 1.2, Thu Feb 17 10:47:04 2000 UTC (24 years, 3 months ago) by tam
Branch: MAIN
CVS Tags: maekawa-ipv6, R_1_3_1-2, RELEASE_1_3_1_13b, RELEASE_1_2_3_12, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, RELEASE_1_1_3, RELEASE_1_1_2, KNOPPIX_2006, HEAD, DEB_REL_1_2_3-9
Changes since 1.1: +12 -5 lines 

I added test routine.

def taylor(F,X,P,N){
	R = 0;
	A = F;

	for(I=0;I<N;I++){
		R += subst(A,X,P) * x^I;
		A = diff(A,X)/(I+1);
	}
	return R;
}

def taylor_test(){
	print(taylor(exp(z),z,0,10));
	print(taylor(sin(z),z,0,10));
	print(taylor(cos(z),z,0,10));
	print(taylor(cos(z),z,1,10));
	print(taylor(1/(1-z),z,0,10));
	print(taylor((1-z)^(1/2),z,0,10));
	print(taylor((1-z)^(1/2),z,2,10));
}

/* taylor_test(); */
quit;