%Generating Triangular RV for (C=0; B=3; A=5)
C=0; B=3; A=5;
Nu = 10000; Nh = 100;
u = rand(1, Nu); % u ~ RV of Uniform(0,1)
% Evaluate y via IDF
y = (u <= B/(A-C)).*sqrt((A-C)*(B-C)*u) + (u > B/(A-C)).*(A-sqrt((1-u)*(A-C)*(A-B)))
hist(y, Nh) % show histogram of y