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;