A=[1 2 0 0;0 0 1 1;2 0 1 0;0 2 1 0;3 0 0 1];
c=[0.5;0.3;0.4;0.35];
si=A*c;
cfores=[0 40 80 0];
d=cfores*c;
p=cfores.*c';
anB=p(2)/c*100