Estoy haciendo pruebas de bondad de ajuste de distribuciones de probabilidad para datos aleatorios dados. Para complementar creo gráficos de este estilo: 

[![Distribución Normal Estandar][1]][1]


Pero me gustaría poder indicar la media y la varianza de la distribución usando sus letras griegas (mu y sigma cuadrado). ¿Es posible hacer eso en Matplotlib para obtener un gráfico con etiquetas como el siguiente? 

[![Imagen de distribuciones normales de wikipedia][2]][2]

Estoy usando matplotlib y scipy, mi código es el siguiente: 


    import matplotlib.pyplot as plt
    import numpy as np
    import scipy.stats as st
    mu, sigma = 0, 1.0 # media y desvio estandar
    datos = np.random.normal(mu, sigma, 10000) #creando muestra de datos
    x=np.linspace(-4,4,num=1000)
    y=st.norm.pdf(x,0,1)

    plt.plot(x,y,'r')
    # histograma de distribución normal.
    plt.hist(datos, 30,color='c',histtype='bar',cumulative=False,edgecolor='black', linewidth=0.7,density=True,label="Datos")
    plt.plot(x,y,'r--',label="PDF distribucion Normal ")
    plt.ylim(0,0.5)
    plt.ylabel('frequencia')
    plt.xlabel('valores')
    plt.title('Histograma - Variable aleatoria normal')
    plt.legend(loc="upper right") 
    plt.show()



  [1]: https://i.sstatic.net/5rQzI.png
  [2]: https://i.sstatic.net/hA624.png