?=

1/pi*integral(-1/2e^(1/4((piI – 4Itheta – 4z)cos(1/2theta) + (pi + 4Iz – 4theta)sin(1/2theta))/(Isin(1/2theta) – cos(1/2theta)))Ei(-(Isin(1/2theta) + cos(1/2theta))/(Isin(1/2theta) – cos(1/2theta)))sin(1/4pi + theta) – 1/2e^(1/4((piI – 4Itheta + 4z)cos(1/2theta) – (pi – 4Iz – 4theta)sin(1/2theta))/(Isin(1/2theta) + cos(1/2theta)))Ei(-(Isin(1/2theta) – cos(1/2theta))/(Isin(1/2theta) + cos(1/2theta)))sin(1/4pi + theta) – 1/2e^(1/4((piI – 4Itheta – 4z)cos(1/2theta) + (pi + 4Iz – 4theta)sin(1/2theta))/(Isin(1/2theta) – cos(1/2theta)))gamma((Isin(1/2theta) + cos(1/2theta))/(Isin(1/2theta) – cos(1/2theta)), z)sin(1/4pi + theta) – 1/2e^(1/4((piI – 4Itheta + 4z)cos(1/2theta) – (pi – 4Iz – 4theta)sin(1/2theta))/(Isin(1/2theta) + cos(1/2theta)))gamma((Isin(1/2theta) – cos(1/2theta))/(Isin(1/2theta) + cos(1/2theta)), z)sin(1/4pi + theta),theta,-pi,pi)