Пользуюсь 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