I'd suggest you reading the material available for R package RSToolbox. A google seach on EVI + Landsat will provide you with tons of other references.
at http://bleutner.github.io/RStoolbox/rstbx-docu/spectralIndices.html
Actually, right now i have Landsat 7 surface reflectance (SR) products along with EVI (GeoTiff). The EVI value ranges upto 0 - 9000 while i recalculated the EVI from the SR bands and clipped to my study area, than the value ranger from -200 to 178, after applying thresholding in ENVI.
I'm not getting why not the value rnages within -1 to 1.
You can try 2- band EVI using this formula EVI2=2.5*((NIR-Red)/(NIR+2.4*Red+1)). Map Algebra in ArcGIS can perform this equation. You can also visit this site and read the article