%主程序
function Planetruss44
%计算平面桁架的单元应力和节点位移
Definemodel %定义模型
Solvemodel %求解模型
Displayresult %显示结果
return
%子程序1
function Definemodel
global E m n YSZS jdls theta YS Element F L A
open('tbb.txt');
V=load('tbb.txt');
E=V(1);m=V(2);n=V(3);YSZS=V(4);jdls=V(5);
theta =zeros(1,m);YS= zeros(YSZS,3);Element=zeros(m,2);F=zeros(jdls,3);
L=zeros(1,m);A=zeros(1,m);
theta=[V(6):V(5+m)];YS=[V(6+m):V(6+m+3*YSZS-1)];
Element=[V(6+m+3*YSZS):V(6+m+3*YSZS+2*m-1)];
F=[V(6+m+3*YSZS+2*m):V(6+m+3*YSZS+2*m+jdls*3-1)];
L=[V(6+m+3*YSZS+2*m+jdls*3):V(6+m+3*YSZS+2*m+jdls*3+m-1)];
A=[V(6+m+3*YSZS+2*m+