code: def mdot(*args):
      
      if len(args) == 1:
          return args[0]
      elif len(args) == 2:
          return _mdot_r(args[0], args[1])
      else:
          return _mdot_r(args[:-1], args[-1])
docs: Multiply all the arguments using matrix product rules.
      The output is equivalent to multiplying the arguments one by one
      from left to right using dot().
      Precedence can be controlled by creating tuples of arguments,
      for instance mdot(a,((b,c),d)) multiplies a (a*((b*c)*d)).
      Note that this means the output of dot(a,b) and mdot(a,b) will differ if
      a or b is a pure tuple of numbers.