diff --git a/org/qmckl_blas.org b/org/qmckl_blas.org
index d42cc5f..7d204d6 100644
--- a/org/qmckl_blas.org
+++ b/org/qmckl_blas.org
@@ -979,7 +979,7 @@ integer function qmckl_dgemm_f(context, TransA, TransB, &
 
   do j=1,n
   do i=1,m
-     C(i,j) = alpha*C1(j,i) + beta*C(j,i)
+     C(i,j) = alpha*C1(j,i) + beta*C(i,j)
   end do
   end do