with(linalg);
a1[11] := 1;
a1[12] := 1;
a1[21] := 0;
a1[22] := 1;
A1:=matrix(2,2,[a1[11],a1[12],a1[21],a1[22]]);
b1[1]:= 223;
b1[2]:= 70.4;
B1:=matrix(2,1,[b1[1],b1[2]]);
sol := linsolve(A1,B1);a1 := 152.6;
b1 := 70.4;
ratio1:=a1/(a1+b1);
1-ratio1;a2[11] := 1;
a2[12] := 1;
a2[21] := 1;
a2[22] := 0;
A2:=matrix(2,2,[a2[11],a2[12],a2[21],a2[22]]);
b2[1]:= 223;
b2[2]:= 30.5;
B2:=matrix(2,1,[b2[1],b2[2]]);
sol := linsolve(A2,B2);a2 := 30.5;
b2 := 152.6;
ratio2:=a2/(a2+b2);
1-ratio2;(78.144+68.4+83.34)/3;100-%;