Rational simplification bug ("Quotient by a polynomial of higher degree" error message)
- Subject: Rational simplification bug ("Quotient by a polynomial of higher degree" error message)
- From: andre maute
- Date: Thu, 29 Nov 2012 10:45:34 +0100
For convenience a self-contained maxima batch file.
-------------- next part --------------
display2d : false;
(2*(bet + gamm)*(G*M_1/sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) + G*M_2/sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2))/c^2 - (vx_0^2 + vy_0^2 + vz_0^2)*gamm/c^2 - 1)*((x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2)) + 1/2*(4*gamm + 3)*G*M_1*a_x1/(sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)*c^2) + 1/2*(4*gamm + 3)*G*M_2*a_x2/(sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)*c^2) + (x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2) - 1/2*(x_0 - x_2)*(2*(gamm + 1)*(vx_2^2 + vy_2^2 + vz_2^2) - 4*(gamm + 1)*(vx_0*vx_2 + vy_0*vy_2 + vz_0*vz_2) - 2*(2*bet - 1)*U_2 - (z_0 - z_2)*a_z2 - (y_0 - y_2)*a_y2 - (x_0 - x_2)*a_x2 - 3*((z_0 - z_2)*vz_2 + (y_0 - y_2)*vy_2 + (x_0 - x_2)*vx_2)^2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2))*G*M_2/(((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2)*c^2) - 1/2*(x_0 - x_1)*(2*(gamm + 1)*(vx_1^2 + vy_1^2 + vz_1^2) - 4*(gamm + 1)*(vx_0*vx_1 + vy_0*vy_1 + vz_0*vz_1) - 2*(2*bet - 1)*U_1 - (z_0 - z_1)*a_z1 - (y_0 - y_1)*a_y1 - (x_0 - x_1)*a_x1 - 3*((z_0 - z_1)*vz_1 + (y_0 - y_1)*vy_1 + (x_0 - x_1)*vx_1)^2/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2))*G*M_1/(((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2)*c^2) + (vx_0 - vx_2)*(2*(gamm + 1)*((z_0 - z_2)*vz_0 + (y_0 - y_2)*vy_0 + (x_0 - x_2)*vx_0) - (2*gamm + 1)*((z_0 - z_2)*vz_2 + (y_0 - y_2)*vy_2 + (x_0 - x_2)*vx_2))*G*M_2/(((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2)*c^2) + (vx_0 - vx_1)*(2*(gamm + 1)*((z_0 - z_1)*vz_0 + (y_0 - y_1)*vy_0 + (x_0 - x_1)*vx_0) - (2*gamm + 1)*((z_0 - z_1)*vz_1 + (y_0 - y_1)*vy_1 + (x_0 - x_1)*vx_1))*G*M_1/(((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2)*c^2) + 1/2*(4*(gamm + 1)*(((z_0 - z_1)*(vz_0 - vz_1) + (y_0 - y_1)*(vy_0 - vy_1) + (x_0 - x_1)*(vx_0 - vx_1))*G*M_1*vx_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2) - G*M_1*a_x1/sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) + ((z_0 - z_2)*(vz_0 - vz_2) + (y_0 - y_2)*(vy_0 - vy_2) + (x_0 - x_2)*(vx_0 - vx_2))*G*M_2*vx_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) - G*M_2*a_x2/sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)) - ((x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*(2*(2*gamm + 1)*(G*M_1/sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) + G*M_2/sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)) + vx_0^2 + vy_0^2 + vz_0^2) - 2*((2*gamm + 1)*(((z_0 - z_1)*(vz_0 - vz_1) + (y_0 - y_1)*(vy_0 - vy_1) + (x_0 - x_1)*(vx_0 - vx_1))*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2) + ((z_0 - z_2)*(vz_0 - vz_2) + (y_0 - y_2)*(vy_0 - vy_2) + (x_0 - x_2)*(vx_0 - vx_2))*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2)) + ((x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vx_0 + ((y_0 - y_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (y_0 - y_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vy_0 + ((z_0 - z_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (z_0 - z_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vz_0)*vx_0)/c^2 + 1/2*((2*gamm + 1)*((x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*(vx_0^2 + vy_0^2 + vz_0^2) - 2*(2*bet - 1)*(G*M_1/sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) + G*M_2/sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2))*((x_0 - x_2)*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2)) + (x_0 - x_2)*(2*(gamm + 1)*(vx_2^2 + vy_2^2 + vz_2^2) - 2*(2*bet - 1)*U_2 - (z_0 - z_2)*a_z2 - (y_0 - y_2)*a_y2 - (x_0 - x_2)*a_x2 - ((z_0 - z_2)*vz_2 + (y_0 - y_2)*vy_2 + (x_0 - x_2)*vx_2)^2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2))*G*M_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*(2*(gamm + 1)*(vx_1^2 + vy_1^2 + vz_1^2) - 2*(2*bet - 1)*U_1 - (z_0 - z_1)*a_z1 - (y_0 - y_1)*a_y1 - (x_0 - x_1)*a_x1 - ((z_0 - z_1)*vz_1 + (y_0 - y_1)*vy_1 + (x_0 - x_1)*vx_1)^2/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2))*G*M_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2) - (2*(x_0 - x_1)*((z_0 - z_1)*vz_1 + (y_0 - y_1)*vy_1 + (x_0 - x_1)*vx_1)^2/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^2 - 2*((z_0 - z_1)*vz_1 + (y_0 - y_1)*vy_1 + (x_0 - x_1)*vx_1)*vx_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) - a_x1)*G*M_1/sqrt((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2) - (2*(x_0 - x_2)*((z_0 - z_2)*vz_2 + (y_0 - y_2)*vy_2 + (x_0 - x_2)*vx_2)^2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^2 - 2*((z_0 - z_2)*vz_2 + (y_0 - y_2)*vy_2 + (x_0 - x_2)*vx_2)*vx_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2) - a_x2)*G*M_2/sqrt((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2) - 4*(gamm + 1)*(((x_0 - x_2)*G*M_2*vz_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1*vz_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vz_0 + ((x_0 - x_2)*G*M_2*vy_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1*vy_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vy_0 + ((x_0 - x_2)*G*M_2*vx_2/((z_0 - z_2)^2 + (y_0 - y_2)^2 + (x_0 - x_2)^2)^(3/2) + (x_0 - x_1)*G*M_1*vx_1/((z_0 - z_1)^2 + (y_0 - y_1)^2 + (x_0 - x_1)^2)^(3/2))*vx_0))/c^2;
/* expand(%); */
fullratsimp(%);