xmin = 0
xmax = 2
dx = 0.1

xlist = np.around(np.arange(xmin, xmax, dx), decimals=4)

ylist = ((8.831 * 10. ** 12.) / xlist ** 3.) ** 1/5.0

fig = plt.figure(figsize=(4, 3))
plt.plot(xlist, ylist, label=r'$y = \sqrt[5]{\frac{8.831 * 10^{12}}{x^3}} $')
plt.legend(loc='best')
plt.xlabel(r'$x$')
plt.ylabel(r'$y$')
# plt.savefig('image-1.pdf')