added formula 1 for cphi, made one big subplot for plot 4, added newest matlabcode version to attachment