Hi
which release do you use? The following code sets the 6 equations == 0 and produces a result on R2018b:
MAx = -0.2376;
MAy = 0.8430;
MAz = -0.7356;
MBx = -0.2376;
MBy = 0.6675;
MBz = -0.3846;
syms r_ax r_ay r_az Fx Fy Fz;
r_bx = r_ax +0.0351;
r_by = r_ay;
r_bz = r_az;
eqn1 = ((r_ay*Fz) - (r_az*Fy))- MAx == 0;
eqn2 = ((r_az*Fx) - (r_ax*Fz))- MAy == 0;
eqn3 = ((r_ax*Fy) - (r_ay*Fx))- MAz == 0;
eqn4 = ((r_by*Fz) - (r_bz*Fy))- MBx == 0;
eqn5 = ((r_bz*Fx) - (r_bx*Fz))- MBy == 0;
eqn6 = ((r_bx*Fy) - (r_by*Fx))- MBz == 0;
Eq = [eqn1, eqn2, eqn3, eqn4, eqn5, eqn6];
vars = [r_ax, r_ay, r_az, Fx, Fy, Fz];
sol = solve(Eq,vars)
A=double([sol.r_ax; sol.r_ay; sol.r_az; sol.Fx; sol.Fy; sol.Fz])
The values for all M are like the ones you commented out. The calculated r-values differ from your ones.
Best regards
Stephan
Best Answer