少し複雑な式が自分の力で書けません。。

Question

基本的で申し訳ないです。

mathematicaでなく、あえてpythonで書きたく、以下の式をplotしたいのですが
エラーが式の範囲外になる的なことが書かれているのですが、確認してもわかりません。

import numpy as np
import matplotlib.pyplot as plt
import math

#極限関数
a=1/math.sqrt(2)
x=np.linspace(-a, a, 100);
plt.figure(0)
plt.gca().set_aspect('equal', adjustable='box')

y=1/((math.pi)*(1-x**2)*(math.sqrt(1 - 2*x**2)))
plt.plot(x,  y)
plt.show()

Traceback (most recent call last):
  File "C:\Users\Rei Sato\Desktop\python�R�[�h\����.py", line 11, in <module>
    y=1/((math.pi)*(1-x**2)*(math.sqrt(1 - 2*x**2)))
TypeError: only size-1 arrays can be converted to Python scalars

申し訳ないです、ご指摘お願いします！！

Comments

[tachikoma]

math.sqrtをnp.sqrtで通ったりしませんかね。

[tachikoma]

あと、3/2乗とかも使えますよ。参考までに。

[tachikoma]

—-

Answer 1

大丈夫です。ほぼ書けてます。
ただ、numpy配列を一括で計算したいのでmath.sqrtではなくnp.sqrtを使う必要があります。
なお、提示されている式と実際のコードにすこし食い違いがあるようです。
コードを式に合わせて修正して得られたグラフは以下のようになりました。

Comments

[Fallout_18]

なるほど、意味がわかりました。
ありがとうございます！！

Answer 2

最近話題のアレ
http://ja.wolframalpha.com/input/?i=y%3D1%2F((pi)*(1-x**2)*(sqrt(1+-+x**2)))

Answer 3

ちょっと整理して、こんなのでどうですか。

y=1/(np.pi*
     ((1-x**2)**(3/2)))

Comments

[Fallout_18]

できました

[Fallout_18]

ありがとうございます！

Answer 4

単純に式の間違いじゃないのかな？（試してません）

y=1/((math.pi)*(1-x**2)*(np.sqrt(1 - x**2)))

ルートの中のxの2乗にかかっているかける2が余計かと。

違ってたらごめんなさい。

※修正：
よく見たら、math.sqrtじゃなくてnp.sqrtですね。エラーはそこででているようです。

Comments

[Fallout_18]

なるほど、意味がわかりました。
ありがとうございます！！