EM demo
http://lcn.epfl.ch/tutorial/english/gaussian/html/index.html