关于计算方法的误差分析

灰袍巫师

如何用Python实现一下C语言代码
#include "stdio.h"
#include "math.h"
float I0=log(6)-log(5);
float I1;
int n=1;
int main()
{
        printf("I0=%f\r\n",I0);
        while(n<=20)
        {
                I1=1.0/n-5*I0;
                printf("I%d=%f\r\n",n,I1);
                I0=I1;
                n++;
        }
        return 0;
}

jackz007

1. import math
   
2. 
   
3. I0 = math . log(6) - math . log(5)
   
4. print('I00 = %f' % I0)
   
5. n = 1
   
6. while n <= 20:
   
7.     I1 = 1 / n - 5 * I0
   
8.     print('I%02d = %f' % (n , I1))
   
9.     I0 = I1
   
10.     n += 1

复制代码

    C 代码的变量 I0、I1 不能是 float，应该是 double：

1. #include <stdio.h>
   
2. #include <math.h>
   
3. double I0 = log(6)-log(5) ;
   
4. double I1                 ;
   
5. int n = 1                 ;
   
6. int main()
   
7. {
   
8.         printf("I0=%lf\n" , I0)                ;
   
9.         while(n<=20)
   
10.         {
    
11.                 I1 = 1.0 / n - 5 * I0          ;
    
12.                 printf("I%d = %lf\n" , n , I1) ;
    
13.                 I0 = I1                        ;
    
14.                 n ++                           ;
    
15.         }
    
16.         return 0;
    
17. }

复制代码

        只有这样，C 和 Python 代码的输出才能较好地一致起来。