算法小题二(求高次方程的一个实根2x^4-4x^3+6X^2-8x-8=0)

 求高次方程的一个实根2x^4-4x^3+6X^2-8x-8=0(精度要求10^-3) 

 算法分析如下:求高次方程的一个实根的算法有很多种,这里介绍一种二分法.

        设高次方程f(x)=0,则用二分法求其一个实根的算法如下

       1,输入初值x1和x2,计算它们对应的函数值y1 = f(x1), y2 = f(x2).若y1 和 y2同号,说明x1和x2之间没有实根转第一步重新输入;否则说明  x1,x2之间有实根,转为第二步.

       2,计算中间位置 x= (x1+x2)/2

       3,计算x1点x点函数值y1=f(x1),y=f(x)。若函数值y=0则算法结束，x就是这个高次方根的一个实根,否则转为4步

       4,若函数值y和y1同号,说明实根在x1和x2之间，则令x1=x,转为第2；否则，说明实根在x1到x之间,x2=x,转第2

#include<stdio.h>
#include <math.h>                       /*程序使用了求绝对值的函数fabs()*/
int main()
{
    float x1, x2, x , y1, y2, y;
    do
    {
        printf("进入do");
        scanf("%f,%f",&x1,&x2);          
        y1= 2*x1*x1*x1*x1-4*x1*x1*x1+6*x1*x1-8*x1-8;
        y2= 2*x2*x2*x2*x2-4*x2*x2*x2+6*x2*x2-8*x2-8;
        printf("y1*y2 = %d y1*y2>0为:%d",y1*y2,(y1*y2)>0);
    }while((y1*y2)>0);
    while(1)                           /*这是个条件恒真的死循环*/
    {
        x=(x1+x2)/2;
        y= 2*x*x*x*x-4*x*x*x+6*x*x-8*x-8;
        y1 = 2*x1*x1*x1*x1 -4 *x1*x1*x1+6*x1*x1-8*x1-8;
        if(fabs(y)<1e-5)break;        /*找到根后强迫退出*/
        else if(y*y>0)x1=x;
        else x2 =x;
    }
    printf("x=%f/n",x);               /*本方程在0-10间有个实根为2.0*/
}