A=[99.23;100.05;91;107.71;104.1];
B=[3441 3441 301720.5;68750 1068750 0;170040 13085020 0;229350 229350 5729350;34194000 0 0];
N=[300000;1000000;13000000;5500000;32800000];
E=[-0.00302;-0.00261;-0.00208];



[c3,r3]=size(A);
[c4,r4]=size(B);

Valorantigo(1:c3,1)=A(1:c3,1).* N(1:c3,1) ./100;


x=sym ('x',[1 c3]);
x=transpose(x);

 for i=1:c3
    Valor(i,1)=symfun(0,x);
    for j=1:r4
        Valor(i,1)=symfun((Valor(i,1)/(1+E(j,1)+x(i,1))^j)+((B(i,j)/((1+E(j,1)+x(i,1))^j))),x);
    end
 end


eqn=Valor(1:c3,1)==Valorantigo(1:c3,1);

[x1,x2,x3,x4, x5, param, cond] = solve(eqn, x, 'ReturnConditions', true);