function answer = RMSE(matrix1, matrix2, n)
matrix2 = matrix2(n+1:length(matrix2)-n);
answer = norm((matrix1-matrix2'))/sqrt(length(matrix1));
% sum = 0;
% for i=1:length(matrix1)
%     sum = sum + ((matrix1(i) - matrix2(i + n+1))^2);
% end
% answer = sqrt(sum/length(matrix1));