一开始的就可以运行，结果如下
l =
(0.5*((2.0*R^2*cos(a)^4 + 2.0*b^2*cos(a)^4 - 2.0*r^2*cos(a)^4 - 1.0*f*cos(a)^2 + 2.0*f*cos(a)^2*sin(a)^2 - 1.0*f*cos(a)^2*sin(a)^4 - 2.0*p*t*cos(a)^2 - 2.0*R^2*cos(a)^4*sin(a)^2 + 2.0*b^2*cos(a)^4*sin(a)^2 + 2.0*r^2*cos(a)^4*sin(a)^2 + 4.0*p*t*cos(a)^2*sin(a)^2 - 2.0*p*t*cos(a)^2*sin(a)^4)*(4.0*p^2*t^2 - 4.0*R^4*cos(a)^4 - 8.0*b*z - 8.0*R^2*sin(a)^2 + 4.0*R^2*sin(a)^4 - 4.0*b^4*cos(a)^4 - 8.0*b^2*sin(a)^2 + 4.0*b^2*sin(a)^4 - 4.0*r^4*cos(a)^4 - 2.0*f^2*sin(a)^2 + f^2*sin(a)^4 + 8.0*r^2*sin(a)^2 - 4.0*r^2*sin(a)^4 + 4.0*R^2 + 4.0*b^2 + f^2 - 4.0*r^2 + 8.0*R^2*r^2*cos(a)^4 + 8.0*b^2*r^2*cos(a)^4 + 16.0*b^2*z^2*cos(a)^4 - 8.0*p^2*t^2*sin(a)^2 + 4.0*p^2*t^2*sin(a)^4 + 16.0*b*z*sin(a)^2 - 8.0*b*z*sin(a)^4 + 4.0*f*p*t - 8.0*R^2*b^2*cos(a)^4 - 8.0*f*p*t*sin(a)^2 + 4.0*f*p*t*sin(a)^4))^(1/2))/(2.0*R^2*cos(a)^4 + 2.0*b^2*cos(a)^4 - 2.0*r^2*cos(a)^4 - 1.0*f*cos(a)^2 + 2.0*f*cos(a)^2*sin(a)^2 - 1.0*f*cos(a)^2*sin(a)^4 - 2.0*p*t*cos(a)^2 - 2.0*R^2*cos(a)^4*sin(a)^2 + 2.0*b^2*cos(a)^4*sin(a)^2 + 2.0*r^2*cos(a)^4*sin(a)^2 + 4.0*p*t*cos(a)^2*sin(a)^2 - 2.0*p*t*cos(a)^2*sin(a)^4)
-(0.5*((2.0*R^2*cos(a)^4 + 2.0*b^2*cos(a)^4 - 2.0*r^2*cos(a)^4 - 1.0*f*cos(a)^2 + 2.0*f*cos(a)^2*sin(a)^2 - 1.0*f*cos(a)^2*sin(a)^4 - 2.0*p*t*cos(a)^2 - 2.0*R^2*cos(a)^4*sin(a)^2 + 2.0*b^2*cos(a)^4*sin(a)^2 + 2.0*r^2*cos(a)^4*sin(a)^2 + 4.0*p*t*cos(a)^2*sin(a)^2 - 2.0*p*t*cos(a)^2*sin(a)^4)*(4.0*p^2*t^2 - 4.0*R^4*cos(a)^4 - 8.0*b*z - 8.0*R^2*sin(a)^2 + 4.0*R^2*sin(a)^4 - 4.0*b^4*cos(a)^4 - 8.0*b^2*sin(a)^2 + 4.0*b^2*sin(a)^4 - 4.0*r^4*cos(a)^4 - 2.0*f^2*sin(a)^2 + f^2*sin(a)^4 + 8.0*r^2*sin(a)^2 - 4.0*r^2*sin(a)^4 + 4.0*R^2 + 4.0*b^2 + f^2 - 4.0*r^2 + 8.0*R^2*r^2*cos(a)^4 + 8.0*b^2*r^2*cos(a)^4 + 16.0*b^2*z^2*cos(a)^4 - 8.0*p^2*t^2*sin(a)^2 + 4.0*p^2*t^2*sin(a)^4 + 16.0*b*z*sin(a)^2 - 8.0*b*z*sin(a)^4 + 4.0*f*p*t - 8.0*R^2*b^2*cos(a)^4 - 8.0*f*p*t*sin(a)^2 + 4.0*f*p*t*sin(a)^4))^(1/2))/(2.0*R^2*cos(a)^4 + 2.0*b^2*cos(a)^4 - 2.0*r^2*cos(a)^4 - 1.0*f*cos(a)^2 + 2.0*f*cos(a)^2*sin(a)^2 - 1.0*f*cos(a)^2*sin(a)^4 - 2.0*p*t*cos(a)^2 - 2.0*R^2*cos(a)^4*sin(a)^2 + 2.0*b^2*cos(a)^4*sin(a)^2 + 2.0*r^2*cos(a)^4*sin(a)^2 + 4.0*p*t*cos(a)^2*sin(a)^2 - 2.0*p*t*cos(a)^2*sin(a)^4)
m =
(4.0*p^2*t^2 + 4.0*R^4*cos(a)^4 - 8.0*b*z - 12.0*R^2*sin(a)^2 + 12.0*R^2*sin(a)^4 - 4.0*R^2*sin(a)^6 - 12.0*b^4*cos(a)^4 - 12.0*b^2*sin(a)^2 + 12.0*b^2*sin(a)^4 - 4.0*b^2*sin(a)^6 + 4.0*r^4*cos(a)^4 - 3.0*f^2*sin(a)^2 + 3.0*f^2*sin(a)^4 - 1.0*f^2*sin(a)^6 + 12.0*r^2*sin(a)^2 - 12.0*r^2*sin(a)^4 + 4.0*r^2*sin(a)^6 + 4.0*R^2 + 4.0*b^2 + f^2 - 4.0*r^2 - 8.0*R^2*r^2*cos(a)^4 + 8.0*b^2*r^2*cos(a)^4 + 16.0*b^2*z^2*cos(a)^4 - 12.0*p^2*t^2*sin(a)^2 + 12.0*p^2*t^2*sin(a)^4 - 4.0*p^2*t^2*sin(a)^6 - 4.0*R^4*cos(a)^4*sin(a)^2 + 24.0*b*z*sin(a)^2 - 24.0*b*z*sin(a)^4 + 8.0*b*z*sin(a)^6 - 4.0*b^4*cos(a)^4*sin(a)^2 - 4.0*r^4*cos(a)^4*sin(a)^2 - 4.0*R^2*f*cos(a)^2 + 4.0*f*p*t + 4.0*b^2*f*cos(a)^2 + 4.0*f*r^2*cos(a)^2 - 8.0*R^2*b^2*cos(a)^4 - 8.0*R^2*p*t*cos(a)^2 + 8.0*b^2*p*t*cos(a)^2 + 24.0*R^2*b^2*cos(a)^4*sin(a)^2 + 8.0*p*r^2*t*cos(a)^2 + 8.0*R^2*r^2*cos(a)^4*sin(a)^2 - 24.0*b^2*r^2*cos(a)^4*sin(a)^2 - 16.0*b^2*z^2*cos(a)^4*sin(a)^2 + 8.0*R^2*f*cos(a)^2*sin(a)^2 - 4.0*R^2*f*cos(a)^2*sin(a)^4 - 12.0*f*p*t*sin(a)^2 + 12.0*f*p*t*sin(a)^4 - 4.0*f*p*t*sin(a)^6 - 8.0*b^2*f*cos(a)^2*sin(a)^2 + 4.0*b^2*f*cos(a)^2*sin(a)^4 - 8.0*f*r^2*cos(a)^2*sin(a)^2 + 4.0*f*r^2*cos(a)^2*sin(a)^4 + 16.0*R^2*p*t*cos(a)^2*sin(a)^2 - 8.0*R^2*p*t*cos(a)^2*sin(a)^4 - 16.0*b^2*p*t*cos(a)^2*sin(a)^2 + 8.0*b^2*p*t*cos(a)^2*sin(a)^4 - 16.0*p*r^2*t*cos(a)^2*sin(a)^2 + 8.0*p*r^2*t*cos(a)^2*sin(a)^4)/(16.0*b^3*cos(a)^4 - 8.0*b*f*cos(a)^2 + 16.0*b^3*cos(a)^4*sin(a)^2 + 16.0*R^2*b*cos(a)^4 - 16.0*b*r^2*cos(a)^4 + 16.0*b*f*cos(a)^2*sin(a)^2 - 8.0*b*f*cos(a)^2*sin(a)^4 - 16.0*b*p*t*cos(a)^2 - 16.0*R^2*b*cos(a)^4*sin(a)^2 + 16.0*b*r^2*cos(a)^4*sin(a)^2 + 32.0*b*p*t*cos(a)^2*sin(a)^2 - 16.0*b*p*t*cos(a)^2*sin(a)^4)
(4.0*p^2*t^2 + 4.0*R^4*cos(a)^4 - 8.0*b*z - 12.0*R^2*sin(a)^2 + 12.0*R^2*sin(a)^4 - 4.0*R^2*sin(a)^6 - 12.0*b^4*cos(a)^4 - 12.0*b^2*sin(a)^2 + 12.0*b^2*sin(a)^4 - 4.0*b^2*sin(a)^6 + 4.0*r^4*cos(a)^4 - 3.0*f^2*sin(a)^2 + 3.0*f^2*sin(a)^4 - 1.0*f^2*sin(a)^6 + 12.0*r^2*sin(a)^2 - 12.0*r^2*sin(a)^4 + 4.0*r^2*sin(a)^6 + 4.0*R^2 + 4.0*b^2 + f^2 - 4.0*r^2 - 8.0*R^2*r^2*cos(a)^4 + 8.0*b^2*r^2*cos(a)^4 + 16.0*b^2*z^2*cos(a)^4 - 12.0*p^2*t^2*sin(a)^2 + 12.0*p^2*t^2*sin(a)^4 - 4.0*p^2*t^2*sin(a)^6 - 4.0*R^4*cos(a)^4*sin(a)^2 + 24.0*b*z*sin(a)^2 - 24.0*b*z*sin(a)^4 + 8.0*b*z*sin(a)^6 - 4.0*b^4*cos(a)^4*sin(a)^2 - 4.0*r^4*cos(a)^4*sin(a)^2 - 4.0*R^2*f*cos(a)^2 + 4.0*f*p*t + 4.0*b^2*f*cos(a)^2 + 4.0*f*r^2*cos(a)^2 - 8.0*R^2*b^2*cos(a)^4 - 8.0*R^2*p*t*cos(a)^2 + 8.0*b^2*p*t*cos(a)^2 + 24.0*R^2*b^2*cos(a)^4*sin(a)^2 + 8.0*p*r^2*t*cos(a)^2 + 8.0*R^2*r^2*cos(a)^4*sin(a)^2 - 24.0*b^2*r^2*cos(a)^4*sin(a)^2 - 16.0*b^2*z^2*cos(a)^4*sin(a)^2 + 8.0*R^2*f*cos(a)^2*sin(a)^2 - 4.0*R^2*f*cos(a)^2*sin(a)^4 - 12.0*f*p*t*sin(a)^2 + 12.0*f*p*t*sin(a)^4 - 4.0*f*p*t*sin(a)^6 - 8.0*b^2*f*cos(a)^2*sin(a)^2 + 4.0*b^2*f*cos(a)^2*sin(a)^4 - 8.0*f*r^2*cos(a)^2*sin(a)^2 + 4.0*f*r^2*cos(a)^2*sin(a)^4 + 16.0*R^2*p*t*cos(a)^2*sin(a)^2 - 8.0*R^2*p*t*cos(a)^2*sin(a)^4 - 16.0*b^2*p*t*cos(a)^2*sin(a)^2 + 8.0*b^2*p*t*cos(a)^2*sin(a)^4 - 16.0*p*r^2*t*cos(a)^2*sin(a)^2 + 8.0*p*r^2*t*cos(a)^2*sin(a)^4)/(16.0*b^3*cos(a)^4 - 8.0*b*f*cos(a)^2 + 16.0*b^3*cos(a)^4*sin(a)^2 + 16.0*R^2*b*cos(a)^4 - 16.0*b*r^2*cos(a)^4 + 16.0*b*f*cos(a)^2*sin(a)^2 - 8.0*b*f*cos(a)^2*sin(a)^4 - 16.0*b*p*t*cos(a)^2 - 16.0*R^2*b*cos(a)^4*sin(a)^2 + 16.0*b*r^2*cos(a)^4*sin(a)^2 + 32.0*b*p*t*cos(a)^2*sin(a)^2 - 16.0*b*p*t*cos(a)^2*sin(a)^4)
扫一扫
7087
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
