/*   PROJEKT  Kettenbruch                            Seite               */
/*   Speichermodell = small                                              */
/*   Kettenbruchzerlegung                                                */
/*-----------------------------------------------------------------------*/

#include <stdlib.h>
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <dos.h>

    main(void)
    {
      // double PI   = 3.141592653589793238462643;
	 double PI   = 3.14159265359;
      int    N     ;      /* Entwicklungstiefe        */
      double z     ;      /* INPUT                    */
      double target;      /* kopie dto                */
      double zf    ;      /* frac(INPUT)              */
      int    zi    ;      /* int (INPUT)              */
      double a[100];      /* Kettenbruchkoeffizienten */
      double z0,n0,z1,n1; /* Rekursive Brueche        */
      int i        ;

      clrscr();
      N=10;

    /*******************************************************/
    /* 1)  ANALYSE Umwandlung Zahl in Kettenbruch          */
    /*******************************************************/

    z=PI;                         /* zu analysierende Zahl */
    z=(sqrt(5)-1)/2.0;

    target=z;
    printf("Kettenbruchzerlegung\n");
    printf("---------------------\n");
    printf("target: %lf\n",z);

    for (i=0; i<N; i++)
    {
      zi=(int) z;   /* int       */
      zf=z-zi   ;   /* frac      */
      a[i]=zi   ;   /* koeff     */
      z=1.0/zf  ;   /* neue Zahl */
      printf("a%d = %d\n",i,zi);
    }
    printf("\n\n");

    /**************************************************/
    /*     SYNTHESE Umwandlung Kettenbruch in Bruch   */
    /**************************************************/
    /**************************************************/
    /* INIT Ermittle Zaehler/Nenner z1,n1 aus a[N-1]   */
    /***************************************************/
    z1=1.0; n1=a[N-1];
    /**************************************************/
    /* Rekursives Bestimmen Nenner/Zaehler            */
    /**************************************************/
    i=N-1;
    while (i--)
    {
    z0=n1;
    n0=a[i]*n1+z1;
    z1=z0;n1=n0;
    }
    /**************************************************/
    /*                  Ergebnis                      */
    /**************************************************/
    z0=n1;n0=z1;
    printf("\nERGEBNIS: %lf/%lf = %.11lf\n",z0,n0,z0/n0);
    printf("\nDIFFERENZ %.11lf",z0/n0-target);
    getch();
    return 0;

    }