>> A = randn(4,3) A = 0.5377 0.3188 3.5784 1.8339 -1.3077 2.7694 -2.2588 -0.4336 -1.3499 0.8622 0.3426 3.0349 >> [m,n] = size(A) m = 4 n = 3 >> A = randn(4) A = 0.7254 -0.1241 0.6715 0.4889 -0.0631 1.4897 -1.2075 1.0347 0.7147 1.4090 0.7172 0.7269 -0.2050 1.4172 1.6302 -0.3034 >> B = randn(4) B = 0.2939 -1.0689 0.3252 -0.1022 -0.7873 -0.8095 -0.7549 -0.2414 0.8884 -2.9443 1.3703 0.3192 -1.1471 1.4384 -1.7115 0.3129 >> X = rh7(A,B) X = 17.1736 1.2079 -11.7209 5.3803 7.8162 2.3658 -6.8693 3.3235 -6.3173 -0.7627 4.9633 -2.5238 -17.8394 -3.7642 14.2929 -7.1062 >> A*X*A-B ans = 1.0000 0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 -0.0000 0.0000 1.0000 -0.0000 0.0000 0.0000 -0.0000 1.0000 >> A*X*A-B-eye(4) ans = 1.0e-14 * 0 0.1110 0.1943 0.0180 0.0999 -0.1110 0.3997 0.0139 -0.1110 0.1776 0.3331 -0.1665 0.1776 0.4885 -0.1332 0.2442 >> p = [1 1 2 2 4] p = 1 1 2 2 4 >> roots(p) ans = 0.5000 + 1.3229i 0.5000 - 1.3229i -1.0000 + 1.0000i -1.0000 - 1.0000i >>