Pythonでcan't convert complex to floatの解決方法

Question

前提・実現したいこと

重積分のグラフを描画したいのですがエラーが出て描画できません。

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

TypeError: can't convert complex to float

該当のソースコード

import numpy as np
from scipy import integrate
import sympy as sym
import math
import matplotlib.pyplot as plt
import cmath

E0 = 1
f = 3*10**9
c = 3*10**8
lum = c / f
k = 2*math.pi / lum

Nf = np.arange(0,20,0.5)

ds = 100
dp = 100

W = []
L = []

x = -1
y = 1

xs = 0
ys = 0

xp = 0
yp = 0

n = 1

rs = np.sqrt(ds**2 - 2*xs*x + x**2 - 2*ys*y + y**2)
rp = np.sqrt(dp**2 - 2*xp*x + x**2 - 2*yp*y + x**2)

Q = (n*rp / (abs(n)*abs(rp))) + (n*rs / (abs(n)*abs(rs))) / 2

intE = []
EP1 = Q*((np.exp(-1j*k*(rs+rp)))/(rs*rp))
for i in range(len(Nf)):
    w = np.sqrt((4*lum*ds*dp*Nf[i])/(ds+dp))
    l = w
    W.append(w)
    L.append(l)
a = sym.symbols('a')
b = sym.symbols('b')
for j in range(len(Nf)):
    inte = sym.integrate(EP1,(a,-W[j]/2,W[j]/2),(b,-L[j]/2,L[j]/2))
    intE.append(inte)

print(intE)

EP = []
for k in range(len(Nf)):
    EP2 = (1j*E0 / lum) * intE[k]
    EP.append(EP2)
print(EP)
print(EP[2])

fig = plt.figure()
ax = fig.add_subplot(111)
ax.plot(Nf,EP)
ax.set_xlabel('Nf')
ax.set_ylabel('E(P)')
plt.legend(loc='best')
plt.grid()

試したこと

複素数を扱っているためcmathモジュールを使ったり、配列を使って数値を格納したりしたのですが、特に変化は起こりませんでした。
EPをprintしたら中身はありました。

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

Python3.7.6
Spyder

Answer 1

matplotlib でプロットできないのは、EPがPythonが扱える数値ではなく、Sympyの式だからです。

type(EP[2])
# <class 'sympy.core.mul.Mul'>

type(1 + 1j*2)
# complex

Sympy ではプロット用のメソッドも実装されているのでそれを使えばよいようです。同じような質問がありましたのでリンクを貼っておきます

https://stackoverflow.com/questions/52574265/sympys-plot-parametric-says-typeerror-cant-convert-complex-to-float/52574358#52574358

Comments

[earths]

ご回答ありがとうございます。
ですが、sympyの中のメソッドで「EP」のプロットをしようとして、
----
p = sym.plot(EP)
p.show()
----
の文を挿入してみたのですが、
----
'NoneType' object is not iterable
----
という新たなエラーが発生してしまいました。
「EP」はお書きして頂いたように、NoneTypeではないですよね？
どういう事なのでしょうか。