function [mserror]=E2(y,fit)

n=length(y);
mserror=sqrt(sum((fit-y).^2)/n);