How to write a function which performs deconvolution to solve for x given y and h where y = h*x (all column vectors) using linsolve

I do know that the deconv() function can do it perfectly. However, I was wondering if anyone can write a function with linsolve() to perform deconvolution. I have tried to make a simple implementation but it seems to be not universal.

function x = mdeconv(h, y)
%Deconvolution
%
hT = toeplitz([h(1) zeros(1,length(y)-length(h))], [h zeros(1,length(y)-length(h)) ]);
%x = linsolve(hT, y);
x = y/hT;
end