python 计算物理_python与计算物理：正负离子势能随距离的变化

正负离子势能随距离的变化曲线，如下图：

直接上代码：

#!/usr/bin/env python

#coding=utf-8

'''

Created on 2013年10月12日

@author: yanghl

正负离子势能随距离的变化

'''

import numpy as np

import matplotlib.pylab as plt

x = np.linspace(1e-10, 1.0, 1000)

E_N = -1.436/x+5.86e-6/x**9

Emin = min(E_N)

x1 = np.linspace(1e-10,0.279,1000)

Emin = [Emin for i in x1]

plt.plot(x,-1.436/x+5.86e-6/x**9,x,-1.436/x,':',x,5.86e-6/x**9,'-.',linewidth=3)

plt.legend(['$E_N$','$E_A$','$E_R$'])

plt.plot(x1,Emin,'--',linewidth=2.5)

E0 = np.linspace(-6,-4.575,1000)

x0 = [0.279 for i in E0]

plt.plot(x0,E0,'--',linewidth=2.5)

plt