load soi88.dat
soi88 = [soi88(1:5,1:2);
         soi88(1:5,3:4)];

for i=1:2;
  x  = soi88((i-1)*5+1:i*5,:);
  r  = [x(2:5,2)./x(2:5,1);
        (x(1,2)-sum(x(2:5,2)))/(x(1,1)-sum(x(2:5,1)))];
  R(:,i) = r;
  I(:,i) = x(:,2);
  A(:,i) = x(:,1);
end;

