/* fib-double.c -- Fibonacci-Zahlen, iterativ (konstanter Platz), double */

#include <stdio.h>

double fibonacci (unsigned int n)
{
  if (n == 0) return 0;
  else if (n == 1) return 1;
  else {
    double a = 0;
    double b = 1;
    unsigned int i = 2;
    while (i <= n) {
      /* Invariante: a == F(i-2) und b == F(i-1) */
      double sum = a + b;
      a = b;
      b = sum;
      i = i + 1;
      /* Wieder korrekt */
    }
    /* Hier ist i == n + 1, also b == F(n) */
    return b;
  }
}

int main ()
{
  unsigned int n = 0;

  while (n <= 1500) {
    double fib_n = fibonacci(n);
    printf("F(%u) = %f (%e)\n", n, fib_n, fib_n);
    n = n + 1;
  }
  
  return 0;
}