function rho_0 = fcn_rho_0(Z) global X; [Nz,~] = size(Z); rho_0 = 2+repmat(sin(X),[1,Nz]);