After successfully defining estimate_gaussian, as an experiment (since earlier I had confusing error messages) I changed the last line to
return mu, var+1
And the error message incorrectly refers to mu instead of var.
AssertionError: Wrong value for mu. Expected: [1, 2, 3] got: [2. 3. 4.]