maple计算问题
> restart:> omega22:=-k*(k^2*J^2*m+k*J*m*sqrt(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)-2*S^2*k*sqrt(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)-2*k*J*k*m^2+4*S^2*k*k*m-k*m^2*sqrt(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)+k^2*m^3)/(-m*J+S^2)/(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)^(1/2)/(-k*m+J*k+sqrt(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)):
> omega:=1/2/(J*m-S^2)*(k*m+J*k-sqrt(k^2*m^2-2*J*k*k*m+J^2*k^2+4*S^2*k*k)):
> simplify(omega-omega22);
0
> l:=244;p:=1.29;theta:=Pi/4;h:=10/1000;c:=-18;c:=3.8;c:=6;d:=-1.02;d:=0.99;d:=-0.21;M:=1.7;a:=0.4;dd:=0.028:M:=h*dd*0.9:dh:=dd+h*sin(theta);U:=y;m:=3*l*(M+M)*(2/Pi)^2;J:=(2/Pi)^2*l*(M+M)*(a)^2;mu:=2*m*0.04*sqrt(k/m):mu:=(2*J*0.06*sqrt(k/J)):a := (-1/2*p*U*dh*c-p*U*dh*c*theta-3/2*p*U*dh*c*theta^2)*l;b := (-1/2*p*U*dh^2*d-p*U*dh^2*d*theta-3/2*p*U*dh^2*d*theta^2)*l; k:=3*4*26000/l:k:=864:S:=3*l*(2/Pi)^2*M*dd/2;a := -1/2*p/U*dh*c*l;b := -1/2*p/U*dh^2*d*l;
> simplify(omega-omega22);
35.10262596
这是我的程序，蓝色部分是显示的计算结果。可以看到直接omega-omega22进行化简，结果是0，这说明omega和omega22是相等的。但是将各个参数的值都给出来以后，omega-omega22=35.10262595，两者竟然不相等了。有哪位高手给解释一下。
回复 楼主 zhanghb 的帖子
中间是不是出现了一些舍入的问题？
页:
[1]