def dep(M):
    x=6
    y= 4.8*x+252.1
    while ... :
        x=...
        y=...
    return x+2013