function f_0 = ID_EQ(Z)
phi_x_0 = fcn_phi_x_0(Z);
rho_0 = fcn_rho_0(Z);
f_0 = EQ(rho_0,phi_x_0);