uva 10341

uva 10341 Solve It

链接:http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=1282

题意:求所给的方程的根

思路:对给定的方程求导,可得其在区间[0, 1] 为单调函数,故在区间[0, 1] 最多只有一个解。当 f(0)*f(1) > 0 时, 无解。

  有解情况时,根据离散牛顿迭代法求方程的根。 离散牛顿法(割线法)公式如下:

xk+1 = xk - f(x) / (f(xk)-f(xk-1))  * (xk - xk-1)

代码:

 1 #include <iostream>
 2 #include <algorithm>
 3 #include <cstdio>
 4 #include <cmath>
 5 using namespace std;
 6 const double eps = 1e-9;
 7 inline int sig(double x){ return (x>eps) - (x<-eps); }
 8 
 9 int p, q, r, s, t, u;
10 double x;
11 inline double f(double x)
12 {
13     return p*exp(-x) + q*sin(x) + r*cos(x) + s*tan(x) + t*x*x + u;
14 }
15 
16 void solve()
17 {
18     double x1 = f(0.0), x2 = f(1.0);
19     if(x1 * x2 > 0)
20     {
21         puts("No solution");
22         return ;
23     }
24     x = 0.5, x1 = x2 = 0.0;
25     while(sig(x-x2))
26     {
27         x2 = x;
28         x = x - (f(x)*(x-x1)/(f(x)-f(x1)));
29         x1 = x2;
30     }
31     printf("%.4lf\n", x);
32     return ;
33 }
34 
35 int main()
36 {
37     //freopen("in.txt", "r", stdin);
38     while(scanf("%d %d %d %d %d %d", &p, &q, &r, &s, &t, &u) != EOF)
39         solve();
40     return 0;
41 }

转载于:https://www.cnblogs.com/Duahanlang/p/3382292.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值