/*************************************************************************************************************************************/
/*                                                                                                                                   */
/*        A P P R O X I M A T I O N   R A T I O N N E L L E   D ' U N   R E E L                                                      */
/*        P A R   D E C O M P O S I T I O N   E N   F R A C T I O N S   C O N T I N U E S  :                                         */
/*                                                                                                                                   */
/*                                                                                                                                   */
/*        Author of '$xtc/ApproximationRationnelle.21$c' :                                                                           */
/*                                                                                                                                   */
/*                    Jean-Francois COLONNA (LACTAMME, 20180828175854).                                                              */
/*                                                                                                                                   */
/*************************************************************************************************************************************/

#include  <stdio.h>
extern    double floor();

#define   N         10

main()
          {
          int       n;

          double    x=3.141592653589;

          double    a0=1,a1,a;
          double    b0=0,b1=1,b;

          double    en,q;

          a1=floor(x);
          en=x-a1;

          for       (n=1 ; n <= N ; n++)
                    {
                    q=floor(1/en);
                    en=(1/en)-q;
                    a=(q*a1)+a0;
                    b=(q*b1)+b0;

                    printf("%.0f/%.0f\n",a,b);

                    a0=a1;
                    b0=b1;

                    a1=a;
                    b1=b;
                    }
          }



Copyright © Jean-François COLONNA, 2021-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2021-2024.