Added function compute E~