非线性方程求根

/*假设本金为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);
#endif
	double 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%