function hsmatlab