泰勒级数展开求log（x)

#include "stdio.h" #include "math.h" double mathlog(double x); void main(void) { double M1=0,M2=0,sum=0,x; printf("Type x(0<x<=2^16)...\nx="); bool req=true,ack=true,dataready=true; while(ack&& dataready) //while循环实现握手协议功能 { dataready=false; scanf("%lf",&x); while(x>0) { if(req) { req=false; M1=mathlog(x); ack=true; } if(ack) { req=true; M2=mathlog(x+5); ack=false; } sum=M1+M2; printf("ln(%f)+ln(%f) = %f.\n",x,x+5,sum); dataready=true; ack=true; break; } } } double mathlog(double x) //mathlog函数功能是通过泰勒级数展开实现求ln(x)； { unsigned int i; double ln1_x,pox,dlt; if(x>-1.0 && x<=1.0) //实现x在(0 1)取值范围的ln(x)值 { for(ln1_x=0.0,dlt=pox=i=1;fabs(dlt)>1.0e-6;i++) { dlt=(pox*=x)/i; i&1 ? ln1_x+=dlt : ln1_x-=dlt; } printf("ln(%f) = %f.\n",x+1,ln1_x); } else if(x>1.0) // //实现x在(1 +00)取值范围的ln(x)值 { x=1/x; for(ln1_x=0.0,dlt=pox=i=1;fabs(dlt)>1.0e-8;i++) { dlt=(pox*=(x-1))/i; i&1 ? ln1_x+=dlt : ln1_x-=dlt; } ln1_x=-ln1_x; printf("ln(%f) = %f.\n",1/x,ln1_x); } else { printf("输入数据不在求值范围内,请重新输入"); } return ln1_x; }