TensorSolve

Added tensorsolve function.
In fact, I think it is good when axes is None, but when it isn't, the broadcastable pattern won't be good. There is to chosse, find a way to compute the broadcast pattern when axes isn't none, or just already return a broadcast patttern of False in that case. I'll go with the second case, as I don't see rapidly a way for the first, when the axes aren't constant.