diff --git a/math/Binomial Function.cpp b/math/Binomial Function.cpp
new file mode 100644
index 00000000..3a042efd
--- /dev/null
+++ b/math/Binomial Function.cpp	
@@ -0,0 +1,19 @@
+//Funcion factorial
+unsigned long int factorial (unsigned long int valor)
+{
+    if (valor<=1)
+        return valor=1;
+    else
+        return valor=valor*factorial(valor-1);
+}
+
+/*///////////////////////////////////////////////////////////////////////////////////////*/
+/*///////////////////////////////////////////////////////////////////////////////////////*/
+
+//Funcion binomial(requiere la funcion factorial)
+float numero_binomio (unsigned long int arriba,unsigned long int abajo)
+{
+    float respuesta;
+    respuesta = factorial(arriba) / (factorial(abajo) * factorial(arriba-abajo));
+    return respuesta;
+}
diff --git a/math/Gauss Jordan.cpp b/math/Gauss Jordan.cpp
new file mode 100644
index 00000000..5559fa76
--- /dev/null
+++ b/math/Gauss Jordan.cpp	
@@ -0,0 +1,32 @@
+
+void Gauss_Jordan (double m[][MAX],int fila,int columna)
+{
+    int i,j,t,k,l;
+    double divisores[MAX];
+    int fila_pivote=0;
+
+    for (t=0; t<fila; t++)
+    {
+
+        //Division unica
+        if(m[t][t]!=1)
+            for (l=columna-1; l>=t; l--)
+            {
+                m[t][l]=m[t][l]/m[t][t];
+            }
+
+        //Divisores
+         for (k=0; k<fila; k++)
+            divisores [k] = 1.0*m[k][t]/m[t][t];
+
+        for (i=0; i<fila_pivote; i++)
+            for (j=0; j<columna; j++)
+                m[i][j] = m[i][j] - m[t][j]*divisores[i];
+
+        for (i=1+fila_pivote; i<fila; i++)
+            for (j=0; j<columna; j++)
+                m[i][j] = m[i][j] - m[t][j]*divisores[i];
+
+        fila_pivote++;
+    }
+}