for k=1:8
   err(k) = RMSE(k, 0, 0, 4); 
end

plot(1:8, err);