Chinaunix首页 | 论坛 | 博客
  • 博客访问: 1592337
  • 博文数量: 354
  • 博客积分: 8137
  • 博客等级: 中将
  • 技术积分: 5137
  • 用 户 组: 普通用户
  • 注册时间: 2009-11-26 15:40
文章分类

全部博文(354)

文章存档

2010年(300)

2009年(54)

分类: C/C++

2010-06-25 07:42:26

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

float f(float x)
{
    float y;
    y = ((x - 8.0)*x - 12.0)*x - 30.0;
    return y;
}

float xpoint(float x1, float x2)
{
    float y;
    y = (x1 * f(x2) - x2*f(x1)) / (f(x2) - f(x1));

    return y;
}
float root(float x1, float x2)
{
    float x, y, y1;
    y1 = f(x1);

    do{
        x = xpoint(x1, x2);
        y = f(x);
        if(y * y1 > 0)
        {
            y1 = y;
            x1 = x;
        }
        else
            x2 = x;
    }
    while(fabs(y) >= 0.0001);
    return x;
}

int main(int argc, char *argv[])
{
    float x1, x2, f1, f2,x;

    do{
        printf("Please input x1, x2:\n");
        scanf("%f, %f", &x1, &x2);
        f1 = f(x1);
        f2 = f(x2);
    }
    while(f1 * f2 > 0);
     x = root(x1, x2);

    printf("A root of equation is %9.6f\n", x);
    system("pause");
    return 0;
}


阅读(551) | 评论(0) | 转发(1) |
给主人留下些什么吧!~~