%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第七章例7－2－1程序
% 靜不定桿系的受力計算
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
P=3000;E=200e9;L=2
A1=200e-6;A2=300e-6; A3=400e-6;
a1=pi/3; a2 = pi/2; a3=3*pi/4;
L1=L/sin(a1); L2=L/sin(a2); L3=L/sin(a3);	% 計算桿長
K1=E*A1/L1; K2=E*A2/L2; K3=E*A3/L3;			% 計算剛度系數
D(1,1:5) = [cos(a1),cos(a2),cos(a3),0,0];   % 給系數矩陣按行賦值，避免語句太長
D(2,1:5) = [sin(a1),sin(a2),sin(a3),0,0];
D(3,1:5) = [1/K1,0,0,-cos(a1),-sin(a1)];
D(4,1:5) = [0,1/K2,0,-cos(a2),-sin(a2)];
D(5,1:5) = [0,0,1/K3,-cos(a3),-sin(a3)];
B = [P;0;0;0;0];
format long, X = D\B