Create copyed code.em

This code is copied from https://github.com/xccds/my-r-code/blob/master/EM.R