If you have deep learning toolbox, then you can use the mse() function: https://www.mathworks.com/help/releases/R2020a/deeplearning/ref/dlarray.mse.html
estimated = [x_estimated y_estimated]; % x_predicted and y_predicted are column vectors
real = [x_real y_real]; % x and y are column vectors
mse_value = mse(estimated, real)
If you don't have the toolbox, then you can write your own function
function mse_value = mse(estimated, real)
mse_value = mean((estimated-real).^2, 'all');
end
