added code for calculating analytical SS probability, no output yet