5-18 二分法求多项式单根 (用不少人竟然用 浮点数 0.0 == 0 的用法,哎)
来源:互联网 发布:万网域名交易 编辑:程序博客网 时间:2024/05/24 06:30
二分法求函数根的原理为:如果连续函数f(x)在区间[a,b]的两个端点取值异号,即f(a)f(b)<0,则它在这个区间内至少存在1个根r,即f(r)=0。
二分法的步骤为:
- 检查区间长度,如果小于给定阈值,则停止,输出区间中点(a+b)/2;否则
- 如果f(a)f(b)<0,则计算中点的值f((a+b)/2);
- 如果f((a+b)/2)正好为0,则(a+b)/2就是要求的根;否则
- 如果f((a+b)/2)与f(a)同号,则说明根在区间[(a+b)/2,b],令a=(a+b)/2,重复循环;
- 如果f((a+b)/2)与f(b)同号,则说明根在区间[a,(a+b)/2],令b=(a+b)/2,重复循环。
本题目要求编写程序,计算给定3阶多项式f(x)=a3x3+a2x2+a1x+a0在给定区间[a,b]内的根。
输入格式:
输入在第1行中顺序给出多项式的4个系数a3、a2、a1、a0,在第2行中顺序给出区间端点a和b。题目保证多项式在给定区间内存在唯一单根。
输出格式:
在一行中输出该多项式在该区间内的根,精确到小数点后2位。
#include <stdio.h>#include <math.h>double a3, a2, a1, a0, a, b, fLeft, fRight, mid, fMid, left, right;double f(double x){return a3*x*x*x + a2*x*x + a1*x + a0;}int main(){int isZero = 0;scanf("%lf %lf %lf %lf %lf %lf", &a3, &a2, &a1, &a0, &a, &b);double eps = 1e-7;while(b - a > eps){mid = (a + b) / 2;fMid = f(mid);if(fabs(fMid - 0) < eps) {printf("%.2f", mid); isZero = 1;break;}if(fMid * f(a) > 0) a = mid;else b = mid;}if(!isZero){printf("%.2f", mid);}return 0;}
0 0
- 5-18 二分法求多项式单根 (用不少人竟然用 浮点数 0.0 == 0 的用法,哎)
- PAT 5-18 二分法求多项式单根 浮点数二分法
- 5-18 二分法求多项式单根
- 5-18 二分法求多项式单根
- 7-18 二分法求多项式单根
- 5-18 二分法求多项式单根 C语言版
- 5-18 二分法求多项式单根 (20分)
- 5-18 二分法求多项式单根 (20分)
- 5-18 二分法求多项式单根 (20分)
- PTA 5-18 二分法求多项式单根
- 5-18 二分法求多项式单根 (20分)
- 5-18 二分法求多项式单根 (20分)
- 7-18 二分法求多项式单根(20 分)
- 二分法求多项式单根
- 二分法求多项式单根
- 5-21 二分法求多项式单根 (20分)
- 循环-08. 二分法求多项式单根
- PAT:二分法求多项式单根
- windows下nginx安装、配置与使用
- 剑指offer——二进制中1的个数
- URL基本介绍
- android:excludefromRecents属性控制在不在recent列表中显示
- oracle三层嵌套分页
- 5-18 二分法求多项式单根 (用不少人竟然用 浮点数 0.0 == 0 的用法,哎)
- 2016.06.15廖雪峰JS__学习笔记(浏览器对象)__P11
- 100多个纯CSS动画图标
- swiftUI基础 滑动手势UISwipeGestureRecognizer
- 数据恢复问题分析及注意事项
- android 6.0上getResources.getColor遭弃用
- ajaxFileUpload ajax文件上传
- OpenCV中的RANSAC详解
- Android 使用animator实现fragment的3D翻转效果