Solve wont return obvious solution

Question

Ouvrir dans MATLAB Online

0 votes

I am trying to solve a symbolic system of nonlinear equations using

solve()

but it will not return obvious solution.

I try to run

sol = solve(K*x0)

with K and x0 being as follows:

x0 =
 x1
 y1
 z1
 x2
 y2
 z2
 x3
 y3
 z3
K =
[              21216.21*(x1 - 1.0*x2)^2,  21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2),  21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2),                                                   -21216.21*(x1 - 1.0*x2)^2,                                       -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2),                                       -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2),                                     0,                                     0,                                     0]
[  21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2),              21216.21*(y1 - 1.0*y2)^2,  21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2),                                       -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2),                                                   -21216.21*(y1 - 1.0*y2)^2,                                       -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2),                                     0,                                     0,                                     0]
[  21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2),  21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2),              21216.21*(z1 - 1.0*z2)^2,                                       -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2),                                       -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2),                                                   -21216.21*(z1 - 1.0*z2)^2,                                     0,                                     0,                                     0]
[             -21216.21*(x1 - 1.0*x2)^2, -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2), -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2),                         21216.21*(x1 - 1.0*x2)^2 + 21216.21*(x2 - 1.0*x3)^2, 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2) + 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2) + 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3),             -21216.21*(x2 - 1.0*x3)^2, -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3), -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3)]
[ -21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2),             -21216.21*(y1 - 1.0*y2)^2, -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2), 21216.21*(x1 - 1.0*x2)*(y1 - 1.0*y2) + 21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),                         21216.21*(y1 - 1.0*y2)^2 + 21216.21*(y2 - 1.0*y3)^2, 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2) + 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3), -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),             -21216.21*(y2 - 1.0*y3)^2, -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3)]
[ -21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2), -21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2),             -21216.21*(z1 - 1.0*z2)^2, 21216.21*(x1 - 1.0*x2)*(z1 - 1.0*z2) + 21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), 21216.21*(y1 - 1.0*y2)*(z1 - 1.0*z2) + 21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3),                         21216.21*(z1 - 1.0*z2)^2 + 21216.21*(z2 - 1.0*z3)^2, -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3), -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3),             -21216.21*(z2 - 1.0*z3)^2]
[                                     0,                                     0,                                     0,                                                   -21216.21*(x2 - 1.0*x3)^2,                                       -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),                                       -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3),              21216.21*(x2 - 1.0*x3)^2,  21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),  21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3)]
[                                     0,                                     0,                                     0,                                       -21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),                                                   -21216.21*(y2 - 1.0*y3)^2,                                       -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3),  21216.21*(x2 - 1.0*x3)*(y2 - 1.0*y3),              21216.21*(y2 - 1.0*y3)^2,  21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3)]
[                                     0,                                     0,                                     0,                                       -21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3),                                       -21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3),                                                   -21216.21*(z2 - 1.0*z3)^2,  21216.21*(x2 - 1.0*x3)*(z2 - 1.0*z3),  21216.21*(y2 - 1.0*y3)*(z2 - 1.0*z3),              21216.21*(z2 - 1.0*z3)^2]