from scipy.interpolate import interp1d

N = [12, 13, 14, 15, 16, 
     17, 18, 19, 20, 21, 
     22, 24, 26, 28, 30, 
     34, 38, 43, 50, 60, 
     75, 100, 150, 300, 400]

Y = [0.245, 0.261, 0.277, 0.290, 0.296, 
     0.303, 0.309, 0.314, 0.322, 0.328, 
     0.331, 0.337, 0.346, 0.353, 0.359, 
     0.371, 0.384, 0.397, 0.409, 0.422, 
     0.435, 0.447, 0.460, 0.472, 0.480]

def fatorFormaLewis(N_dentes):
    fator = interp1d(N, Y)
    return fator(N_dentes)
