I am trying to rewrite some codes from Mathematica into Python. Is there a way to describe Findroot function in Python 3.x? here is the code from wolfram mathematica!

```
[Mu]FindRoot[CDF[GammaDistribution[1/b, b], p] - .95 == 0, {b, 1, If[p == 1, 1.1, p]}][[1, 2]]]
#b is for beta
```