Als ik scipy laat starten (achterliggend wordt de Powell hybrid method uit MINPACK gebruikt) met zoeken van een kleinere waarde van de stroom
def solve_for_i(t, initial_guess=0.0001):
# Define a wrapper function that fixes t
def equation_to_solve(i):
return f(i, t)
# Use fsolve to solve the equation
solution = fsolve(equation_to_solve, 0.000000001)
return solution[0]
Ziet het er al veel beter uit zonder dat ik iets aan de formules heb veranderd.