RAS

# initialize
data = load("data.txt");
a = data(1:4,1:6);
b = data(5:8,1:6);
c = ones(4,6);
c(1:4,6:6) = b(1:4,6:6);
c(4:4,1:6) = b(4:4,1:6);
cc = ones(4,6);
ep = 1;

#loop
while (ep > 0.001)
# R
temp = b(1:4,6:6) ./ a(1:4,6:6);
a = a .* temp;

#epsilon check


cc(1:4,6:6) = a(1:4,6:6);
cc(4:4,1:6) = a(4:4,1:6);

ep = abs(sum(sum(cc - c)));
if (ep < 0.001)
    printf ("a is final matrix \n");
    break;
endif

# C
temp = a(4:4,1:6) ./ b(4:4,1:6);
b = b .* temp;

#epsilon check
cc(1:4,6:6) = b(1:4,6:6);
cc(4:4,1:6) = b(4:4,1:6);

ep = abs(sum(sum(cc - c)));
endwhile