In [40]:

    

import numpy as np
import matplotlib.pyplot as plt
import math

x = np.arange(0.97,1.03,0.001)

x0 = 1
v0 = 5

harm = (36*v0/x0**2)*(x-x0)**2 - v0
anharm = v0*((x0/x)**12-2*(x0/x)**6)

plt.plot(x,harm)
plt.plot(x,anharm)


print(x)


    

    




[0.97  0.971 0.972 0.973 0.974 0.975 0.976 0.977 0.978 0.979 0.98  0.981
 0.982 0.983 0.984 0.985 0.986 0.987 0.988 0.989 0.99  0.991 0.992 0.993
 0.994 0.995 0.996 0.997 0.998 0.999 1.    1.001 1.002 1.003 1.004 1.005
 1.006 1.007 1.008 1.009 1.01  1.011 1.012 1.013 1.014 1.015 1.016 1.017
 1.018 1.019 1.02  1.021 1.022 1.023 1.024 1.025 1.026 1.027 1.028 1.029
 1.03 ]






    

In [ ]: