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