MATLAB command Daquan (below)