diff --git a/openmp/dgemm_omp.c b/openmp/dgemm_omp.c
new file mode 100644
index 0000000..197bf1a
--- /dev/null
+++ b/openmp/dgemm_omp.c
@@ -0,0 +1,62 @@
+#include <stdio.h>
+#include <omp.h>
+
+#define N 2000
+
+double second()
+{
+        struct timeval tp;
+        int i;
+        i = gettimeofday(&tp,NULL);
+        return ( (double) tp.tv_sec + (double) tp.tv_usec * 1.e-6 );
+}
+
+int main(){
+
+	double ** A;
+	double ** B;
+	double ** C;
+	int mysize, myrank;
+	int chunk;
+	int i,j,k;
+
+	double t1,t2;
+
+	A = (double**)malloc(N*sizeof(double*));
+	B = (double**)malloc(N*sizeof(double*));
+	C = (double**)malloc(N*sizeof(double*));
+
+	for (i=0;i<N;i++){
+		A[i]=(double*)malloc(N*sizeof(double));
+		B[i]=(double*)malloc(N*sizeof(double));
+		C[i]=(double*)malloc(N*sizeof(double));
+	}
+
+	for (i=0;i<N;i++){
+		for (j=0;j<N;j++){
+		        A[i][j] = 1.0;
+		        B[i][j] = 2.0;
+		        C[i][j] = 0.0;
+		}
+	}
+
+	t1 = second();
+	#pragma omp parallel shared(A,B,C) private(i,j,k,myrank)
+	{
+	myrank=omp_get_thread_num();
+	mysize=omp_get_num_threads();
+	chunk=(N/mysize);
+	#pragma omp for schedule(static, chunk)
+		for (i=0;i<N;i++){
+			for (j=0;j<N;j++){
+				for (k=0;k<N;k++){
+					C[i][j]=C[i][j] + A[i][k]*B[k][j];
+				}
+			}
+		}
+	}
+	t2 = second();
+
+	printf("Time = %E [s], GF/s = %E \n",(t2-t1), (pow(N,3)/(t2-t1))/10E9 );
+
+}