Fitting a Weibull distribution using Scipy

Question

I am trying to recreate maximum likelihood distribution fitting, I can already do this in Matlab and R, but now I want to use scipy. In particular, I would like to estimate

Answer 1

the order of loc and scale is messed up in the code:

plt.plot(x, weib(x, scale, loc))

the scale parameter should come first.