Theoretical and Numerical Study of SIR with variable susceptibility