私は機能を持っています:
def g(R, r):
return (np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r) -
(1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R))
- 1))
関数は、delta v_B = delta v_H場所delta v_Bを設定することによって定義されます
np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r)
そしてdelta v_H_
1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R)) - 1
ということで、 と書きgましdelta v_b - delta v_Hた。
これが私の関数であり、以下で使用しているコードです。
import pylab
import numby as np
def g(R, r):
return (np.sqrt(2.0 * (R + r) / (r * R)) - (1 + np.sqrt(R)) / np.sqrt(R) -
np.sqrt(2.0 / (r * (1 + r))) * (1 - r) -
(1.0 / np.sqrt(R) - np.sqrt(2.0) * (1 - R) / np.sqrt(R * (1 + R))
- 1))
r = np.linspace(11.9, 16, 500000)
R = np.linspace(1, 20, 500000)
fig2 = pylab.figure()
ax2 = fig2.add_subplot(111)
ax2.plot(R, g(R, r), 'r')
pylab.xlabel('$R_1 = \\frac{r_C}{r_A}$')
pylab.ylabel('$R_2 = \\frac{r_B}{r_A}$')
pylab.xlim((0, 25))
pylab.ylim((0, 100))
pylab.show()
関数は、約で無限大に漸近し、約で11.94線と交差する必要がありますy = x15.58
どうすればそのようなプロットを作成できますか? 私はこれを行う方法に慣れておらず、そのような関数をプロットする方法も知りません。
私の定義はgasには不適切g(R, r)ですか? そうでない場合、どのように定義する必要がありますか?
