重庆分公司,新征程启航
为企业提供网站建设、域名注册、服务器等服务
二分法的基本思路是:任意两个点x1和x2,判断区间(x1,x2)内有无一个实根,如果f(x1)与f(x2)符号相反,则说明有一实根。接着取(x1,x2)的中点x,检查f(x)和f(x2)是否同号,如果不同号,说明实根在(x,x2)之间,如果同号,在比较(x1,x),这样就将范围缩小一半,然后按上述方法不断的递归调用,直到区间相当小(找出根为止)!
成都创新互联公司是一家专业提供德城企业网站建设,专注与网站设计、成都网站制作、HTML5建站、小程序制作等业务。10年已为德城众多企业、政府机构等服务。创新互联专业的建站公司优惠进行中。
代码如下(已调试):
#include "math.h"
main()
{
float x,x1,x2;
float F(float x,float x1,float x2);
printf("请输入区间[x1,x2]\n");
scanf("%f%f",x1,x2);
printf("x=%f\n",F(x,x1,x2));
}
float F(float x,float x1,float x2)
{
float f,f1,f2;
do
{
f1=pow(x1,3)-x1-1.0;
f2=pow(x2,3)-x2-1.0;
}while(f1*f20); //确保输入的x1,x2使得f1,f2符号相反
do
{
x=(x1+x2)/2; //求x1,x2的中点
f=pow(x,3)-x-1.0;
if(f1*f0) //当f与f1符号相同时
{x1=x;f1=f;}
else if(f2*f0) //当f与f2符号相同时
{x2=x;f2=f;}
}while(fabs(f)1e-6); //判断条件fabs(f)1e-6的意思是f的值非常0
return x;
}
输入:1 1.5
则输出:x=1.324718
重点改动:
scanf("%lf,%lf",x1,x2);
#include stdio.h
#include math.h
int main()
{ double calculate(double y1,double y2,double fy1,double fy2);
double x1,x2,fx1,fx2;
double a;
printf("enter x1 x2:");
scanf("%lf,%lf",x1,x2);
fx1=(((x1+3)*x1-8)*x1+12)*x1-10;
fx2=(((x2+3)*x2-8)*x2+12)*x2-10;
printf("%f %f %f %f\n",x1,x2,fx1,fx2);
a=calculate(x1,x2,fx1,fx2);
printf("x=%lf\n",a);
return 0;
}
double calculate(double y1,double y2,double fy1,double fy2)
{ double fy0,y0;
double a;
y0=(y1+y2)/2;
fy0=(((y0+3)*y0-8)*y0+12)*y0-10;
printf("%f ",y0);
if(fabs(fy0)=1e-10)
{ if ((fy0*fy1)0)
{ y2=y0;
fy2=fy0;
}
else
{ y1=y0;
fy1=fy0;
}
a=calculate(y1,y2,fy1,fy2);
}
else
{ a=y0;
printf("%6.2f",a);
}
return(a);
}
#include cmath double func(double x){
return x * (x - 3.0) - 5.0;
}
double binary_search(double x1, double x2, double (*f)(double)){
static const double epsilon = 1.0e-14;
double y1 = f(x1); double y2 = f(x2);
bool b1 = (y1 0.0 y2 0.0);
bool b2 = (y1 0.0 y2 0.0);
while (b1 || b2)
{
double x = (x1 + x2) / 2;
double y = f(x);
if (fabs(y) = epsilon) { return x;}
((b1 == y 0.0) ? x1 : x2) = x;
}
return 0.0;
}
int main(int, char*[]){
double f = binary_search(1.5, 10, func);
return 0;
}
按照你的要求只实现了这个方程的,不能泛用,不过很容易修改,代码如下:
#includestdio.h
#includemath.h
float function(float x)
{
float result = 0.0;
result = ((x*x*x) - (3*x) + 1);
return result;
}
int main()
{
float a = 0;
float b = 1;
float e = (float)0.001;
float m = (a+b)/2;
while(1)
{
if (fabs((function(m)) = e))
{
printf("result is %lf ", m);
return 1;
}
else if (function(m) 0)
{
a = m;
}
else if (function(m) 0)
{
b = m;
}
m = (a+b)/2;
}
return 1;
}