Something like this
A = [-5000 1200 750 2000 1200; -200 150 60 40 50];
A = fliplr(A); % Do this in order to use the roots function
R1 = roots(A(1,:)); % roots of first row: 4 values from quartic polynomial
R2 = roots(A(2,:)); % roots of second row: ditto
% R = 1./(1 + x) so x = 1/R - 1
x1 = 1./R1 - 1;
x2 = 1./R2 - 1;
x1(abs(imag(x1))>10^-8)=[]; % delete complex values
x2(abs(imag(x2))>10^-8)=[]; % delete complex values
disp(x1)
disp(x2)
