#define LL long long
#define SZ(x) ((int)(x).size())
#define ALL(x) (x).begin(),(x).end()
#define REP(i,n) for(int i=0;i<n;i++)
#define REV(i,n) for(int i=n-1;i>=0;i--)
#define FOR(i,a,b) for(int i=a;i<=b;i++)
#define pri(a) cout<<a<<endl
#define prii(a,b) cout<<a<<" "<<b<<endl
#define priii(a,b,c) cout<<a<<" "<<b<<" "<<c<<endl
#define hlw printf("Hello World\n");
using namespace std;
const int INF = 1<<29;
const int MX = 1e5+10;
double p,q,r,s,t,u;
double f(double x)
{
return p * exp(-x) + q*sin(x) + r*cos(x) + s*tan(x) + t * x * x + u;
}
double bisection()
{
double hi = 1, lo = 0;
while(lo+1e-7 < hi)
{
double mid = (lo + hi) / 2;
if(f(lo) * f(mid) <= 0)
{
hi = mid;
}
else
{
lo = mid;
}
}
return (lo + hi) / 2;
}
int main()
{
while(cin>>p>>q>>r>>s>>t>>u)
{
if(f(0) * f(1) > 0)
{
printf("No solution\n");
}
else
{
printf("%.4lf\n", bisection());
}
}
return 0;
}
No comments:
Post a Comment