y0=-1/sqrt(3); y1=1/sqrt(3); gy0=(1+y0)*exp(2*(1+y0)); gy1=(1+y1)*exp(2*(1+y1)); integral=gy0+gy1