一下午了,找不出错在哪里,有求高手了
实在是找不到错在哪里了,就是调试不过去。先谢了#include "math.h"
#include "stdio.h"
double r0=100.0,a=3.9083e-3,b=-5.775e-7,r=138.51;
double rcp(double x)
{
double r300,ra,k,ta;
r300=r0*(1+300*a+300*300*b);
k=(r300-r0)/300;
ra=x;
while(1)
{
ta=(ra-r0)/k;
ra=r0*(1+a*ta+b*ta*ta);
if((fabs(x-ra))<=0.0001)break;
}
return ta;
}