% function err = rmse(X,EST)
function err = rmse(X,EST)
  err = sqrt(mean(sum((X-EST).^2)));
  
