Update matrixfactorizations.mac

Fixed a typo in the name of `diagonalize`, thanks @sangwinc!

Author: LukeLongworth

stack/maxima/contrib/matrixfactorizations.mac

@@ -69,7 +69,7 @@ get_Jordan_form(M):= block([jordan_info,J,P],
  * @param[matrix] M a diagonalizable matrix
  * @return[list] A list of two matrices: [P, D] such that D is diagonal and M = P . D . P^^-1. Returns empty list if M is not diagonalizable
  */
-ddiagonalize(M):= block([P,D],
+diagonalize(M):= block([P,D],
   if not(squarep(M)) then return([]),
   [P, D]: get_Jordan_form(M),
   if ev(zeromatrixp(M - transpose(M)),simp) then P:ev(transpose(apply('matrix, map(lambda([ex], ex/sqrt(ex.ex)), gramschmidt(args(transpose(P)))))), simp),