function [f_approx, f_coef] = Proj(f_0,f_0_select)
f_coef = f_0_select'*f_0;
f_approx = (f_coef'*f_0_select')';

