前言
用下列方法求方程f(x)=e^x+10x-2的近似根,要求误差不超过0.5(e-3),并比较计算量:
在区间[0,1]上用二分法;
二分法:对于区间[a,b]上连续不断且f(a)·f(b)<0的函数y=f(x),通过不断地把函数f(x)的零点所在的区间一分为二,使区间的两个端点逐步逼近零点,进而得到零点近似值的方法(来源于百度百科)
代码如下
#include运行结果如下//输入与输出 #include #include ** **//数学函数调用 #include double f(double x)//定义数据类型 { return exp(x) + x * 10 - 2;//幂指数函数的调用 } int main() { double a = 0, b = 0, e = 0; printf("input a b e: "); scanf_s("%lf%lf%lf", &a, &b, &e); e = fabs(e); if (fabs(f(a)) <= e) { printf("solution: %lgn", a); } else if (fabs(f(b)) <= e) { printf("solution: %lgn", b); } else if (f(a) * f(b) > 0) { printf("f(%lg)*f(%lg) > 0 ! need <= 0 !n", a, b); } else { while (fabs(b - a) > e) { double c = (a + b) / 2.0;//取中间值 if (f(a) * f(c) < 0)//判断是否存在一个零点 b = c; else a = c; } printf("solution: %lgn", (a + b) / 2.0); } return 0; }
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)