File: [local] / OpenXM / src / asir-contrib / packages / sample / taylor (download)
Revision 1.2, Thu Feb 17 10:47:04 2000 UTC (24 years, 7 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;