Hi kilik, I'm a stats phd student myself though I haven't done much with the game.
However, I can vouch for your equation. It follows straightforwardly from application of combinatorics and exchangeability. However, it needs to be noted that this probability is in fact the long-run probability of drawing two or more treasure maps simultaneously per hand drawn. Practical application has a couple of problems: firstly you need to apply this formula to the remaining cards in your deck, not the entire deck. Secondly, there's an issue with draws that force a reshuffle - you need to condition on the fact that in theory at least, you know exactly what the first few cards you draw will be.
But for example, if you write this probability as 1-q, then it's true to say that if your deck will last you k hands with possibly some cards spare left over, then the probability of drawing two treasure maps in those k draws is indeed 1-q^k. Wait no it isn't, due to the correlation effect. The average number of successes in those k draws is however (1-q)*k. If it's possible to only have one success (because you have 2 or 3 treasure maps) then this will equal the probability of drawing the right cards per deck cycle. David707's calculation is also correct, however. All of this can be validated with simulation code.
You shouldn't start your plots at 10 cards. Use of TM seriously without +cards generally requires lots of trashing. The real question is what having 1 haven or courtyard or similar does to these probabilities - but that's likely very difficult to calculate.