function?x=followup（Ab）
n?=?rank（A）;
for（i=1:n）
????if（A（ii）==0）
????????disp（‘Error:?對角有元素為0！‘）;
????????return;
????end
end;

d?=?ones（n1）;
a?=?ones（n-11）;
c?=?ones（n-1）;

for（i=1:n-1）
????a（i1）=A（i+1i）;
????c（i1）=A（ii+1）;
????d（i1）=