Trho=[1:0.1:50]
TSinf=[]

for i=1:length(Trho)

    rho=Trho(i)
    
    function y=f(Sinf)
        y=Sinf-50*exp(-(51-Sinf)/rho)
    endfunction
    
    R=dichoto(f,0,50,100)
    TSinf=[TSinf,R($)]
end

scf(1)
clf()

plot2d(50*Trho^(-1),TSinf)
