390(1): 名前は開発中のものです。 [sage] 02/04/04 15:41 ID:??? AAS
こんなかんじ
#include <stdio.h>
#include <math.h>

#define L0 1.0 //振子長さ

int main ()
{
　　double x=L0, y=0.0, vx=0.0, vy=0.0;
　　double dt=0.01, endTime=2.0, t=0.0;

　　int c=0;
　　while ( t < endTime )
　　{
　　　　vy -= 9.8 * dt;

　　　　double l = sqrt ( x*x + y*y );
　　　　double dx = x / l * ( L0 - l );
　　　　double dy = y / l * ( L0 - l );

　　　　vx += dx / dt;
　　　　vy += dy / dt;

　　　　x += vx * dt;
　　　　y += vy * dt;

　　　　if ( c % 10 == 0 )
　　　　{
　　　　　　printf ( "%3.1f[sec], (%3.2f, %3.2f), l=%f\n", t, x, y, l );
　　　　}
　　　　c++;
　　　　t += dt;
　　}
　　return 0;
}

411(7): 名前は開発中のものです。 [sage] 02/04/06 08:11 ID:??? AAS
>>410
Rは水平軸との角度で，vは角速度ですかね？
�冲=1.0で省略てな感じかな？

まっとうなのは以下のようになると思います

振子の中心のモーメント M ＝ ｍ ｇ cosθ × L
振子の中心の慣性モーメント I ＝ ｍ L^2
よって角加速度 ωdot ＝ M / I ＝ ｇ cosθ / L
あとは積分

while ( )
{
　　ω += �冲 * G * cosθ / L;
　　θ += dt * omg;

　　x = L * cosθ;
　　y = L * sinθ;

　　t += �冲;
}
これだと逆三角関数いりませんぜ

＞無視してるようでゴメンなさい
オナニーして公開しただけですの気にしないで下さい

418(2): 411 [sage] 02/04/06 12:23 ID:??? AAS
でこんな感じ

while ( )
{
　　//モーメント
　　M = mass * G * cosθ * L;

　　//外力Fが重心において回転方向に加わる場合
　　M += F * L;

　　//（Lを含まないこともあろうとあえてLは約分せず）
　　ω += �冲 * M / （ L^2 * mass ）;
　　θ += dt * omg;

　　x = ax + L * cosθ;
　　y = ay + L * sinθ;
}

と書きながら>>410が今やっと解りました(ｱﾌｫ
でもやっぱ角速度と角度で管理したほうがよくないですか？
�冲＝1.0にしてθ＝π/2 - R とすればほぼ同じ形になります

で，一つ違う所は，
×　v += -(G/L)*sin(R);
○　v += -G*sin(R);
じゃないかな・・・
次元も合うし．�冲＝1.0で考えると前者のｖは角速度ですよ（L＝1.0なら別ですが）

423(1): 名前は開発中のものです。 [sage] 02/04/07 08:24 ID:??? AAS
>>422
なんも始まらんよ
playerが古いのかな．．．

426(1): 名前は開発中のものです。 [] 02/04/07 10:32 ID:QkW4z/Hg(1) AAS
ひもと直角な方向の速度を何回も足してるわけでしょ
円を描く様にはなりませんよ
移動するたびにらせん状にひもは長くなる