matlab解符号二元方程组和符号积分

符号方程组（分析力学）：

>> 1

ans =

     1

>> syms a b p k e w
>> clear
>> syms m1 m M xa xb g rou r
>> eqns = [m1*xa+m1*(rou/r)^2*(xa+xb)-m1*g==0,m*xb+M*xb+m1*(rou/r)^2*(xa+xb)-m*g==0]
 
eqns =
 
[m1*xa - g*m1 + (m1*rou^2*(xa + xb))/r^2 == 0, M*xb - g*m + m*xb + (m1*rou^2*(xa + xb))/r^2 == 0]
 
>> res=solve(eqns,xa,xb)

res = 

  包含以下字段的 struct:

    xa: [1×1 sym]
    xb: [1×1 sym]

>> res.xa
 
ans =
 
(g*(M*r^2 + m*r^2 - m*rou^2 + m1*rou^2))/(M*r^2 + M*rou^2 + m*r^2 + m*rou^2 + m1*rou^2)
 
>> 

符号积分（有限元节点内力分配求解）：

>> clear
>> syms kesi q1 q2
>> int((1+kesi)*kesi*((q2-q1)*kesi+(q1+q2)),kesi,-1,1)
 
ans =
 
(4*q2)/3
 
>> int((1+kesi)*kesi*((q2-q1)*kesi+(2*q1+q2)),kesi,-1,1)
 
ans =
 
(2*q1)/3 + (4*q2)/3
 
>> int((1-kesi)*(-kesi)*((q2-q1)*kesi+(q1+q2)),kesi,-1,1)
 
ans =
 
(4*q1)/3
 
>> int((1-kesi*kesi)*((q2-q1)*kesi+(q1+q2)),kesi,-1,1)
 
ans =
 
(4*q1)/3 + (4*q2)/3