GE(A[0..n-1,0..n])
// Input an n × (n + 1) matrix A[0 . . n − 1, 0 . . n] of real numbers
for i = 0 to n − 2
for j = i + 1 to n − 1
for k = i to n
A[j, k] = A[j, k] − A[i, k] ∗ A[j, i]/A[i, i]
How to compute for its running time formally?