r=exp(-theta^2/16) from theta=-10 to 10