/**
 * 1ステップの時間発展を行います。
 */
void update(){

	// 質点間の距離
	float dx = LENGTH / N;

	// 質点1個あたりの質量
	float m = DENSITY * LENGTH / (N+1);

	// 端以外の頂点にかかる力の計算
	for(int i=1; i<N-1; i++){

		vertexFY[i] =
			(vertexY[i+1]+vertexY[i-1]-2.0*vertexY[i])
			* TENSION/dx - FRICTION*vertexVY[i];
	}

	// 両端にかかる力の計算
	if(FREE_END){ // 自由端

		// 隣接する1個の質点からのみ力を受ける
		vertexFY[0] = (vertexY[1]-vertexY[0])
			* TENSION/dx - FRICTION*vertexVY[0];

		vertexFY[N-1] = (vertexY[N-2]-vertexY[N-1])
			* TENSION/dx - FRICTION*vertexVY[N-1];

	}else{ //固定端

		// 壁から固定する抗力を受けるので、合力は0で動かない
		vertexFY[0] = 0.0;
		vertexFY[N-1] = 0.0;
	}

	// 位置と速度の計算
	for(int i=0; i<N; i++){

		// 頂点の速度変化の計算
		vertexVY[i] += vertexFY[i] / m * DT;

		// 頂点の座標変化の計算
		vertexY[i] += vertexVY[i] * DT;
	}
}