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;
}