solve z*(1/(sqr(pi)*a^(3/2))*exp(-r/a))=-h^2/(2*m)*exp(-r/a)/(a^(7/2)*sqr(pi))+(-(k*c^2)/r)*(1/(sqr(pi)*a^(3/2))*exp(-r/a)) for z