用二分法求方程的根

总时间限制: 1000ms    内存限制: 65536kB

题目链接：http://ica.openjudge.cn/function1/4/

描述

用二分法求下面方程在(-10, 10)之间的一个根。 
2x³- 4x²+ 3x- 6 = 0

输入

一个小于1的非负实数e，它的值表示所能允许的误差

输出

一个实数，其值为求得的一个根，要求精确到小数点后8位。
若该区间上没有根，则输出“No Solution”

样例输入

0

样例输出

2.00000000

提示对于一个连续函数f(x)，若f(a)*f(b) <= 0，则f(x)在区间[a, b]内至少有一个根。
特别的，对于一个单调的连续函数，上述定理得逆定理也成立
若[a, b]上有根，则可进一步考察根是否在 [a, (a+b)/2]内，或者在[(a+b)/2, b]内。

若b-a <= e 则可终止迭代，并认为(a+b)/2是一个近似解，将它输出

请使用double类型！

#include<stdio.h>
#include<math.h>
double f(double x);
double fun(double a,double b,double e);
int main()
{
    double e,x1,x2,x;
    x1=-10;x2=10;
    scanf("%lf",&e);
    x=fun(x1,x2,e);
    printf("%.8lf\n",x);
    return 0;
}
double f(double x)
{
    double y;
    y=2*x*x*x-4*x*x+3*x-6;
    return y;
}
double fun(double a,double b,double e)
{
    double fa,fb,fc,c;
    fa=f(a);
    fb=f(b);
    c=(a+b)/2;
    while(fabs(b-a)>e)
    {
        fc=f(c);
        if(fc==0)
        {
            break;
        }
        else if(fc*fa<0)
        {
            b=c;
            fb=fc;
        }
        else
        {
            a=c;
            fa=fc;
        }
        c=(a+b)/2;
    }
    return c;
}

这个题目要用while语句实现才可以通过。下面的代码不能通过。（一直没懂什么原因……）下面这段是不行的。

#include<stdio.h>
#include<math.h>
double f(double x);
double fun(double a,double b,double e);
int main()
{
    double e,x1,x2,x;
    x1=-10;x2=10;
    scanf("%lf",&e);
    x=fun(x1,x2,e);
    printf("%.8lf\n",x);
    return 0;
}
double f(double x)
{
    double y;
    y=2*x*x*x-4*x*x+3*x-6;
    return y;
}
double fun(double a,double b,double e)
{
    double fa,fb,fc,c;
    fa=f(a);
    fb=f(b);
    c=(a+b)/2;
    do
    {
        c=(a+b)/2;
        fc=f(c);
        if(fc==0)
        {
            break;
        }
        else if(fc*fa<0)
        {
            b=c;
            fb=fc;
        }
        else
        {
            a=c;
            fa=fc;
        }
    }while(fabs(b-a)>e);
    return c;
}

 

下面的代码是通过了的。

#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
double f(double x);
double fun(double a,double b,double e);
int main()
{
    double e,x1,x2,x;
    x1=-10;x2=10;
    cin>>e;    
    x=fun(x1,x2,e);
    cout<<setprecision(8)<<setiosflags(ios::fixed)<<x<<endl; 
    return 0;
}
double f(double x) 
{
    double y;
    y=2*x*x*x-4*x*x+3*x-6;
    return y;
}
double fun(double a,double b,double e)
{
    double xm,y1,y2,ym;
    y1=f(a);
    y2=f(b);
    xm=(a+b)/2;
    while(fabs(a-b)>e)
    {
        ym=f(xm);
        if (ym==0) return xm;
        else  if(y1*ym<=0)
        {
               b=xm;
               y2=ym;
        }
        else
        {
               a=xm;
               y1=ym;
        }
        xm=(a+b)/2;                
    }
    return xm;
}