/*********已知f(x) = cos(x) - x。 x的初值为3.14159/4,用牛顿法求解方程f(x)=0的近似值**/
/*********要求精确到10E－6************/
#include"stdio.h"
#include"math.h"
double f1(double x);//要求解的函数
double f2(double x);//要求解的函数的一阶导数函数
double newton(double x0,double e);//通用牛顿迭代子程序
int main()
{
    double x0=3.14159/2;
    double e=10E-6;
    printf("x=%f\n",newton(x0,e));
    getchar();
    return 0;

}
double f2(double x)
{
    return -sin(x)-1;

}
double f1(double x)
{
    return cos(x)-x;
}
double newton(double x0,double e)   
{
    double x1;
    do
    {
        x1=x0;
        x0=x1-f1(x1)/f2(x1);
}
    while(fabs(x0-x1)>e);
    return x0;
}

大家帮帮忙啊！！谢谢

答案对的啊，你看看

#include"stdio.h"
#include"math.h"
double f1(double x);//要求解的函数
double f2(double x);//要求解的函数的一阶导数函数
double newton(double x0,double e);//通用牛顿迭代子程序
int main()
{
    double x0=3.14159/2;
    double x;
    double e=10E-6;
    printf("x=%f\n",x=newton(x0,e));
    printf("%f\n",cos(x)-x);
    printf("%f\n",cos(0.739085)-0.739085);
    return 0;

}
double f2(double x)
{
    return -sin(x)-1;

}
double f1(double x)
{
    return cos(x)-x;
}
double newton(double x0,double e)   
{
    double x1;
    do
    {
        x1=x0;
        x0=x1-f1(x1)/f2(x1);
    }
    while(fabs(x0-x1)>e);
    return x0;
}

我刚才又运行了一次，结果还是0.739085，cos0.739085接近1，怎么可能是0.739085？？谢谢了

cos0.739085
如果0.739085是用degree(度)表示就是接近1，用rad(弧度)表示就是0.739085啦
，我用计算器算的。

谢谢，你不说我还真没有想到呢！！看了几天程序都没发现错误，原来是对的啊，多谢各位了！