#include #include int main (void) { float g, L, angle, theta, x, u, dxdt, dudt, t; /* 変数を宣言する */ g = 9.8; /* 定数を代入する */ printf ("Length [m], angle [deg.]?\n"); scanf ("%f,%f", &L, &angle); /* 条件を入力する */ theta = 3.14159 / 180.0 * angle; x = L * theta; /* 初期値を計算する */ u = 0.0; printf ("%9.3f %9.3f\n", 0.0, angle); /* 結果を出力する */ dxdt = u; /* 時間微分を計算する */ theta = x / L; /* 角度を計算する */ dudt = - g * sin (theta); /* 時間微分を計算する */ x = x + 0.01 * dxdt; /* 次の時刻の値を計算する */ u = u + 0.01 * dudt; t = 0.01; /* 経過時間を計算する */ theta = x / L; /* 角度を計算する */ angle = 180.0 / 3.14159 * theta; printf ("%9.3f %9.3f\n", t, angle); /* 結果を出力する */ return 0; }