import numpy as np
def distribution_invariante(lambd,mu):
    A=np.array([[-lambd,mu],[1,1]])
    B=np.array([0,1])
    pi=np.dot(np.linalg.inv(A),B)
    return pi