Joint/Pooled Weibull Analysis

Does MATLAB have a built-in function for calculating the MLE joint/pooled shape parameter, with different scale parameters, for multiple samples?

I'm analyzing fatigue test data. I have a few sample groups of different loads, each of which has a few data points of cycles to failure or run-out (right censored). I'd like to calculate the shared Weibull shape, of the groups with different scales.