python code to calculate expected prob of finding twin stars