maple计算问题

zhanghb

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

无水1324

中间是不是出现了一些舍入的问题？