Woods-Saxon potentialをPythonで描画する

import numpy as np
from matplotlib import pyplot as plt
v0 = 14
A = 15
a = 0.5
r0 = 1.25
R = r0*A**(1/3)
x = np.linspace(0.0, 10.0)
y = -v0/(1+np.exp((x-R)/a))
plt.plot(x, y, "b-")
plt.title(r"Woods-Saxon potential")
plt.text(3, -v0*0.8, 
         r"$V(r)=\frac{-V_0}{1+{\rm exp}(r-R)/a}$"+
         "\n"+
         r"Parameters: A$={}$, $V_0={}$ MeV".format(A,v0)+
         "\n"+
         r"$a={}$ fm, $r_0={}$ fm, $R={:.1f} $fm".format(a,r0,R),
        fontsize=12)
plt.xlabel(r"$r$ [fm]")
plt.ylabel(r"$V(r)$ [MeV]")
plt.show()

オリジナルコード

http://nucleartalent.github.io/Course2ManyBodyMethods/doc/pub/intro/pdf/intro-print.pdf