Большая СЛАУ в символьном виде. Mathematica.

[Андрей123]

Здравствуйте. Может быть, кто-нибудь подскажет, как лучше решить СЛАУ с такой матрицей в Mathematica:

$$ \begin{array}{lllllllllllllll} -i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & -i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & -\frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & \frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ -i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & -i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & \frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & -\frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ -i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & \frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w1} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -\frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & 0 & 0 & 0 \\ -i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & -\frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w2} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & \frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & 0 & 0 & 0 \\ \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 & 0 & 0 & 0 \\ -i \text{w1} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 \\ -i \text{w2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 \\ \alpha (1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (1,1) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,1) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & -i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w1} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & \frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -i \text{c66}(3) \text{w1} \text{k3}(1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(3) \text{w1} \text{k3}(2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(3) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,3)} \\ 0 & 0 & 0 & 0 & 0 & 0 & -i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w2} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -\frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -i \text{c66}(3) \text{w2} \text{k3}(1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(3) \text{w2} \text{k3}(2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(3) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,3)} \\ 0 & 0 & 0 & 0 & 0 & 0 & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & \text{c66}(3) \text{k2}(1,3) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(3) \text{k2}(2,3) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & -i \text{w1} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}} \\ 0 & 0 & 0 & 0 & 0 & 0 & -i \text{w2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}} \\ 0 & 0 & 0 & 0 & 0 & 0 & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 & \alpha (1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 \end{array}$$

[lel0lel]

А в чём проблема?? Если вам в символьном виде решить нужно, то здесь действительно математика задумается, возможно и навсегда)). Но с численным решением никаких проблем не будет, решали с матрицами в разы побольше. Решение можно оформить как функцию ваших параметров и тогда символьное решение ни к чему совсем. Поэтому, если уж здесь Вам было не лень ентого крокодила набирать, то теперь тоже в математике, ну и функция LinearSolve[m,b]. Посмотрите в хэлпе, у неё есть опция Method, она может быть Вам полезна.

[Андрей123]

Решение нужно в символьном виде. После некоторых доработок размерность системы стала меньше. И ее в явном виде удается решить методом Краммера, а LinearSolve, почему-то, надолго задумывается и выдает сообщение, что не хватает памяти. Я думал, что, может быть, есть что-то эффективнее.

P.S. В Mathematica есть функция TeXForm, поэтому здесь систему мне набирать не пришлось.

[lel0lel]

Андрей123 в сообщении #317582 писал(а):

Решение нужно в символьном виде.

Честно признаться, сложно представить для чего. Вы бы могли написать функцию которая бы возвращала столбец решений, а её аргументами были бы параметры матрицы. Но воля ваша, главное обязательно проверьте совпадение численного решения при определенных параметрах с решением аналитическим при тех же параметрах, так, на всякий случай.

Андрей123 в сообщении #317582 писал(а):

P.S. В Mathematica есть функция TeXForm, поэтому здесь систему мне набирать не пришлось.

Понятно)

[Андрей123]

Комбинацию решений с некоторыми функциями надо потом обращать по Фурье. С численным решением, кстати, тоже не все хорошо. При некоторых заначениях параметров пишет, что матрица плохо обсловлена.