pythonのプログラムで作った動画を再生できない

Question

pythonのプログラムで作った動画が再生されない

橋の地震波に対する減衰振動の様子を示すグラフを作成したのですが、完成した動画の再生ができません。どのような点を直せばいいか教えてほしいです。下のような動画ができ、本当は、振動するのですが再生ボタンを押しても動きませんでした。

作成した動画が再生できない(エラーメッセージはありません。)

``

import matplotlib.pyplot as plt

m = 360000 #質量
k = 57000000 #バネ定数
c = 450000 #減衰定数
dt = 0.01
calsize = 500

t = [0]*calsize
osc = [0]*calsize
t[0] = 0
t[1] = dt
osc[0] = 1
osc[1] = 1
for i in range(1,calsize-1):
 t[i+1]=t[i]+dt
 osc[i+1]=(2 - c/m*dt - k/m*dt**2)*osc[i] + (c/m*dt - 1)*osc[i-1]


l=4
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
from IPython.display import HTML
fig=plt.figure(figsize=[5,5])
plt.xlim(-1,1)
plt.ylim(0,6)

ims = []

for i in range(len(osc)):
 a =  osc[i]
 b = np.sqrt(l**2-osc[i]**2)
line = plt.plot([0, a], [0, b], "-",linewidth=1, color="black")
point = plt.plot(a, b, ".", markersize=10, color="red")
im = line + point
ims.append(im)
ani = animation.ArtistAnimation(fig, ims, interval=50)
plt.close()
ani.save("animation.gif", writer="pillow")
HTML(ani.to_jshtml())

### 試したこと

プログラムを何回か作り直した

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Comments

[etherbeg]

ソースコードは以下を参考にコードブロックとして挿入してください。
https://teratail.com/help/question-tips#questionTips3-5-1
https://teratail.com/help/question-tips#questionTips3-7

[meg_]

> 完成した動画の再生ができません
どのような方法で動画の再生を試みて、どのような結果になったのですか？

[jbpb0]

pythonのコードの一番最初の行のすぐ上に
```python
だけの行を追加してください
また、pythonのコードの一番最後の行のすぐ下に
```
だけの行を追加してください

または、
https://teratail.storage.googleapis.com/uploads/contributed_images/56957fe805d9d7befa7dba6a98676d2b.gif
を見て、そのようにしてみてください

現状、コードがとても読み辛いです
質問にコードを載せる際に上記をやってくれたら、他人がコードを読みやすくなり、コードの実行による現象確認もやりやすくなるので、回答されやすくなります

Answer 1

そのプログラムはJupyter notebookで実行することを前提としているので、Jupyter notebookで実行してください。
通常のターミナルで実行する場合は、plt.close()を削除して、HTML(ani.to_jshtml())をplt.show()に変更してください。

追記：
すみません、Google Colaboratoryでの実行でしたね。Google Colaboratoryで試してみましたが、ちゃんと動画が表示されました。

Answer 2

for i in range(len(osc)):
 a =  osc[i]
 b = np.sqrt(l**2-osc[i]**2)
line = plt.plot([0, a], [0, b], "-",linewidth=1, color="black")
point = plt.plot(a, b, ".", markersize=10, color="red")
im = line + point
ims.append(im)

↓ インデントを修正

for i in range(len(osc)):
 a =  osc[i]
 b = np.sqrt(l**2-osc[i]**2)
 line = plt.plot([0, a], [0, b], "-",linewidth=1, color="black")
 point = plt.plot(a, b, ".", markersize=10, color="red")
 im = line + point
 ims.append(im)

で、動きました