python - Python のアンダーフローの問題

Question

次の関数を描画する必要があります。これを私は理論データと呼んでいます

a=z-(L/2) b=z+(L/2) 以下の場合

• N=75
• 私=1
• L=0.16
• R=0.013

z は -0.011 から 0.011 の間です。私の推測では、アンダーフローの問題があると思います。コードは次のとおりです。

import numpy as np
import matplotlib.pyplot as plt

class Solenoid():

def __init__(self,N,I,L,R):
    self.N=N
    self.I=I
    self.L=L
    self.R=R

def field(self,z):

    mu=4*(np.pi)*(10 ** -7)

    coef=(mu*self.N*self.I)/(2*self.L)

    a=z+(self.L/2)
    b=z-(self.L/2)

    sum_a=a/((self.R ** 2 + a ** 2) ** 0.5)

    sum_b=b/((self.R ** 2 + b ** 2) ** 0.5)

    return (coef*(sum_a - sum_b))

sol=Solenoid(75,1,0.160,0.013)
x=np.linspace(-0.011, 0.011, 10000)
y=[]

for point in x:
    y.append(sol.field(point))

plt.plot(x,y)
plt.show()  

出力は

これは明らかに間違っています。予想される出力は、実験データに沿った多かれ少なかれ滑らかな線です。

UPDATE 新しい出力は