def inv_prob(n):
        if n == 0:
                return 1
        if n > 365:
                return 0
        return inv_prob(n-1)*(1-(n-1.)/365)

def probarray():
        arr=[]
        for people in range(364):
                arr.append(1-inv_prob(people))
        return arr

def qrobarray():
        arr=[0]
        for people in range(399):
                arr.append(1-((1-arr[people])*(1-1./365)))
        return arr

plot(qrobarray())
plot(probarray(),'g.')
text(200,.5,'q(n)', color='blue')
text(35,.6,'p(n)', color='green')