1076-内部收益率

#include<stdio.h>
#include<math.h>
#include<string.h>
#define n 100
double f(double a[],int T,double IRR)
{
	double NPV=0;
	for(int i=0;i<=T;i++)
	{
		NPV=a[i]/pow(IRR+1,i)+NPV;
	}
	return NPV;
}
int main()
{
	int T;
	double a[n],b,IRR;
	double  x=-1+1e-5,y=1e6;
	scanf("%d",&T);
	while(1)
	{
		if(T==0)
		break; 
		for(int i=0;i<=T;i++)
		{
			scanf("%lf",&a[i]);
		}
		while(1)
		{
			IRR=(x+y)/2;
	        if(f(a,T,x)*f(a,T,IRR)<0)  //二分法//
	        {
	        	y=IRR;
			}else if(f(a,T,y)*f(a,T,IRR)<0)
			{
				x=IRR;
			}else
			{
				printf("No\n");
				break;
			}
			if(fabs(x-y)<=0.000001)
			{
				printf("%.2lf\n",IRR);
				break;
			}
	    }
	    x=-1+1e-5,y=1e6;   //重新定义两个边界//
		scanf("%d",&T);
		memset(a,0,sizeof(double));
	}
	return 0;
}