docs: def computePi(self,method='power'):
          """
          Calculate the steady state distribution using your preferred method and store it in the attribute `pi`. 
          By default uses the most robust method, 'power'. Other methods are 'eigen','linear', and 'krylov'
          
          Parameters
          ----------
          method : string, optional(default='power')
              The method for obtaining ``pi``. The possible options are 'power','eigen','linear','krylov'. 
          
          Example
          -------
          >>> P = np.array([[0.5,0.5],[0.6,0.4]])
          >>> mc = markovChain(P)
          >>> mc.computePi('power') 
          >>> print(mc.pi)
          [ 0.54545455  0.45454545]        
          
          See Also
          --------
          For details about the specific methods see
          :func:`powerMethod`,
          :func:`eigenMethod`,
          :func:`linearMethod`, and
          :func:`krylovMethod` .       
          """
          methodSet = ['power','eigen','linear','krylov']
          assert method in methodSet, "Incorrect method specified. Choose from %r" % methodSet
          method = method + 'Method'         
          return getattr(self,method)()