用二分法求解x=cos(x)

[aozhi]

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

float countValue()
{
    float x=0.0;
    int n;
    for(n=2;;n=n<<1)
    {
        if(cos(x)-x>1e-6)
        {
            x=x+1.0/n;
        }
        else if(x-cos(x)>1e-6)
        {
            x=x-1.0/n;
        }
        else
            break;
    }
    return x;
}

int main(void)
{
    printf("%f\n",countValue());
    getch();
    return 0;
}

二分法：结果0.739085，循环次数20

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

float countValue() 
{
    float x0,x1=0.0;
    while(1)
    {
        x0=x1;
        x1=cos(x0);
        if(fabs(x0-x1)<1e-6)  
            break;
    }
    return x1;
}

int main(void)
{
    printf("%f\n",countValue());
    getch();
    return 0;
} 

普通迭代：结果0.739085，循环次数36
看来还是二分法好啊。

[aozhi]

[Quote=引用 4 楼 wangdeqie 的回复:]
楼主好强啊，学习～～～～
[/Quote]

大家都别听他的，他才是高手。

[zgjxwl]

=_=

[青葱岁月不可追]

aozhi

[aozhi]

干嘛?

[luxiaoxun]

up

[wangdeqie]

楼主好强啊，学习～～～～

[aozhi]

不要麻痹我.嘿嘿
有时间再搞个黄金分割版.呵呵

[青葱岁月不可追]

没干嘛，打个招呼！~

嘿嘿！~要往测试方向发展，当然得整几个测试方面的朋友啊！！~

嘿嘿！~

代码很强！~

[aozhi]

谢谢吉尔.晚点给分,多秀会儿.呵呵.

[青葱岁月不可追]

else if(x-cos(x)>1e-6)

是不是不用if了？？？？效率更高

[青葱岁月不可追]

没看清楚，请无视我8楼的

[aozhi]

厚厚.8楼的意思是?

[aozhi]

主要是算法效率提高,语句效率关系没怎么变。

[aozhi]

ls各位再见，上班去了。