function phi_x_0 = fcn_phi_x_0(Z) global X; [Nz,~] = size(Z); [Nx,~] = size(X); %phi_x_0 = repmat(cos(X),[1,Nz]) + repmat(Z',[Nx,1])/10; star = floor(Nz/2); phi_x_0 = repmat(cos(X),[1,Nz]) + repmat([Z(1:star)',1+Z(star+1:Nz)'],[Nx,1])/10;