This requires Symbolic toolbox -
%Defining symbolic variables
syms x n
%summation using symsum()
out = symsum(ExpressionToBeSummed, n, 0, Inf)
You will obtain the sum in terms of the Hypergeometric function - hypergeom
%define final expression
y = 2*exp(-x^2)*out/sqrt(sym(pi))