I am trying to create this function in MATLAB. From a previous question, I have solved for a as "agrid2" and f(a) as "fx", whose dimensions are 600 by 1 and 600 by 2- these two vectors are absolutely correct:
m=600;
for i=1:m
mu= sum(agrid2.*fx(m,:));
gini= sum(fx(m,1)*(fx(m,2))*abs(fx(m,1)-fx(m,2)))./(2*mu);
end
However, my code is returning all blank answers- what did I do wrong? Any advice is greatly appreciated!