syms T(i);
n=1;
k01=0.445;
T0=60;
Le = 0.004;
Q = zeros(n+1,1);
T_L = [T(1);T(2)];
c1(i) = 2*k01*T0/(T(i)+T(i+1));
k1(i) = Bar1D2Node_Stiffness(c1(i),Le);
mm=zeros(n-1,n-1);
KK_L = blkdiag(k1(1),mm);
for i=2:n
KK_L = Bar1D2Node_Assembly(KK_L,k1(i),i,i+1);
end
KK_L(1,1)=5*10^8;
KK_L(n+1,n+1)=5*10^8;
Q(1)=5*10^8*100;
Q(n+1)= 5*10^8*40;
S = solve(KK_L*T_L-Q==zeros(2,1) );
出错 dadj (line 23)
S = solve(KK_L*T_L-Q==zeros(2,1) );
更改为这样也是错的[T(1),T(2)] = solve(KK_L*T_L-Q==zeros(2,1) );,求大佬指教!!!!!