非线性方程求根

来源：互联网发布：百度云管家有mac版吗编辑：程序博客网时间：2024/04/30 05:21

/*假设本金为2000元，分2个月还，每个月还1100，月利率为x，则(2000(1+x)-c)(1+x)-c=0,令f(x)=(2000(1+x)-c)(1+x)-c，当f(x)>0时，说明月利率太高了；当f(x)=0时说明月利率正确；当f(x)<0时，月利率过低。结果要求输出三位小数。 */#include<stdio.h>#include<string.h>#include<stdlib.h>const int maxn=1000;int main(){#ifndef ONLINE_JUDGE  freopen("in.txt","r",stdin);#endifdouble a,c;int b;scanf("%lf%d%lf",&a,&b,&c);double x=0,y=100;//扩大一百倍处理只需算小数点后几位，充分利用整数部分，使运算精确，避免浮点误差 while(y-x>1e-5)//虽然题目说保留三位小数，但为结果不受浮点误差影响，算多两位小数 {double m=x+(y-x)/2;double f=a;for(int i=0;i<b;i++) f+= f*m/100-c;//计算当月利率为m时，函数值为多少if(f<0) x=m;//小了，变大些，在[m,y]找 else y=m;//所猜数m太大了，在[x,m]找 }printf("%.3lf%%\n",x);//以下用来探究浮点误差,mingw5编译器 printf("\n%lf\n",x);printf("%.3lf\n",(int)(0.79450*1000+0.5)/1000.0);//输出结果为0.795，为四舍五入保留三位小数的结果 printf("%.3lf\n",0.79450);//输出为0.794，产生浮点误差 printf("%.3lf\n",0.79451);//输出为0.795，也为四舍五入的结果 return 0;}

样例：

2000 4 510

输出：

0.797%