function [ RE1, RE2 ] = RE12( IC0, UT, R1, R2, RC, Rt, Rg, Beta, Au, RE )
alfa = Beta/(Beta+1);
A = RePlus(RC,Rt);
B = RePlus(R1,R2);
rE= UT/IC0;
rB= (B+1)*rE;


a0= Au*Rg*B*rE + Au*Rg*rB*rE + Au*B*rB*rE + alfa*A*B*rB;
a1= Au*Rg*B + Au*Rg*rB + Au*B*rB + Au*Rg*rE + Au*B*rE + alfa*A*B;
a2= Au*(Rg+B);

p = [a2 a1 a0];
RE1 = roots(p);
RE2 = RE - RE1;
disp(['RE1= ', num2str(RE1')]);
disp(['RE2= ', num2str(RE2')]);

end