function [img_, Cn]=rec_big_silo(eps_simul, C_low, C_high, S_forw, S_inv)

img_=[];
Cn=[];
for i=1:min(size(eps_simul))
    i
    C_calc=S_forw*eps_simul(:,i);
    Cn1=(1./C_calc-1./C_low)./(1./C_high-1./C_low);
    img_=[img_ S_inv*Cn1];
    Cn=[Cn Cn1];
end
