Changed doit() in HadamardProduct to be used by hadamard_product() (fixes #8557 )

fixes #8557.`hadamard_product()` now returns elementwise multiplication for `Matrix` objects and returns unevaluated hadamard product in case of `MatrixSymbol`.
```python
>>> a = Matrix([1, 1])
>>> b = Matrix([2, 2])
>>> h = hadamard_product(a, b)
>>> h
Matrix([
[2],
[2]])
>>> X = MatrixSymbol('X', 2, 1)
>>> Y = MatrixSymbol('Y', 2, 1)
>>> hadamard_product(X, Y)
X.*Y
>>> hadamard_product(a, X, b, Y)
Matrix([
[2],
[2]]).*X.*Y
```