function f_approx = Proj_solu(f_z,f_0_select,f_T_select)
[~, f_coef] = Proj(f_z,f_0_select);
f_approx = (f_coef'*f_T_select')';

