%利用 fsolve 函数求解二元非线性方程组的数值解%计算参数m=641.46; %弯矩设计值，单位：千牛米fc=1.43e4; %混凝土抗压强度，单位：千牛每平方米r=0.5; %桩截面半径，单位：米A=0.785; %桩截面面积，单位：平方米fy=3.6e5; %箍筋抗拉强度，单位：千牛每平方米rs=0.44; %纵筋重心所在圆周的半径，单位：米pi=3.1416; %圆周率%构建内联函数，其中 x(1)表示受压区圆心角与 360 度的比值,x(2)表示全部纵筋面积（单位：平方米）fun=(x)x(1)*fc*A*(1-sin(2*pi*x(1)/(2*pi*x(1)-(1.25-3*x(1)*fy*x(2);2/3*fc*A*r*(sin(pi*x(1)3/pi+fy*x(2)*rs*(sin(pi*x(1)+sin(pi*(1.25-2*x(1)/pi-m;x0=0.2;0.004;%函数解初始值x=fsolve(fun,x0)
