-1

私は機能を持っています:

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)ですか? そうでない場合、どのように定義する必要がありますか?

ここに画像の説明を入力

4

1 に答える 1