Interpolation error

Id x(i), y(i), z(i) are constrained in an exponential cone for each value of i, use a for loop

for i=1:N
{x(i),y(i),z(i)} == exponential(1)
end

And follow the advice at CVXQUAD: How to use CVXQUAD's Pade Approximant instead of CVX's unreliable Successive Approximation for GP mode, log, exp, entr, rel_entr, kl_div, log_det, det_rootn, exponential cone. CVXQUAD's Quantum (Matrix) Entropy & Matrix Log related functions