function [select_i, error_max] = Max_Proj(f_0,f_0_select)
global delx; global delv;
[temp, ~] = Proj(f_0,f_0_select);
[error_max, select_i] = max(sum((f_0 - temp).^2)*delx*delv);
%[error_max, select_i] = max(max(abs(f_0 - temp)));

