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

