l2p3-spd

clc
clear
help p3.m

n=20;
p=1/2;

q=1-p;
X=zeros(2,n+1);
F=zeros(1,n+1);

for k=0:n
    x(1,k+1)=k;
    x(2,k+1)=(factorial(n)/(factorial(n-k)*factorial(k)))*p^k*q^(n-k);
end
%x';

F=cumsum(x(2,:));
[x',F']

subplot(2,1,1);
stem(x(1,:),x(2,:));
subplot(2,1,2);
stairs(x(1,:),F);