Solution.
First, formally express
f(z)=n∈Z∑ane2πinz,
Then f(z+1)=f(z) and
n∈Z∑ane2πinze−2πn=f(z+i)=eczf(z)=n∈Z∑ane2πi(n+2πic)z
Let 2πik=c, and suppose k∈Z. Then
n∈Z∑ane−2πne2πinz=n∈Z∑ane2πi(n+k)z=n∈Z∑an−ke2πinz,
so we require ane−2πn=an−k. If, say, k=−2, then if n=2m is even,
a2m=e−4πma2(m−1)=e−2πm(m+1)a0
and if n=2m+1 is odd,
a2m+1=e−4πma2(m−1)+1=e−2πm(m+1)a1
Thus, if we split the sum into even and odd terms, we see
f(z)=a0n∈Z∑e−2πn2−2πn+2πi(2n)z+a1n∈Z∑e−2πn2−2πn+2πi(2n+1)z=a0n∈Z∑e−2πn2−2πn+4πinz+a1n∈Z∑e−2πn2−2πn+4πinz+2πiz
satisfies f(z+1)=f(z) and f(z+i)=e−4πizf(z). Hence, if we pick a0 and a1 accordingly, we can also ensure that f(0)=0. Finally, the sum converges locally uniformly since e−2πn2 is summable, so such an f exists.