solve it--非线性方程求根

Solve It

Solve the equation:
        p*e^-x+ q*sin(x) + r*cos(x) + s*tan(x) + t*x² + u = 0
        where 0 <= x <= 1.

Input

Input consists of multiple test cases and terminated by an EOF. Each test case consists of 6 integers in a single line: p, q, r, s, t and u (where 0 <= p,r <= 20 and -20 <= q,s,t <= 0). There will be maximum 2100 lines in the input file.

Output

For each set of input, there should be a line containing the value of x, correct upto 4 decimal places, or the string "No solution", whichever is applicable.

Sample Input

0 0 0 0 -2 1

1 0 0 0 -1 2

1 -1 1 -1 -1 1

Sample Output

0.7071

No solution

0.7554

函数是单调递减函数,已知m在范围[0,1]内,运用二分得出一个值,代入方程,如果结果小于0,则m太大,如果结果大于0,则m太小如果等于0,则正好.

不知道为什么y-x>1e-5就错了....改成y-x>1e-10就对了

#include<stdio.h>
#include<algorithm>
#include<math.h>
using namespace std;
int main()
{
    double p,q,r,s,t,u,x,y,m,temp;
    while(scanf("%lf%lf%lf%lf%lf%lf",&p,&q,&r,&s,&t,&u)!=EOF)
    {
        if(p*exp(-0)+q*sin(0)+r*cos(0)+s*tan(0)+t*0*0+u<0||p*exp(-1)+q*sin(1)+r*cos(1)+s*tan(1)+t*1*1+u>0)
        {
            printf("No solution\n");
            continue;
        }
        x=0;
        y=1;
        while(y-x>=1e-10)
        {
            m=(x+y)/2.0;
            temp=p*exp(-m)+q*sin(m)+r*cos(m)+s*tan(m)+t*m*m+u;
            if(temp<0)
            {
                y=m;
            }
            else
                x=m;
        }
        printf("%.4lf\n",x);
    }
    return 0;
}