In [3]:

    

import numpy as np
import matplotlib.pyplot as plt
import matplotlib as mpl
mpl.rcParams["mathtext.fontset"] = "stix"
from matplotlib import rc
rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})
## for Palatino and other serif fonts use:
#rc('font',**{'family':'serif','serif':['Palatino']})
rc('text', usetex=True)


plt.rc('font', family='serif')
plt.rc('xtick', labelsize=30)
plt.rc('ytick', labelsize=30)



x = np.linspace(0.5, 20, 1000)

cou = -1/x
x = np.linspace(0, 20,1000)
mod = -1/np.sqrt(x**2+1)

plt.plot(x, cou,"--", color="red",label=r"Coulomb potential $-\frac{1}{r}$",linewidth=3)
plt.plot(x, mod, label=r"Model potential $-\frac{1}{ \sqrt{r^2+1}}$",linewidth=3)
plt.xlabel("$r \ (a. u.)$",size=40)
plt.ylabel("$V(x)  (a.u.)$",size=40)
plt.legend(loc=4,prop={'size':40},frameon=False)
plt.show()


    

In [ ]:

    




    

In [ ]: