from math import*
def rongeur():
    t=0
    p=1
    while ......
        t=t+1
        p=3*exp(0.5*t)/(exp(0.5*t)+2)
    return(.....)