% 《工程線性代數（MATLAB版）》第七章例題7.4程序ea704
%
clear
syms q					% 規定符號變量
Q(1,2)=q; Q(2,3)=3/8*q-1/4; Q(3,1)=1; 	% 給非零元素賦值
Q(3,3)=0; 	% 給右下角元素Q(3,3)賦值后，矩陣中未賦值元素都自動置零
P=[2;1/4;0]			% 給P賦值
W=inv(eye(3)-Q)*P	% 用信號流圖求傳遞函數的公式
pretty(W(3))