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