var('r')
x=2*r/a
probability=0.5*integrate(x**2*e**(-x),x,0,2)