Hi everyone, this forum seems very useful and this is my first post. I decided to explore Mathematica after using Python. I've generated a figure in Python and decided to replicate it in Mathematica. But I get completely different results. I'm fairly confident that the results from Python are correct.
Here is the Mathematica code I used to generate the figure.
I am also attaching the image from my Python code. I suspect that Mathematica is having an issue with Precision but I'm not sure exactly what's going on.

The code is essentially identical except that I used Python's fsolve and looped over an array. The function is identical.