Lets assume I have an array of shape (36,100) and I would like to solve (using scipy.optimize.fsolve) each of these functions individually. I could do it with a for loop very easy but this is too slow for me… Is there a way I can do this in one single step?