import numpy as np
import matplotlib.pyplot as plt
# プロット範囲のxを用意
x = np.linspace(1,96)
# xに対応するyの値を用意, y:MW x:hour もとは秒なのを時間に変換する
# 一号機は0.5GWで熱出力はその3倍、その6.6%が最初の崩壊熱、-0.2乗で減衰

#y=150*3600**(-0.2)*x**(-0.2)

y=1.50*6.6*3600**(-0.2)*x**(-0.2)
fig, ax = plt.subplots()
ax.set_xlabel("hour")
ax.set_ylabel("decay Heat(MW)")
ax.set_title(r'Fukushima reactor number one')
ax.plot(x, y) 
plt.ylim(0)
plt.show()