%  This solves the birthday problem assuming 365 days 
%  in the year, each equally likely of arising as a birthday
%  Here n is the number of people in a room and the project 
%  computes the probability that no two have the same birthday and 
% subtracts this from 1 to find the probability at least two people 
%  have the same birthday, and graphs this
%
a1=[ ];
a2=[ ];
for i=2:60
  c=1.;
  for j=1:i-1
    c=c*(365.-j)/365.;
  end;
  a1=[a1 i];
  a2=[a2 1.-c];
end;
plot(a1,a2),pause


