%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第八章例8－1－1程序
% 直流電路的穩態計算
% 西安電子科技大學出版社出版  陳懷琛著 2007年5月
%
r1=2;r2=4;r3=4;r4=4;r5=2;r6=12;r7=12;		%為給定元件賦值
a11=1/(r1+r2)+1/r3+1/r7;a12=-1/r3;a13=1/(r1+r2);	%將系數矩陣各元素賦值
a21=1/r3;a22=-1/r3-1/(r4+r5)-1/r6; a23=0; 
us=input('us=');                  % 輸入解(1)的已知條件
A1=[a11,a12;a21,a22] 				% 列出系數矩陣A1
u=A1\[a13*us;0]							% u=[ub;uc]
ubc=u(1)-u(2),ude=u(2)*r5/(r4+r5),i=u(1)/r7	% 解(1)
ude=input('ude='),						% 輸入解(2)的已知條件
A2=[A1,-[a13;a23];0,r5/(r4+r5),0] 			% 列出系數矩陣A2
u=A2\[0;0;ude],							%u=[ub;uc;us]
ubc=u(1)-u(2),i=u(1)/r7,
us=(r1+r2)*(u(1)*a11+u(2)*a12) 	% 解(2)