page2
次に、さらに微小時間Δt秒後の速度と高さを計算する。
(8)
(9)
上記の計算を繰り返せば、任意の時刻（但し、離散値）での物体の位置を計算することができる。
時間の刻み幅を微小値Δtとし、速度vjと位置yjを漸化式で表す。また、t=jΔtとする。
(10)
(11)
j=1から順に(10)式、(11)式を繰り返すことにより、物体の高さを経時的に計算することができる。
これらの式を、Calcシートに記述し、セルをドラッグ&ドロップすることで、順次計算が可能になる。
下記にCalcの計算式を示す。
値を記入する
=A8+$B$1
=0.5*$B$4*A9^2+B$8+$B$3*A9
=C8+(E9+E8)*$B$1/2
=$B$4*A9+D$8
=$B$1*$B$4+E8
=$B$2
=$B$3
最初の行(第8行：A8～E8)に初期値を記入し、次の行(第9行：A9～E9)に計算式を記入する。
A9～E9をハイライト(選択)し、E9セルの右下にカーソルを合わせ、現れた十字マークをクリックし適当な行までドラッグ&ドロップすることにより、計算を実行できる。
B,D列(exact)は、直接解法による計算値、C,E列(numerical)は、数値解法による計算値を示す。直接解法では、時間の関数として計算するので、その時刻においても誤差はコンピュータの処理誤差だけであるが、数値解法においては、時刻Δt以降順次計算するので、誤差が蓄積される。この場合、100ステップまで計算したが、小数点以下7桁までは、誤差は無かった。