maxwell関数内で,受け取ったσを二次元にreshapeすること, 返り値を一次元で返す必要があることに注意します.ギリシャ文字はjupyter labで打っています. 読むときに楽で良いですね.成長せん断粘度ηs=σxy/κxyを見てましょう.plt.plot(sol.t,sol.y[1]/κ[0,1]) ...