>>14
> /* 理論上の微分方程式の一般解 */
> double differentiation_function(double t){return ( theta_i*cos(sqrt(g/r)*t*RADIAN) );}
DEG → RAD 化係数の場所が違うんでない？

return (theta_i *RADIAN) * cos(sqrt(g/r)*t);