function d=da(a,b)
c=length(a);
y=[a,b];
format rat
for n=1:c-1
for v=n+1:c
f=y(v,n)/y(n,n);
y(v,:)=y(v,:)-f*y(n,:);
end
end
a=y(1:c,1:c);
b=y(:,c+1);
d=inv(a)*b;
end
在這裡插入程式碼片

相關推薦