Thank you for your attention!
I have found the solution to the problem in CVX Forum.
Revised Code:
for e=1:E
I3(e,:)=r_ie(e,:)-(V(e)./(hmin_eu(e)+sum_e(e,:))).*(square_pos(norm(y-v(:,e))+lamb)-(norm(y_f-v(:,e))+lamb).^2);
end
Thank you for your attention!
I have found the solution to the problem in CVX Forum.
Revised Code:
for e=1:E
I3(e,:)=r_ie(e,:)-(V(e)./(hmin_eu(e)+sum_e(e,:))).*(square_pos(norm(y-v(:,e))+lamb)-(norm(y_f-v(:,e))+lamb).^2);
end