Пользуюсь Maple 10. 
Возникает ошибка при решении системы из 10 уравнений.
Думал, что ошибка возможна из-за синтаксиса) , поэтому всё набирал руками (а до этого копипастом).
Но всё равно та же проблема.
Возможно я допускаю тут какую-то глупую ошибку, но всё же...
![$
> restart;
>   a[11] := 0.176;>   a[22] := 0.011;>   a[33] := 35.06;>   a[12] := 0.07380000000;>   a[13] := 0;
>   a[23] := 0;>   a[14] := 2.926650000;>   a[24] := 1.323050000;>   a[34] := -61.35500000;>   a[44] := 155.53;
>   u := l l[x] + m m[x] + n n[x];>   v := l l[y] + m m[y] + n n[y];>   w := l l[z] + m m[z] + n n[z];
>   l[x] := \cos(\alpha) \cos(\beta);>   l[y] := \sin(\alpha) \sin(\omega) - \cos(\alpha) \sin(\beta) \cos(\omega);
>   l[z] := \sin(\alpha) \cos(\omega) - \sin(\alpha) \sin(\beta) \sin(\omega);
>   m[x] := \sin(\beta);>   m[y] := \cos(\beta) \cos(\omega);>   m[z] := -\cos(\beta) \sin(\omega);
>   n[x] := -\sin(\alpha) \cos(\beta);>   n[y] := \sin(\alpha) \sin(\beta) \cos(\omega) + \cos(\alpha) \sin(\omega);
>   n[z] := \cos(\alpha) \cos(\omega) - \sin(\alpha) \sin(\beta) \sin(\omega);$ $
> restart;
>   a[11] := 0.176;>   a[22] := 0.011;>   a[33] := 35.06;>   a[12] := 0.07380000000;>   a[13] := 0;
>   a[23] := 0;>   a[14] := 2.926650000;>   a[24] := 1.323050000;>   a[34] := -61.35500000;>   a[44] := 155.53;
>   u := l l[x] + m m[x] + n n[x];>   v := l l[y] + m m[y] + n n[y];>   w := l l[z] + m m[z] + n n[z];
>   l[x] := \cos(\alpha) \cos(\beta);>   l[y] := \sin(\alpha) \sin(\omega) - \cos(\alpha) \sin(\beta) \cos(\omega);
>   l[z] := \sin(\alpha) \cos(\omega) - \sin(\alpha) \sin(\beta) \sin(\omega);
>   m[x] := \sin(\beta);>   m[y] := \cos(\beta) \cos(\omega);>   m[z] := -\cos(\beta) \sin(\omega);
>   n[x] := -\sin(\alpha) \cos(\beta);>   n[y] := \sin(\alpha) \sin(\beta) \cos(\omega) + \cos(\alpha) \sin(\omega);
>   n[z] := \cos(\alpha) \cos(\omega) - \sin(\alpha) \sin(\beta) \sin(\omega);$](https://dxdy-02.korotkov.co.uk/f/1/9/7/197d64b06bd6cc5eb9fd82ab27cf335b82.png)
>  [/math]$ fsolve({ a[11] = (1 - 2 epsilon) b^2 c^2  l[x]^2  + a^2 c^2  l[y]^2  + a^2  b^2  l[z]^2 , 
>   
![$a[22] = (1 - 2 epsilon) b^2  c^2  m[x]^2  + a^2  c^2  m[y]^2  + a^2  b^2  m[z]^2 , 
>   a[33] = (1 - 2 epsilon) b^2  c^2  n[x]^2  + a^2  c^2  n[y]^2  + a^2  b^2  n[z]^2 , 
>  a[12] = (1 - 2 epsilon) b^2  c^2  l[x] m[x] + a^2  c^2  l[y] m[y] + a^2  b^2  l[z] m[z], 
>    a[13] = (1 - 2 epsilon) b^2  c^2  l[x] n[x] + a^2  c^2  l[y] n[y] + a^2  b^2  l[z] n[z], 
>   a[23] = (1 - 2 epsilon) b^2  c^2  n[x] m[x] + a^2  c^2  n[y] m[y] + a^2  b^2  n[z] m[z], 
>   a[14] = -(1 - 2 epsilon)b^2  c^2  l[x] u + a^2  c ^2 l[y] v + a^2  b^2  l[z] w, $ $a[22] = (1 - 2 epsilon) b^2  c^2  m[x]^2  + a^2  c^2  m[y]^2  + a^2  b^2  m[z]^2 , 
>   a[33] = (1 - 2 epsilon) b^2  c^2  n[x]^2  + a^2  c^2  n[y]^2  + a^2  b^2  n[z]^2 , 
>  a[12] = (1 - 2 epsilon) b^2  c^2  l[x] m[x] + a^2  c^2  l[y] m[y] + a^2  b^2  l[z] m[z], 
>    a[13] = (1 - 2 epsilon) b^2  c^2  l[x] n[x] + a^2  c^2  l[y] n[y] + a^2  b^2  l[z] n[z], 
>   a[23] = (1 - 2 epsilon) b^2  c^2  n[x] m[x] + a^2  c^2  n[y] m[y] + a^2  b^2  n[z] m[z], 
>   a[14] = -(1 - 2 epsilon)b^2  c^2  l[x] u + a^2  c ^2 l[y] v + a^2  b^2  l[z] w, $](https://dxdy-04.korotkov.co.uk/f/f/d/9/fd90d5fdb3d97a1e640e1369ee69755b82.png)
>    a[24] = -(1 - 2 epsilon) b^2  c^2  m[x] u + a^2  c ^2 m[y] v + a^2  b^2  m[z] w, 
>    a[34] = -(1 - 2 epsilon) b^2 c^2  n[x] u + a^2  c ^2 n[y] v + a^2 b^2  n[z] w,
>   a[44] = (1 - 2 epsilon) b^2  c^2  u^2  + a^2  c^2  v^2  + a^2  b^2  w^2  + (epsilon - 1) a^2  b^2  c^2  }, 
>   complex);$[/math] 
Error, (in fsolve) {alpha, beta, omega, a, l, m, n, xi, b, c} are in the equation, and are not solved for