H = hilb(3)                                    %生成三階希爾伯特數值矩陣
H = sym(H)                                   %將數值矩陣轉換成為符號矩陣
inv(H)                                        %求符號矩陣的逆矩陣
det(H)                                        %求符號矩陣的行列式
V = vpa(H,16)                                          %用16位精度轉換成數值矩陣
inv(V)                                                %對帶有誤差的矩陣進行求逆
det(V)                                                %求帶有誤差的矩陣的行列式
syms s;
H(1,1) = s                                            %定義基本符號變量并替換H(1,1)
Z = det(H)
sol = solve(Z)
H = subs(H,s,sol)
inv(H)