Octave1

        % x = L3_segitseg("fel_1", [0,1],10000);
        case "fel_1"
            U = UMersenneTwisterRNG(n);
            for i = 1:n
                    X(i) = (sqrt(576*U(i)+49)-1)/6;
            end

PDF
    case "fel_1"
        f = zeros(1,n);
        for i = 1:n
            if(x(i)>=1 && x(i)<=4)
                f(i) = 2*(1/16)*(x(i)+1/6);
            else
                f(i) = 0;
            end
        end
https://gyazo.com/565fed98bf252dee921b3c960ed9364e