二分法计算函数f(x)=x*x*x*x+2*x*x*x-x-1;
创新互联建站长期为千余家客户提供的网站建设服务,团队从业经验10年,关注不同地域、不同群体,并针对不同对象提供差异化的产品和服务;打造开放共赢平台,与合作伙伴共同营造健康的互联网生态环境。为青海企业提供专业的成都网站设计、成都网站制作,青海网站改版等技术服务。拥有10多年丰富建站经验和众多成功案例,为您定制开发。
本程序在turbo c或c++下编译:
#include "stdio.h"
#include math.h
float f(float x)
{float y;
y=x*x*x*x+2*x*x*x-x-1;
return y;
}
void main()
{float a=0,b=0,h,y,x;
int k,n0;
printf("please input qujian a and b");
scanf("%f%f%d",a,b,n0); /*输入含根区间a,b,循环次数n0 */
for(k=0;k=n0;k++)
{ x=(a+b)/2;
h=(b-a)/2;
y=f(x);
if(h10e-6||fabs(y)10e-6)
{ printf("k=%d,x=%f,y=%f",k,x,y);
break; } /*输出分半次数k,函数的根x,及x对应的函数值.*/
else
{if(f(a)*f(x)0)
b=x;
else a=x;
}
}
}
没听说过二分法能够求根,二分法只是对于有序数组的查找方法,并不能去求解任意方程的根(只可以求解在某个区间单调递增或递减的函数曲线对应方程的根),较为通用的方程求根方法是牛顿法。
用二分法求方程在(-10,10)之间的根:2x^3-4x^2+3x-6=0.
解:x1=x0=(x1+x2)/2=x2
程序:
#includestdio.h
#includemath.h
int main()
{
float x0,x1,x2,fx0,fx1,fx2;
do
{
printf("输入x1,x2的值:");
scanf("%f,%f", x1, x2);
fx1 = 2*x1*x1*x1 - 4 * x1*x1 + 3 * x1 - 6;
fx2 = 2 *x2*x2*x2 - 4 *x2*x2 + 3 * x2 - 6;
} while (fx1*fx20);
do
{
x0 = (x1 + x2)/2;
fx0 = 2 * x0*x0*x0 - 4 * x0*x0 + 3 * x0 - 6;
if (fx0*fx1 0)
{
x2 = x0;
fx2 = fx0;
}
else
{
x1 = x0;
fx1 = fx0;
}
} while (fabs(fx0)= 1e-5);
printf("x=%5.2f\n",x0);
return 0;
}
if(c*d0)
{
a=m;
c=n;
}
else
{
b=m;
d=n;
}
改为
if(m*d0)
{
a=m;
c=m;
}
else
{
b=m;
d=m;
}