p=[...], G=[...], v=[...], M1=[...], N1=[...], N2=[...], z=[...], h=[...], v=[...], R=[...] eq=表达式 syms x r theta u1=int(int(eq,theta,-1,1),r,0,R)
