format long
error=data2;
n=length(error);
errsqr=error.^2;
sumerrsqr=sum(errsqr);
meanerrsqr=sumerrsqr/n;
rmserr=sqrt(meanerrsqr)
clear