E A NU t        已知條件
ki              stiffness 單元剛陣
K               zeros(n) assemble 組合總剛
T               eye(n) T(T,i,alpha)  K0=T*K*T'
k               邊界條件 取部分K
f               力邊界條件
u               u=k\f
U               位移邊界加u
F               F=K*U
ui              單元對應位移
sigmai          單元應力
sheardiagram    畫剪力圖
momentdiagram   畫彎矩圖