RMS error

n=500;
for j=1:n;
    rmsn(j)=sqrt((Reference(j)-b(j))^2);
end

rms=sum(rmsn)/n;