Matlab - Algoritmo en matlab

function [b,tp] = beta34v2(n)

c = 1.25;

b = []; tp = [];

for i = 1:n,                            % 

    v = rand;

	y=(-3/2)*ln(v);

    u = c * rand;                      % 

    while (u <  (( (2^(1/2)) * (exp(1/2))) / (3^(1/2)))* (y^(1/2)) * (exp(-y/3) )                    )   % 

        tp = [tp ; u1];                 % 

       v = rand;

	   y=(-3/2)*ln(v);

       u = c * rand;

    end;

    x = y;                             % 

    b = [b ; x];                        % 

    tp = [tp ; u1];                     % 

end;