r=1/(e^(50(theta-3*pi/4))+1)*1/(e^(-50(theta+3*pi/4))+1) from theta=-pi to pi