function RMSE=RMSE(predicted_time_series,observed_time_series)

RMSE=sqrt((sum(((predicted_time_series-observed_time_series).^2)))/(length(observed_time_series)));

end