Нужно сделать генератор случайных чисел. "треугольный шум" через метод усечения

close all; clear all; clc;
 
n = 10000;
f = 2*(1-abs(-n/2+1:n/2) + n/2)/n;
% plot(f)
 
seq = zeros(1, n);
k = 1;
count = 0;
while(k<n)
    count = count + 1;
    x = 1+round(rand*(n-1));
    y = 1-2*rand;
    if abs(y)>f(x); continue; end
    seq(k) = y;
    k = k+1;
end
count
% plot(seq)
% figure; hist(seq)

%Генератор случайных чисел. на дельта y_noise. и N=1000.
%Random number generator. Треугольный шум.
function RNG_Triangle(delta_y, POINT_AMOUNT)
 
    j = 1;
    
    for i = 1:POINT_AMOUNT
        
        %a - по горизонтальной прмяой - аргумент.
        a = -delta_y + (delta_y-(-delta_y)).*rand();  
        %b - по вертикальной прмяой - функция.
        b = 0 + (1/delta_y-0).*rand();
        
        %fprintf('a = %d\n', a);
        %fprintf('b = %d\n', b);
        fprintf('Value_of_Function(a, delta_y)  = %d\n', Value_of_Function(a, delta_y));
 
        if Value_of_Function(a, delta_y) > b
            %точка внутри
            y_noise(j) = b;
            x_noise(j) = a;
            fprintf('y_noise(j) = %d\n', y_noise(j));
 
            j = j + 1;
        end
 
    end;
 
    figure;
    hist(x_noise);
end

function y_triangle = Value_of_Function(argument, delta_y)
 
    x1 = 0;
    y1 = 1/delta_y;
    y2 = 0;
 
    %fprintf('argument = %d\n', argument);
    %значение по горизонтали (аргумента) больше или меньше нуля 
    if argument > 0
        x2 = delta_y;       
    else
        x2 = -delta_y;       
    end;
    
    y_triangle = (argument - x1)*(y2 - y1)/(x2 - x1) + y1; 
     
    %fprintf('x2 = %d\n', x2);
    %fprintf('y_triangle = %d\n', y_triangle);
end