calcular bimodal (dos modas)

[biancai]

Hola a todos, estoy haciendo un programa que calcula la moda, mediana, varianza y desviación estándar de un arreglo, hasta aquí todo bien. El problema que tengo es que no he podido hacer que muestre cuando hay mas de dos modas (bimodal), pues hasta ahora solo muestra una moda. Si alguien puede ayudarme con esto muchas gracias.
el codigo...

private void button2_Click(object sender, EventArgs e)
{
int[] aux = new int[arreglo.Length];

for (int i = 0; i < arreglo.Length; i++)
{
for (int j = 0; j < arreglo.Length; j++)
{
if (arreglo == arreglo[j])
aux++;
}
}
int mayor = int.MinValue;
int index = 0;

for (int i = 0; i < aux.Length; i++)
{
if (aux > mayor)
{
mayor = aux;
index = i;
}
textBox2.Text = Convert.ToString(arreglo[index]); // aqui muestra cual es la moda, pero si existen dos modas, solo muestra 1, por ejemplo si los numeros 5 y 8 se repiten 10 veces, los dos son moda.
textBox3.Text = Convert.ToString(mayor);
}

}