wblfit(vX, a, b) uses Newton-Raphson method to calculate Maximum Likelihood Estimates (MLE) of a (scale) and b (shape) for Weibull distributed dataset vX.
The probability density function of Weibull distribution is
where a > 0 is scale and b > 0 is shape, with and .