function phi_x = fcn_phi_x(z)
global X;
[Nx,~] = size(X); [Nz,~] = size(z); phi_x = zeros(Nx,Nz);
%phi_x = repmat(cos(X),1,Nz)+ repmat(z'/10,Nx,1);
for i = 1:Nz
if z(i)>3
    phi_x(:,i) = cos(X)+ (z(i)+1)/10;
end
if z(i)<= 3
    phi_x(:,i) = cos(X)+ z(i)/10;
end
end
