% Program to Graph Joint & Marginal PDFs 

clc;
n = 10;
joint = 0;
marginal = 1;

x1 = zeros(n+1,1);
x2 = zeros(n+1,1);
f = zeros(n+1,n+1);
f1 = zeros(n+1,1);

for i  = (1:1:n+1)
for j  = (1:1:n+1)
    x1(i)=(i-1)/n;
    x2(j)=(j-1)/n;
    if j>i;
       f(i,j)=2;
    end
    f1(i)=2*(1-(i-1)/n);
end
end

if (joint == 1)
    surf(x1',x2,f);
    xlabel('x2');
    ylabel('x1');
end

if (marginal == 1)
    plot(x1,f1);
    xlabel('x1');
    ylabel('f1');
end