某IT外企面试题-二分法求方程！看看大家的C++功底

//f(x)=RMAX-R+(R-r)*double (cos(0.5*x))+double (sqrt(r*r-(R-r)*(R-r)* sin(0.500*x)* sin(0.500*x)))
请用二分法或其他方法求出解x.要求精度esp=1e-5
其中RMAX,R，r从键盘输入。


#include
#include
void get(double RMAX,double R,double r,double x,double y)
{

double m,n,c,p;
double esp=1e-5;

do
{
m=RMAX-R+(R-r)*double (cos(0.5*x))+double (sqrt(r*r-(R-r)*(R-r)* sin(0.500*x)* sin(0.500*x)));
n=RMAX-R+(R-r)*double (cos(0.500*y))+double (sqrt(r*r-(R-r)*(R-r)* sin(0.500*y)*sin(0.500*y)));
c=0.500*(x+y);
p=RMAX-R+(R-r)*double (cos(0.500*c))+double (sqrt(r*r-(R-r)*(R-r)* sin(0.500*c)*sin(0.500*c)));
if(fabs(p)
break;


else if(m*p<0)
{
x=x;
y=c;

}
else if (m*p>0)
{
y=y;
x=c;

}

}while(fabs(y-x)

cout<<\”The value of angle:tt\”<<(x+y)*0.5<
cout<<<\” \”<<<\” \”<<

}

void main()
{
double RMAX,R,r,x,y;
cout<<\”input RMAX,R,r,x,y:\”<

cin>>RMAX;
cin>>R;
cin>>r;
cin>>x;
cin>>y;

get(RMAX,R,r,x,y);
}