function [sqrt_M,A_diag] = fcn_M_A(Z)
phi_x = fcn_phi_x(Z);
[sqrt_M, sqrt_M_plus] = M_np1(phi_x);
A_diag = VPFP_A_diag(sqrt_M, sqrt_M_plus,Z);