万年ダイエッターにささげる、“停滞期の正体”：世界を「数字」で回してみよう（22） ダイエット（10/10 ページ）

// gcc -g weight.c -o weight -std=c99
#include <stdio.h>
#include <stdlib.h>
int main ()
{
  // 体重変動シミュレーションプログラム
  //基礎代謝量と求め方	?ハリス・ベネディクト方程式(日本人版)　計算式  //http://www.kintore.info/kisotaisya_mass/index.cgi				
  //【計算式】
  // 男性 66.5+(体重kg×13.8)+(身長cm×5.0)-(年齢×6.8)
  // 女性 665.1+(体重kg×9.6)+(身長cm×1.9)-(年齢×7.0)							
  // 江端智一 2015日11月8日現在のデータを使って計算
#if 0   // 増量シミュレーション
  double weight = 66.5;  // 体重 66.5kg
  double lengthw = 172.0;  // 身長 172cm
  double age = 50.0 + 7.0/365;  // 年齢 50歳と7日目 (2015年11月8日現在)
#endif
  // 江端智一 2014日2月4日現在のデータを使って計算
#if 1   // 100日ダイエットシミュレーション
  double weight = 78.0;  // 開始時の体重 78.0kg
  double lengthw = 172.0;  // 身長 172cm
  double age = 49.0 + 94/365;  // 開始時の年齢 49歳と94日 (2014年2月4日現在)
#endif
  // 最初の基礎代謝カロリー 男性  66.5+(体重kg×13.8)+(身長cm×5.0)-(年齢×6.8)
  double consumption_calorie =  66.5 + weight * 13.8 + 172.0 * 5.0 - age * 6.8;
  // このカロリーに運動消費カロリーを加えて、トータルの一日の消費カロリーを算出する
  // (現在の体重 66.5kgにおいて)一日運動消費カロリー 708kcalと推定し、このカロリーは、体重に比例するものと仮定する。
  consumption_calorie += weight/66.5 * 708.0;
  printf("consumption_calorie = %f\n",  consumption_calorie);
#if 0 // 増量シミュレーション
  // 毎日、以下の量を食べると仮定する。
  // 私がマクドナルドで、ビッグマック、マックフライポテト(M)、チキンマックナゲット(5個)、コカ・コーラ(M)を頼むとすると、合計1,382kcal。ここに吉野家の牛丼(並盛)660Kcalを追加
  double source_calorie = 1382.0 * 2 + 660.0; 
#endif
#if 1 // 100日ダイエットシミュレーション
  // 江端のデータと一致する摂取カロリー
  double source_calorie = 1720;    // この値をなんどもいじって、実測データに近づける
#endif
  printf("source_calorie = %f\n",  source_calorie);
  for (int day=0; day < 1460; day++){  // 1460日は、ちょうど4年分の日数
	// 7kcal = 1g = 0.001kg となる
	double diff_weight = (source_calorie - consumption_calorie)/7.0 / 1000.0;
	weight += diff_weight;
	age += 1.0/365.0;
	consumption_calorie =  66.5 + weight * 13.8 + 172.0 * 5.0 - age * 6.8;
	consumption_calorie += weight/66.5 * 708.0;
	printf("day:%d\tweight = %f\tconsumption_calorie =%f\n", day,weight, consumption_calorie);
 }
}