how do I find 7*3 roots using bisection , where roots are supposed to in s_range. I need complete code and it has to be Bisection method

 

import numpy as np import scipy.optimize as optimize A = 2e-5 nsample = 6 B = np.array ( [0.000215, 0.001424, 0.001274, 0.002568, 0.000748, 0.000273]) L = np.array([0.0124, 0.0305, 0.111, 0.301, 1.14, 3.011) rhos = [0, -.01950, .0052] rho=A*s+np.sum (B/s+L) s_range= [[-L[0], np.inf], # S0 # s1 [-L[1], -L[0]], [-L[2], -L[1]], #s2 [-L[3], -L[2]], #s3 [-L[4], -L[3]], #s4 [-L[5], -L[4]], #s5 [-np.inf, -L[5]] # s6 ]