function elts=single_gaussian_matrix_elts(alph,alphp)
pprod=alph.*alphp;
pprod34=pprod.^(3/4);
ssum=alph+alphp;
ssum32=ssum.^(3/2);

pe_analytic=-2^(5/2)*pprod34./(sqrt(pi)*(alph+alphp));
ke_analytic=3*2^(3/2)*pprod.*pprod34./(ssum.*ssum32);
s1s=2^(3/2).*pprod34./ssum32;
e1s=pe_analytic+ke_analytic;
elts=[s1s e1s];
% 