how to calculate the gamma distribution in matlab

main.m
a_values = [1 2 3];     % shape parameter
b_values = [2 2 4];     % scale parameter

x = linspace(0, 10, 1000);  % generate 1000 points between 0 and 10

% calculate PDF for each value of a and b
pdfs = zeros(length(a_values), length(x));
for i = 1:length(a_values)
    pdfs(i,:) = gampdf(x, a_values(i), b_values(i));
end

% plot the PDFs
figure;
hold on;
plot(x, pdfs(1,:), 'r-', 'LineWidth', 2);
plot(x, pdfs(2,:), 'g--', 'LineWidth', 2);
plot(x, pdfs(3,:), 'b-.', 'LineWidth', 2);
xlabel('x');
ylabel('PDF(x)');
legend('a=1, b=2', 'a=2, b=2', 'a=3, b=4');
title('Gamma PDFs for different a and b');
603 chars
22 lines

main.m
a_values = [1 2 3];     % shape parameter
b_values = [2 2 4];     % scale parameter
x_values = [2 4 6];     % upper endpoint

% calculate CDF for each value of a, b, and x
cdfs = zeros(length(a_values), length(x_values));
for i = 1:length(a_values)
    for j = 1:length(x_values)
        cdfs(i,j) = gamcdf(x_values(j), a_values(i), b_values(i));
    end
end

% plot the CDFs
figure;
hold on;
plot(x_values, cdfs(1,:), 'r-', 'LineWidth', 2);
plot(x_values, cdfs(2,:), 'g--', 'LineWidth', 2);
plot(x_values, cdfs(3,:), 'b-.', 'LineWidth', 2);
xlabel('x');
ylabel('CDF(x)');
legend('a=1, b=2', 'a=2, b=2', 'a=3, b=4');
title('Gamma CDFs for different a, b, and x');
664 chars
23 lines