グラフプロットを行いたい。

Question

前提・実現したいこと

まだ、慣れておらず何をどうしたらいいかわからない
グラフプロットしたいが、途中のエラーでプロットの項まで進めない。

発生している問題・エラーメッセージ

NameError                                 Traceback (most recent call last)
<ipython-input-1-18ad85d14bbd> in <module>
      5 
      6 
----> 7 y = integrate.quad(lambda z: special.kv(5/3,z),x,np.inf)
      8 F=x*y
      9 plot((F),(x,0,10),title='test',xlabel='x',ylabel='F')

NameError: name 'x' is not defined

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt
import scipy.integrate as integrate
import scipy.special as special


y = integrate.quad(lambda z: special.kv(5/3,z),x,np.inf)
F=x*y
plot((F),(x,0,10),title='test',xlabel='x',ylabel='F')

試したこと

x=symbols('x')を行える。
from sympy import *
を入れるとエラー
AttributeError: 'function' object has no attribute 'quad
が出てしまい積分自体ができない。

補足情報（FW/ツールのバージョンなど）

windowsのWSLでanacondaを通して、jupyter notebookを使用。

Comments

[jeanbiego]

積分範囲（x,inf）を変えて結果をプロットしたいのだと思うのですが、積分したいのは
∫ (special.kv(5/3,z)) dz ということですか？　更にその結果とxの積を取りたいのでしょうか。
具体的な数式を、手書きでもいいので載せてもらったほうが回答がつくと思います。

[meg_]

｢該当のソースコード｣のコードにはxが定義されてませんので、｢NameError: name 'x' is not defined｣のエラーが出るのは当然かと思いますが、どこか別の場所で定義済なのでしょうか？

[hataki7]

ご指摘ありがとうございます。
x=変数
y=∫ (special.kv(5/3,z)) dz 　範囲[x~∞]
F=ｘ*ｙ　　
最終的にはグラフ横軸　x　縦軸　F での描画を行いたいのでｘは定義していないですね。

Answer 1

たとえば、下記のような感じでしょうか？

from scipy.special import kv
import matplotlib.pyplot as plt
from scipy.integrate import quad
import numpy as np

xs = np.linspace(0, 10, 11)
f = lambda z: kv(5/3,z)
fs = [quad(f,x,np.inf)[0]*x for x in xs]

fig = plt.figure()
ax = fig.add_subplot(1,1,1)
ax.plot(xs,fs)
ax.set_title('test')
ax.set_xlabel('x')
ax.set_ylabel('F')
plt.show()

Comments

[hataki7]

流れとしてはそうですね。
F=ｘ*（∫ (special.kv(5/3,z)) dz 　範囲[x~∞]）
なので、関数としては少し違う感じがします。

[jeanbiego]

このコードは、kv(5/3,z)をx~∞まで積分した結果を、x = 0,1,...,10の11個得るものです。
具体的にはどの辺が違っていましたか？　あと、計算結果とは合っていましたか？

[hataki7]

具体的には、kv(5/3,z)をx~∞まで積分した結果[y]　にｘを掛けたもの　[F]とｘとの関係グラフを書きたいのですが、ブロードキャストのエラーが出てしまいます。

[jeanbiego]

xとの積をとったものに変更してみました。

[hataki7]

これです。おぉ。すごい。
こういう書き方なんですね。ありがとうございます。
pythonにまだ自分が慣れていないのでエラーメッセージでてんてこ舞いだったんですが、jeanbiegoさんのおかげ勉強になりました。