a1=860;a=2594900;
b1=322;c=286200;
b=991200;
c1=96;t=a1+b1+c1;
r1=t*a/(a+b+c);
r2=t*b/(a+b+c);
r3=t*c/(a+b+c);r=[r1 r2 r3];co=[35 16 10];w=[41 46 56];ww=sum(w);
for i=1:3
   ff(i)=co(i)*w(i)/ww;
   end
a2=35;
b2=43;

c2=70;
aa=[a b c];
bb=[a2 b2 c2];
e=aa./bb;
d1=3872;
d2=2083;
d3=2512;d=[d1 d2 d3];
s=2805;
g1=d1*s/(d1+d2+d3);
g2=d2*s/(d1+d2+d3);
g3=d3*s/(d1+d2+d3);g=[g1 g2 g3];
ss=1/s.*d;
