function phi_x_0 = fcn_phi_x_0(Z)
global X;
[Nz,~] = size(Z);
phi_x_0 = repmat(cos(X),[1,Nz]);
