import numpy as np

from msibi.optimize import R


def mie(eps, sig, m=12, n=6):
    """Mie pair potential. """
    return 4 * eps * [(sig / R) ** m - (sig / R) ** n]


def morse(D, alpha, r0):
    """Morse pair potential. """
    return D * [np.exp(-2 * alpha * (R - r0)) -
                2 * np.exp(-alpha * (R - r0))]


if __name__ == "__main__":
    import matplotlib.pyplot as plt

    y = mie(1.0, 1.0)
    plt.plot(R, y)
    plt.show()
@@ -126,17 +126,17 @@
 * eps * 
-%5B
+(
 (sig / R
@@ -158,17 +158,17 @@
  R) ** n
-%5D
+)
 %0A%0A%0Adef m
@@ -239,9 +239,9 @@
 D * 
-%5B
+(
 np.e
@@ -317,127 +317,6 @@
 r0))
-%5D%0A%0A%0Aif __name__ == %22__main__%22:%0A    import matplotlib.pyplot as plt%0A%0A    y = mie(1.0, 1.0)%0A    plt.plot(R, y)%0A    plt.show(
 )
+%0A

fix potentials