数値計算と誤差 (1)

/* gcd.c
   3つの整数の最大公約数を求めるプログラム */

#include <stdio.h>

int gcd(int x, int y) {
  int r;
  r = x % y;
  while (r > 0) {
    x = y;
    y = r;
    r = x % y;
  }
  return y;
}

int main(void) {
  int s, t, u;
  scanf("%d", &s);
  scanf("%d", &t);
  scanf("%d", &u);
  printf("%d\n", gcd(gcd(s, t), u));
}

/* int_overflow.c
   整数計算での変な現象 */

#include <stdio.h>     // 入出力のライブラリを使う宣言

int main(void) {       // C 言語は関数 main から実行が始まる．
  int i;               // int で変数 i が整数型であることを示す．
  i = 1;               // 文の区切りを ; で表す．
  while (i > 0) {
    i = i + 1;
  }
  printf("%d \n", i);  /* 関数 printf の第1引数は，文字列で出力の形式を表す．
                            %d は10進数による数の表現．
                            \n は改行を表す．
			  第2引数以後は，%d などで表示される変数を現れる順番に記述する．*/
}

/* round-off_error.c
   丸め誤差の例 */

#include <stdio.h> 

int main(void) {
  double d1 = 1.0 - 0.9;
  double d2 = 0.1;
  if (d1 == d2) {
    printf("d1 is %g and d2 is %g, then d1 == d2 is true.\n", d1, d2);
  } else {
    printf("d1 is %g and d2 is %g, but d1 == d2 is false.\n", d1, d2);
    printf("d1 - d2 = %g\n", d1-d2);
  }
}

/* taylor.c
   情報落ちの例: log 1.5 の計算 */

#include <stdio.h>
#include <math.h>

int main(void) {
  int i, j;
  double s, x;
  s = 0;
  for (i = 1; i <= 100; i++) {
    x = i;
    for (j = 0; j < i; j++) {
      x = x * 2;
    }
    if (i % 2 == 0) {
      s = s - 1/x;
    } else {
      s = s + 1/x;
    }
  }
  printf("%.20f\n", s);
}