Added functions for computing stationary distribution, fixation probability and transition matrix