Question

求用c语言编写一个函数二分法求根的算法

或插值法也行

Answer 1

二分法计算函数f(x)=x*x*x*x+2*x*x*x-x-1;

本程序在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(h<10e-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;

}

}

}

Answer 2

二分法计算函数f(x)=x*x*x*x+2*x*x*x-x-1;
本程序在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(h<10e-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;
}
}
}

Answer 3

二分法计算函数f(x)=x*x*x*x+2*x*x*x-x-1;
本程序在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(h<10e-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;
}
}
}