[已解决]牛顿法求解

h-b-z-d-j-s-m

[code]#include<stdio.h>
#include<math.h>
#define PI 3.1415
double f(double x0)
{
        double x1;
        x1=(cos(x0)-x0)/(sin(x0)+1);
        x1=x1+x0;
        return x1;
}
int main()
{
        int t=0,k=100,n=0;
        double x0=PI/4,x1;
        while(n<k)
        {
                x1=f(x0);
                if(fabs(x0-x1)<0.001)
                {
                        t=1;
                        break;
                }
                else
                {
                        x0=x1;
                        n=n+1;
                }
                if(t==1)
                        printf("%10.5f",x1);
                else
                        printf("not found");
                return 0;
        }
}
这个答案是不是有问题啊，中间符号都不对

最佳答案

月排行榜 / 总排行榜

jackz007

2022-11-23 17:11:31

        验证计算结果是否正确的方法很简单，在输出 x1 的值后，多添加一句：

1.         printf("%.4f\n" , x1 - cos(x1)) ;  // 看看输出是不是 0 值就可以了。

复制代码

        这是我写的代码

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. #define PI 3.1415
   
5. 
   
6. double f(double x)
   
7. {
   
8.         return x + (cos(x) - x) / (sin(x) + 1) ;
   
9. }
   
10. 
    
11. int main(void)
    
12. {
    
13.         int t = 0 , k = 100 , n = 0 ;
    
14.         double x0 = PI / 4 , x1     ;
    
15.         while(n < k) {
    
16.                 x1 = f(x0)          ;
    
17.                 if(fabs(x0 - x1) < 1e-14) {
    
18.                         t = 1       ;
    
19.                         break       ;
    
20.                 } else {
    
21.                         x0 = x1     ;
    
22.                         n = n + 1   ;
    
23.                 }
    
24.         }
    
25.         if(t == 1) printf("\nFengcheng geng is %10.5f\n" , x1) ;
    
26.         else printf("\nSorry,not found!\n")                    ;
    
27. }

复制代码

       编译运行实况：

1. D:\[00.Exerciese.2022]\C>g++ -o x x.c
   
2. 
   
3. D:\[00.Exerciese.2022]\C>x
   
4. 
   
5. Fengcheng geng is    0.73909
   
6. 
   
7. D:\[00.Exerciese.2022]\C>

复制代码

跳转到最佳答案楼层

jackz007

        验证计算结果是否正确的方法很简单，在输出 x1 的值后，多添加一句：

1.         printf("%.4f\n" , x1 - cos(x1)) ;  // 看看输出是不是 0 值就可以了。

复制代码

        这是我写的代码

1. #include <stdio.h>
   
2. #include <math.h>
   
3. 
   
4. #define PI 3.1415
   
5. 
   
6. double f(double x)
   
7. {
   
8.         return x + (cos(x) - x) / (sin(x) + 1) ;
   
9. }
   
10. 
    
11. int main(void)
    
12. {
    
13.         int t = 0 , k = 100 , n = 0 ;
    
14.         double x0 = PI / 4 , x1     ;
    
15.         while(n < k) {
    
16.                 x1 = f(x0)          ;
    
17.                 if(fabs(x0 - x1) < 1e-14) {
    
18.                         t = 1       ;
    
19.                         break       ;
    
20.                 } else {
    
21.                         x0 = x1     ;
    
22.                         n = n + 1   ;
    
23.                 }
    
24.         }
    
25.         if(t == 1) printf("\nFengcheng geng is %10.5f\n" , x1) ;
    
26.         else printf("\nSorry,not found!\n")                    ;
    
27. }

复制代码

       编译运行实况：

1. D:\[00.Exerciese.2022]\C>g++ -o x x.c
   
2. 
   
3. D:\[00.Exerciese.2022]\C>x
   
4. 
   
5. Fengcheng geng is    0.73909
   
6. 
   
7. D:\[00.Exerciese.2022]\C>

复制代码