符号运算求助

lyf1225

下面程序的最终目的是求解方程 得出beta值
syms d sigmas beta real  
miu1=284.2; sigma1=24.5; miu2=0.03; sigma2=0.003;  
g=500*pi*d^2*sigmas-9.8e+5;  % 函数 %
pd1=-diff(g,'sigmas')*sigma1; % 函数对各变量求偏导 %
pd2=-diff(g,'d')*sigma2;
sigmasx=miu1;dx=miu2;  % 初值取随机变量均值 %
subs(pd1,sigmas,sigmasx);   
subs(pd1,d,dx);
subs(pd2,sigmas,sigmasx);
subs(pd2,d,dx);
cos1=pd1/(pd1^2+pd2^2)^0.5;
cos2=pd2/(pd1^2+pd2^2)^0.5;
sigmasxp=miu1+beta*sigma1*cos1;
dxp=miu2+beta*sigma2*cos2;
solve('500*pi*dxp^2*sigmasxp-9.8e+5=0');