きっと明日もいい天気

メモ用ブログ

レポート用pythonスクリプト

グラフ

import numpy as np
import matplotlib.pyplot as plt

x = np.arange(0, 2000, 1)
plt.figure(figsize=(10, 10),dpi=100)
# 計算式
y = 0.981 * (x**2)+2.07*10**7 / (x**2)
z = 40191

plt.plot(x, y)
plt.plot(x, z)
plt.xlabel("V [m/s]")
plt.ylabel("T [N]")
plt.xlim(0, 250)
plt.ylim(0, 50000)
plt.grid()
http://plt.show ()

 

行列の累乗

 行列の固有値がすべて異なる場合
import numpy as np
import numpy.linalg as LA 

A=[[2,0],[1,1]]
v,w=LA.eig(A)
print(np.dot(np.dot(w,np.diag(np.power(v,4))),LA.inv(w)))