from math import *
import matplotlib
from matplotlib import pyplot as plt
a=52.9*10**-9
pi=3.1415
def phiOA(Z,a,x):
    return((1/(81*sqrt(3*pi)))*((Z)**1.5)*(27-18*((Z*x))+2*(Z*Z*x*x))*exp((-Z*x)))

P=[phiOA(12,a,i*0.001) for i in range (1,10000)]
A=[j*0.001 for j in range(1,10000)]
plt.plot(A,P)
matplotlib.pyplot.axis([0,1,-0.2,1])
plt.show()