scipy クックブックのコードを変更しようとしています
scipy クックブックは、いくつかの状態方程式を解いてから出力します。
wsol = odeint(two_springs.vectorfield, w0, t, args=(p,),
atol=abserr, rtol=relerr)
for t1, w1 in zip(t, wsol):
print t1, w1[0], w1[1], w1[2], w1[3]
scipy クックブックでは、プログラムを実行し、印刷物をファイルに保存する必要があります。
program.py > file.dat
次に、別のプロット関数でそのファイルからデータをロードする必要があります
t, x1, xy, x2, y2 = loadtxt('two_springs.dat', unpack=True)
作図前
plot(t, x1, 'b', linewidth=lw)
plot(t, x2, 'g', linewidth=lw)
t、x1、x2とプロットを設定したいだけです
誰でもその方法を教えてもらえますか?
これは私が変更しようとしているものです: