Rick,
At the end of this article, the code have the same value for Qntl and y . I think it was not right ?

Here
u = rand("Uniform");
did not used in your code( random from Burr ) ,right?

I think the code should be like this ?
u = rand("Uniform");
y = Burr_Quantile(u, Theta, Alpha, Gamma);